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CHAPTER  ONE 
INTRODUCTION 


How  best  to  detect  known  signals  in  noise  is  the  object  of 
much  study  and  effort.  Among  the  forms  in  which  this  problem  occurs  is 
that  of  deciding  when  a signal  is  present  and  from  what  direction  it  is 
arriving.  In  this  study,  we  are  concerned  with  the  capabilities  of 
sonar  arrays  to  perform  these  functions. 

1.1  BACKGROUND  OF  THE  PROBLEM 

The  emphasis  to  be  found  here  is  upon  analytical  rigor,  and 
so  the  class  of  array  models  treated  is  restricted.  However,  in  the 
presentation  of  this  work  we  have  striven  to  facilitate  extensions  to 
related  problems  by  being  quite  general  in  the  earlier  stages  of  the 


analysis. 


Specifically,  we  shall  derive  the  detection  and  bearing  esti- 


mation performance  of  a "multiplicative  linear  array"  in  comparison 
with  that  of  the  "standard  linear  array." 

1.1.1  LINEAR  ARRAYS 

By  the  term  "linear  array"  we  refer  to  a set  of  2M  isotropic 
point  receivers  (hydrophones),  arranged  in  a straight  line  with  ehual 
spacing  d.  In  the  literature,  this  configuration  is  often  called  a 
"broadside  array,"  particularly  in  connection  with  transmitting  systems. 
In  receiving  applications  this  configuration  makes  use  of  the  greater 
coherence  or  correlation  among  the  samples  of  the  signal  present  at  the 
hydrophones  over  those  of  the  noise,  when  this  condition  exists.  In 
effect,  it  is  a method  of  simultaneously  obtaining  2M  samples  or  obser- 
vations of  the  signal  plus  noise  process,  a spatial  analog  to  time 
sampling.  H* 


When  these  samples  are  summed,  the  signal  components  tend  to 
reinforce  each  other  while  the  noise  components,  under  certain  conditions, 
tend  to  cancel  or  "average  out."  Also,  when  the  signal  is  in  the  form 
of  a sinusoidal  plane  wave,  the  individual  signal  components  sum  in  such 
a way  as  to  give  the  array  a highly  directive  receiving  pattern. 

In  this  plane  wave  case,  which  shall  be  our  signal  model 
throughout  the  study,  the  signal  components  at  the  hydrophones  can  be 
modeled  as  sk(t)  = Ak(t)cos[u>t-0-(k-l )esinct],  k = 1,2,...2M,  where  a is  the 
bearing  of  the  (distant)  source  relative  to  the  perpendicular  of  the 
array,  0 is  an  arbitrary  initial  phase,  and  6 = 2Tid/x.  The  angular 
frequency  is  equal  to  2ttc/X:  where  c and  \ are  the  velocity  of  sound 
and  the  wavelength  in  water,  respectively. 

For  a symmetrical  gain  array  (Ak  = A2M+1_k,  k = 1,2,...M), 
the  sum  of  the  signal  components  is 


2M  M 

s(t)  ■ ^ sk(t)  = 2cos[ut-0-|(2M-l)0sina]^  AM+1_kcos[-^(2k-l  )0sina] 
k=l  k=l 

(1-1) 


and  for  a uniform  gain  array  (Ak  = A,  all  k),  the  sum  is 

• M 

s(t)  * 2Acos[u)t-0-i(2M-l  )6sina]^  cos[i(2k-l  )0s1na] 

k=l 


(1-2) 


* 2MAD(iBsina;  2M)cos[ut-0-^(2M-'l)Bsina], 


where  D(x ; N)  = sinNx/Nsinx. 


WWW W». 


Tne  sum  of  the  noise  components  meanwhile,  to  illustrate  with 

the  simple  case  in  which  the  components  are  independent,  zero-mean,  and 

2 2 
have  equal  variances  a , would  have  a variance  equal  to  2Mo  . Thus,  in 

this  ideal  case,  at  zero  bearing  the  uniform  linear  array  sum  would 

achieve  a signal  to  noise  ratio  "gain"  of 

<3,«W/(SNR>in  * + A •*  2M  (1-3) 

. 1 ,n  2Ho  2c 

In  general  there  exist  proportional  relationships  between  the 
number  of  elements  and  the  array  processing  gain,  and  between  the  number 
of  elements  and  the  array  directivity,  considered  separately.  However, 
as  will  be  discussed  below,  the  relationship  between  directivity  and  SNR 
gain  is  subject  to  other  considerations. 

Schelkunoff  [1],  Dolph  [2],  and  others,  by  considering  the 
directivity  function  D(j  Bsin  <r,  N)  as  a polynomial  in  cos  (5  B sin  a)> 
have  shown  that  the  directivity  of  a linear  array  whose  elements  are 
summed  can  be  optimized  as  a function  of  the  inter-element  spacings 
and/or  the  individual  weightings  of  the  elements.  In  this  study  we 
shall  restrict  our  attention  to  uniform  gain  linear  arrays  with  equal 
spacings. 

1.1.2  SQUARE-LAW  ARRAY  PROCESSING 

In  his  textbook  an  acoustic  signal  processing,  Horton  [3] 
terms  "standard  detector"  that  configuration  in  which  the  sum  of  the 
hydrophone  outputs  is  squared  and  averaged  by  means  of  a lowpass  filter, 
usually  with  a variable  time  constant  or  bandwidth  to  accommodate 
tracking  conditions.  The  output  of  the  filter  is  then  observed  as  the 
array  is  steered  either  physically  or  electronically.  When  the 
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observation  exceeds  a certain  threshold,  the  decision,  "signal  present," 
is  made.  When,  as  the  steering  takes  place,  the  observation  passes 
through  a peak  value,  then  the  bearing  of  the  signal  source  is  taken 
to  be  the  direction  in  which  the  array  was  steered  at  the  time.  A 
diagram  representing  the  standard  or  square-law  processor  is  shown  in 
Figure  1-1. 

1.1.3  MULTIPLICATIVE  ARRAY  PROCESSING 

Berman  and  Clay  [4]  showed  that  methods  of  generating  poly- 
nomials can  be  found  which  require  fewer  array  elements  than  the  uniform 
linear  array  whose  outputs  are  summed.  One  such  method  was  shown  to 
be  taking  time  averages  of  products  of  various  array  element  signals, 
rather  than  summing  them,  and  manipulating  the  relative  spacings  among 
elements.  Thus  a 4-element  TAP  (time-averaged-product)  array  achieved 
the  direction  pattern  of  an  8-element  additive  array,  and  an  11-element 
TAP  array  matched  the  pattern  of  a 1024-element  additive  array. 

Tucker  and  Welsby  [5-8]  and  Cath  [9]  considered  systems  in 
which  the  N elements  of  a uniform  gain  linear  array  are  divided  into  two 
groups;  the  summed  output  of  one  group  is  multiplied  by  that  of  the  other, 
and  the  product  is  smoothed  by  lowpass  filtering--in  essence,  the  corre- 
lation of  two  arrays  of  smaller  size.  The  beamwidth  of  this  configuration, 
shown  diagrammatical ly  in  Figure  1-2,  was  snown  to  be  about  one  half  that 
of  the  corresponding  N-element  standard  array.  Also,  the  relative 
magnitudes  of  the  unwanted  sideiobes  are  greatly  reduced,  the  largest 
sidelobes  being  of  negative  polarity  and  thus  removable  by  rectification. 
These  advantages  are  illustrated  in  Figure  1-3  for  a specific  case. 
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Steering 


FIGURE  1-2  MULTIPLICATIVE  ARRAY  PROCESSOR 


NORMALIZED  DIRECTIVITY  FUNCTION 


FIGURE  1-3.  ARRAY  DIRECTIVITY  FUNCTIONS 


While  multiplier  configurations  can  offer  more  efficient  beam- 
forming,  in  the  processing  area  they  represent  a departure  from  conven- 
tional detection  toward  correlation  techniques.  Thus  there  come  into  play 
additional  tradeoffs,  earlier  stated  by  Faran  and  Hills  [10],  and  subse- 
quently confirmed  by  these  array  designers  in  their  calculations  and 
experiments. 

It  was  found  that  any  improvement  in  the  directional  pattern 
over  the  standard  array  by  multiplicative  processing  is  accompanied  by 
degradation  of  the  noise  factor,  or  decreased  SNR  gain.  Also,  the  multi- 
plication tends  to  mar  the  performance  of  discrimination  between  two 
targets  or  sources  of  unequal  strength,  while  for  equal  strengths  the 
performance  was  demonstrated  to  be  greater. 

The  degradation  in  array  gain  for  TAP  arrays  relative  to  stan- 
dard was  examined  by  Fakley  [11]  for  narrowband  noise  and  coherent 
narrowband  signals.  He  shows  that  in  a specific  instance,  for  the  same 
probability  of  correct  detection,  TAP  operations  require  a higher  input 
SNR.  Fakley  also  shows  the  occurrence  of  the  small  signal  suppression 
effect  in  the  case  of  two  sources,  as  well  as  the  interesting  case  of  a 
line  source.  In  the  latter  he  presents  a calculation  implying  that  the 
multiplicative  process  output  will  not  track  the  angular  energy  distri- 
bution of  the  extended  source,  but  rather  peaks  at  its  center  bearing. 

Thus,  among  the  experimenters  and  researchers  in  the  field  of 
sonar  arrays,  there  exists  a running  debate  on  the  relative  merits  and 
tradeoffs  between  standard  and  multiplicative  array  processing. 
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DESCRIPTION  OF  THE  PRESENT  WORK 


This  "debate"  over  the  relative  merits  of  standard  and  multipli- 


cative arrays  has  been  handicapped  by  a significant  consideration.  Tucker  [5 j 


alluded  to  it  as  "...  the  difficult  question  of  whether  the  better 


signal/noise  performance  ...  is  really  indicative  of  a higher  probability 


of  detection."  That  is  to  say,  the  present  methods  of  characterizing  the 


array  processing  performance  are  inadequate.  • Except  in  a few  instances 


the  performance  criteria,  because  of  the  analytical  difficulty  of  obtaining 


statistical  detection  parameters,  have  been  SNR  gains,  calculated  usually 


for  rather  specialized  cases  such  as  white  noise. 


1.2.1 


OBJECTIVES  > ND  PROCEDURES 


The  object  of  this  study  is  to  make  a positive  contribution  to 


the  array  processing  debate  by  deriving  exact  probability  density  functions 


for  the  filter  outputs  of  reasonably  general  narrowband  models  of  both  the 


standard  array  and  the  multiplicative,  with  the  major  emphasis  on  the 


latter.  In  doing  so,  we  have  been  able  to  develop  expressions  for  detec- 


tion performance  criteria  in  which  appear  all  the  various  parameters  which 


go  into  the  tradeoffs.  This  material  is  presented  in  Chapters  Two  and 


Three. 


Also,  we  have  utilized  these  probability  distributions  to 


calculate  the  theoretical  best  bearing  estimation  performance  (minimum 


variance,  error)  attainable  due  to  the  information  contained  in  the 


array  filter  output  signal.  Forms  for  maximum  likelihood  estimators 


are  derived  and  evaluated  with  respect  to  the  predicted  minimum  variance. 


This  work  is  found  in  Chapters  Four  and  Five. 
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In  Chapter  Six  the  study  is  concluded  by  summarizing  the 


performance  of  the  multiplicative  array  model,  particularly  in  comparison 


wth  the  standard  model.  Suggestions  also  are  given  for  extensions  of 


the  work. 


1.2.2  THE  ARRAY  CONFIGURATIONS  TO  BE  STUDIED 


The  standard  array  configuration,  as  shown  in  Figure  1-1,  and 


the  multiplicative  array  configuration  mentioned  by  Welsby  and  Tucker, 


given  in  Figure  1-2,  are  the  models  to  be  studied— primarily  the  multi- 


plicative one.  That  is,  the  multiplicative  array  studied  shall  be 


considered  to  have  a total  of  2M  isotropic  point  receivers  or  hydro- 


phones equally  spaced  d units  along  a straight  line  and  each  with  unity 


gain.  The  outputs  of  these  hydrophones,  indexed  sequentially  from  one 


end  of  the  array  to  the  other,  are  denoted  by  xk(t),  k = 1,2,...,2M. 
The  standard  array  shall  be  considered  to  have  N hydrophones. 


In  the  case  of  the  standard  array,  the  output  of  the  processor 


is  the  filtered  version  of 

r n i2 


(t)  = nT  xk(t)  = [NA(t)D(£6sina;  N)  cos(wt-0-£(N-l )Gsina)  + n(t)F 


(1-4) 


For  the  multiplicative  array,  the  output  is,  before  filtering, 


H IT  2M 

I xk(t>  I xk(t) 

k=l  J Lk=M+l 


i!4A { t ) D { i 6s i na ; M)cos[ut-0-r(M-l )6sina]  + n,(t)} 


X (MA(t)D(-ksina;  M)cos[wt-P-i(3M-l )psina]  + n9(t)}.  (1-5) 
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The  form  of  the  signal  used  here  and  throughout  the  study  is  that  of  a 

known  narrowband  source  whose  bandwidth  is  centered  about  o>/2tt,  and 

located  at  bearing  sufficiently  far  away  that  the  array  experiences 

It  as  a plane  wave.  This  form  is  emphasized  diagrammatical ly  in 

Figures  1-1  and  1-2  by  the  presence  of  narrowband  filters.  For  the  sake 

of  avoiding  ambiguities,  our  attention  will  be  mainly  on  bearings  which 

fall  within  the  major  lobe  of  the  array.  That  is,  |Ng  sin  ot | <tt  for  the 

standard  array  and  |Mg  sin  a|<ir  for  the  multiplicative.  Any  steering 

done  by  the  array  shall  be  assumed  ideal,  that  is  unity  gain,  distor- 

* 

tionless,  time  delays. 

1.2.3  THE  NOISE  PROCESS  MODEL 


Although  the  conditions  do  not  always  warrant  it,  for  conveni- 
ence the  noise  process  in  the  locality  of  the  array  shall  be  considered 
stationary  Gaussian  with  zero  mean.  In  Appendix  A a discussion  is 
presented  to  motivate  the  selection  of  values  for  the  elements  of  the 
following  covariance  matrices  used  to  characterize  the  noise  components 
of  the  array  summations: 

SQUARE-LAW  ARRAY:  n(t)  = nc ( t ) cos  wt  + ns ( t ) sin  u)t 


Cov[n(.(t1),  n _ ( t-j ) ] = 


a ' 0 

0 a1 

That  is,  adopting  the  usual  quadrature  expansion  of  narrowband  Gaussian 
processes,  the  sine  and  cosine  components  of  the  noise  are  onsidered 
independent  at  the  same  time  instant. 


(1-6) 


We  shall  not  consider  steering  and  beamforming  as  such,  but  mention 
these  operations  only  to  indicate  that  the  models  we  employ  can 
accommodate  them,  and  to  support  our  practice  of  considering  only 
relative  bearings  (i.e.,  relative  to  the  steering  direction'  within 
the  major  lobe  of  the  array  pattern. 
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MULTIPLICATIVE  ARRAY:  n1(t)  * njc(t)  cos  ut  + n1$(t)  sin  ut  i = 1,2 

Cov[»1e(t1).nls(t1),i<2c(t1),n2s(t1)]  = 


a? 


pof2 


-™f2 


nfz  vfz 


-mlc2 


pa,  a, 


12 


4 


A K 


0-7) 


That  is,  the  noise  correlation  at  the  array  sums  is  considered  to  be  such 

that 


E(n1cn2c)  = E(nlsn2s>  = 0»1°2 


Elnlcn2s>  * -E("lsn2c>  = rola2‘ 


(1-8) 
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CHAPTER  TWO 


ARRAY  PROBABILITY  DENSITY  FUNCTIONS 

The  most  basic  representation  of  the  multiplicative  array 
model  we  are  considering  is  that  of  the  product  of  two  nonzero-mean 
Gaussian  random  variables.  This  description  encompasses  a large  number 
of  cases  which  are  of  interest,  which  cases  are  often  most  readily 
approached  individually  because  one  or  more  simplifying  assumptions  are 
traditionally  made. 

In  this  chapter  we  develop  the  probability  density  function 
for  the  generalized  narrowband  case  of  the  multiplicative  array  configu- 
ration, showinq  also  the  related  function  for  the  square-law  array. 
These,  in  turn,  are  specialized  to  the  cases  which  are  anticipated  to 
be  of  most  interest  and,  moreover,  are  relatively  amenable  to  the 
analysis  in  the  succeeding  chapters. 

The  situation  in  which  the  noise  is  symmetric  across  the 
array  and  has  a narrowband  spectrum  which  is  even  about  the  center 
frequency  is  selected  as  the  "main  case."  The  probability  density 
functions  for  the  multiplicative  and  square-law  array  filter  outputs 
are  applied  to  this  case,  yielding  forms  in  which  all  the  array  para- 
meters appear  explicitly. 

Computed  results  are  shown  for  the  probability  density 
functions,  displaying  the  effects  of  the  various  distribution  para- 


meters. 


The  mean  and  variance  of  the  filter  output  are  calculated,  and 
used  to  express  the  signal  to  noise  ratio  (SNR)  at  the  output  of  the 
array  processors. 


2.1  PROBABILITY  DENSITY  FUNCTION  FOR  THE  GENERALIZED  NARROWBAND 

MULTIPLICATIVE  ARRAY  PROCESSOR 


We  now  give  our  attention  to  a rather  general  case  in  which 
the  multiplicative  array  summations  of  Figure  1-2  are  modeled  as  two 
narrowband  processes,  u-j(t)  and  u2(t),  whose  product  z(t)  is  passed 
through  a zonal  lowpass  filter  to  yield  the  output  y(t),  as  illustrated 
in  Figure  2-1.  The  processes  u^t)  and  u2(t)  are  said  to  consist  of 
the  super-position  of  known  signals  s-j(t)  and  s^C t ) and  noise  processes 
n-j (t)  and  n2(t),  respectively.  At  the  same  instant  of  time,  n^(t)  and 
n2(t)  are  assumed  to  be  jointly  normal  random  variables  with  variances 


2 2 
- 1 and  o2  , 

1.2.3. 


correlation  coefficient  p,  and  zero  means,  as  in  Section 


The  multiplier  is  assumed  to  be  instantaneous,  so  that  its 
output  can  be  written 

z(t)  = Uj(t)  x u2(t) 

= Cs1  (t)  + n-j (t)]  x [s2(t)  + n?'t)] 

■ { lu,(t)  + U2(t)P  - l[u,(t)  - U2(t)]2 


= [s3(t)  + n3(t)J2  - [s4(t)  + n4(t)]2, 

(2-1) 

where 

s3>4(t)  = ^{s1 (t)  ± s2(t)] 

(2-la) 

and 

n3>4(t)  = ^fn^t)  ± n2(t)]. 

(2-lb) 

This  arrangement  may  be  recognized  as  the  old  "quarter-square  multiplier" 
idea,  used  in  analog  computation. 
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FIGURE  2-1  MATHEMATICAL  MODEL  OF  MULTIPLICATIVE  ARRAY 
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For  the  new  variables  we  have  defined,  we  have  the  following 


moments : 


E(n3)  = 0 
E{n4)  = 0 

E (n| ) - = ^ (c^  + 2po-|02  + o|) 

E(n|)  = o|  = ^ (a^1  - 2po-|02  + op  ' 

al). 


(2- 2a) 


E(n,n,)  8 Ra3c4  - i (of 


S S' 


(2-2b) 


Note  that  if  o-j  = o2»  then  R = 0;  that  is,  and  n^  are  uncorrelated. 

Specialization  to  what  we  are  calling  the  narrowband  case 
consists  in  the  conventional  assumption  that  the  following  quadrature 
representations  apply: 

s3(t)  = S3(t)  cos[ut-e3(t)],  s4(t)  = s4(t)  cos[iot-e4(t)] 

n3 ( t ) = n3c(t)  cos  wt  + n3s(t)  sin  ut  = N3(t)  cos[ut-*3(t)] 

n4(t)  = n4c(t)  cos  cat  + n4s(t)  sin  wt  = N4(t)  cos[ut-*4(t)],  (2-3) 

where  if  Sj(t)  = S.(t)  cos[ut-6j(t)] 

1 = 1,2 

ni  (t)  = ( t ) cosTiot-^^  ( t)] 


(2-3a) 


then 


^3,4  " 4 ^1+^2±^^1^2  cos(0i"02^ 


tane 


sino1±S2  sine2 


3,4  S1  coso^+C^  cose2 


N3,4  = \ [Ni+N2±2N1N2  cos(^-^2^ 
sin^iN^  sin^ 


tan$ 


3,4  cos<j>1±N2  cos$2 
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(2-3b) 


-*•  
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With  this  narrowband  representation,  we  have  for  the  output 
of  the  multiplier, 

* 

2 - (s3  + n3)2  - (s4  + n4)2 

c ^ tS3+2S3N3  cos(<^3-g3)  + Nj  - Sj  - 2S4N  cos^-e^-Np 

+ (terms  with  frequency  2u).  (2-4) 

What  is  meant  by  "zonal  lowpass  filter"  is  that  the  filter  output  y is 
equal  to  the  multiplier  output  z,  less  the  terms  of  frequency  2u>.  Thus 
if  we  define  the  sum  and  difference  terms 

*1  (t)  = s3(t)  + n3(t)  = j [u-,(t)  + u2(t)]  = 2,(0  cosfj)t-^,(t)] 

a z,c(t)  cos  '•  + z,s(t)  sin  ut  (2-4a) 

z2(t)  = s4(t)  + n4(t)  = J [u, (t)  - u2(t)]  = Z2(t)  cos[ut-*2(t)] 

= z2c(t)  cos  ojt  + z (t)  sin  ojt  (2-4b) 

we  find  that 

y(t)  * j [zf(t)  - z|(t)]  (2-5) 

where 

l\  2 = } [Sj  + S2  + N2  + W\  ± 2S,S2  005(02-6,) 

+ 2S,N,  cos (4.,-e-, ) ± 2S,N2  cos($2-e,) 

± 2S2N,  cosU,-e2)  + 2S2N2  cos(<j»2-e2)  ± 2N,N2  cos(*,-*2)] 

= S3,4  + N3,4  + 2S3,4N3,4  cos ^3,4"03,4^  ^2“5a^ 
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and 


Sj  sin6j  ± Sp  sinQg  + 
tan^i *2  - Sj  cose^  ± S2  cos82  + N^ 


si n4>i  ± Ng  sinijig 
cos<j>1  ± N2  ccs<f>2 


_ S3,4  s1n03,4  * N3.4  S1‘n»3.4 
S3,4  C(S03,4  + N3,4  COS*3,4 


(2-5b) 


Our  object  is  to  find  the  probability  density  function  of  the  filter 
output  y(t)  at  a given  instant  in  time.  (From  this  point  on,  reference 
to  time  will  be  suppressed.) 

2.1.1  DISTRIBUTION  OF  THE  SUM  AND  DIFFERENCE  TERMS 

If,  to  use  vector  notation,  we  refer  to  the  input  noise  com- 
ponents as  x'  = (n1c>  n^s>  n2c,  n2s),  where  the  prime  (')  is  used  to 
indicate  the  transpose,  then  we  may  express  the  joint  probability  den- 
sity function  of  the  input  noise  components  as 


~~~  e*pL- 

47iVdetK 

A 


where  Kx,  the  covariance  matrix,  is  postulated  to  be 


Kx  ■ Cov[nk(t),  n,s(t),  n2c(t),  n2s(t)] 
= E[x/x] 

? 0 


P0-.0 


■ 2 


rcla2 


2 

°1 


',01°2 


po-,0 


12 


pU'j  1^2 


rala2 


°2 

0 


ro^a2 


• pala2 
0 

2 
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That  is,  the  noise  correlation  at  the  inputs  to  the  multiplier  is  con 
sidered  to  be  such  that,  for  the  same  time  instant, 


E(nlc"2c>  - E("ls"2s>  ' °°\°2 
E("lcn2s>  ‘ -E<n1sn2c)  = ral°2* 


We  may  consider  the  sum  and  difference  noise  components 
- = ^n3c*  n3s’  n4c’  n4s^  as  a ^near  transofrmation  n_  - hx,  where 
the  matrix  A is  given  by 


1/2 

0 

0 

1/2 

1/2 

0 

0 

1/2 

1/2 

0 

0 

1/2 

-1/2 

0 

0 

-1/2 

We  may  then  write  the  probability  density  function  of  the  variables  n 
as 

P](n)  • 4^0  exp[-  \ n'icVl, 

where  the  new  covariance  matrix  «n  is  easily  computed  to  be 
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I 

I 


I 

I 


I 


0 ^°3a4 


= AV 


0 

Ro3°4 

“R'°3a4 


2 o i 

°3  ‘K  a3a4 

f)  i 2 

-K  a3a4  a4 

R a3a^  0 


R03O4 


0 


(2-7) 


and  D = /detKn  = a|a|[i  _ r2  _ (R1)2].  (2-7a) 

The  variances  o2  and  a2,  and  R are  as  given  in  (2-2),  and  R'o^  ~ ro^o^. 

It  follows  then  that  the  quadrature  expansions  of  the  signal 
plus  noise  variables  z1  and  z 2 as  defined  above  have  the  density  function 


P2(z)c.z,s.z2c,z2s)  * 4^5  expt-  \ (z-s.) • K"1  (z-s.) ] (2-8) 

whc-e  (z-s.)'  = (zi c“s3  cos03’zis"S3  sine3’z2c"S4  cose4’z2s'S4  sin04^  (2“9) 


2.1.2  DISTRIBUTION  OF  THE  FILTER  OUTPUT,  y 

Recall  that  the  filter  output  is  y =^-  (Z^2  - Z|).  The  density 
function  of  Z1  and  Z?  is 


•2tt  —2tt 


ft  it 

P4 ( ^ 1 , Z 2 ) = J d<{>i  J ^^2  P3 (^i  ,!^i  ’^2*^2^ 


ft.  TI  fC.T\ 

= Z^Z^  J I dtf>2  P2(zi  cos4>1 , Z^  si nt^1 , Z 2 cos <{>2,Z2  sinf2) 

7Q  •'0 

(2-1 Oa) 
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j 


t 


is 


< 


5 - ; 


Z,Z*  /*2n 


“2^  d<t>2  ®xpC"  2^— ^ ' *Si  ^ (2-1  Ob) 


4ir 


where  (Z  - s_) ' = (Z^  cos^  - S3  ccse3,  Z1  si - S3  sine3, 

Z2  cos4>2  - S4  cose4,  Z2  sin^2  - S4  sine4).  (2-11) 
The  quadratic  form  appearing  in  the  integral's  exponent  reduces  to 


Q = 


. -1 


§0  [o|(Z{  + Sp  + 02(Z2  + «J)  . 2a3o4XS3S4  cos(e4-e3-x) 
-2O3O4XZ1Z2  COS (♦, -♦,-*) 

-2o4VZ^  cos(<^-v)  - 2o3WZ2  cos(<j>2-w)]  (2-12) 


where  X2  = R2  + (R ' ) 2 , tanx  = R'/R 


V2  = aJS2  + a2X2S2  - 2a3a4XS3S4  cos(03-e4+x) 


W2  = cj2S4  + a|X2S3  - 2a3o4XS3S4  cos(q3-64+x) 


tanv  = 


o4S3  sino3  - a3XS4  sin(e4-x) 


o4S3  cose3  - o3XS4  cos(e4-x) 


o3S4  sin04  - o4XS3  sin(o3+x) 
^anW  o3S4  cos04  - o4XS3  cos(o3+x) 


(2-13) 


Using  the  relationship  e 


xcosw 


l cosnw 


, . 4''n  = 

n (2,n  > 


= 0 


and 


n=0  n n n (2,n  > 0 

performing  the  integration  with  respect  to  <t>1  and  4> 2 , in  a manner  quite 

similar  to  Middleton^16^  (chapter  9)  we  obtain 
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P4(2,.Z2)  ^expjj  [o’Zf  * „’z|  * 


X l CJ 

L _ m 


m=0 


m m 


(^)  >m  f$£)  >m  f^)  cosm(v-w-x)  (2-14, 


where  we  have  used  U2  = o2S|  + a2S|  - 2a3a4XS3S4  cos(e4-e3-x), 

and  I ( ) is  the  m:th  order  modified  Bessel  function  of  the  first  kind, 
m 

To  obtain  the  density  for  y,  we  use  the  transformations 


(2-15) 


and 


Zj  = /2y  coshu 
72  = /5y  sinhu 

Zj  = tty  sinhu 
Z2  = tty  cos  hu 


for 


Z^  > Z2,  or  y > 0 and  u > 0 


for  Z-|  < Z2,  or  y < 0 and  v > 0. 


(2-16) 


Then  we  find  that  the  density  for  y is  given  by  the  expression 

00  00 

P6(y)  = du  P5(y,u)  =j£  du  p4(»'2y  coshu  , /2y  sinhu  ),  y > 0 


■jf du  P4(/-2y  sinhu  , tty  coshu  ),  y < 0. 


= M 


exp 


-(4«4)b!-u2 


2D 


y c cos  m(v-w-x) 
L-  m 


m=0 


i 

1 

1 

i 


f du  sinh2u  e"a  cosh2u  lm(b  coshu  )Im(c  sinhu  )Im(d  sinh2u) 
JQ 

(2-17) 
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itftMMHMiikttaHitoiii 


with 


a = |y| (a|+CTj)/2D,  d = o3c4X|y|/D,  and 
(b,c)  = (o^fty/D,  ojW^y/D),  y > 0 

= (o3 M *^2y/D,  a4V  /-2y/D),  y < 0. 


(2-18) 


By  writing  x = 1 + cosh2u  t we  may  change  the  above  integral  to 
J e*£"  dx  e'axIn[^7)]  IB  [cJT]  Im  [d/5TFZ)] 


m 


1 _a  r r (b/2v^)<‘"  v2^;  r . -OA/  pUiTRTin  ^ 

2 e Jo  Jo  HIFfiT!  JQ  dx  e (x+2>  x 2 


(b/2/2)2ritm  (5d)2Hm 


1. 


-ax,..  ,„xn+k+m  k-t^m 


x vsf* 


1 .N2k+m 


la""  n+k+m  (b/2/2)2n+m  ^ ,n+k+mN  „n+k+m-r 
= 2e  lQki0  lQ  n!(n+m)!  kfrk+mJT  v r '2 


,1 


-ax  k^m+r 
dx  e axx  2 I 


m 


lc$l 


(2-19) 


Making  use  of  Gradshteyn  and  Ryzhik^2^  (formulas  6,643.4  and  9.220.2), 
we  obtain  for  (2-19)  the  following  expression: 


00  00 


/l.  \2n+m/l.x2k+m  ,,  „ 

n+k+m  (^b)  ty)  ,n+k+Iru2k"r 


(4}  2a  exp(a+  8a)j0  JQ  n!  (n+m ) ! k ! ( k+m ) ! ( r ^k+r  (k+r)! 


x i m (i£-) 

x Lk+rl  8a ' 


(2-20) 
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ii 


i 
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where  the  L^+r(  ) are  the  Laguerre  polynomials  of  order  k+r  and  parameter 
m.  Substituting  in  (2-17)  and  for  a,  b,  c,  and  d,  and  rearranging  terms 
for  clarity,  we  have  finally 


p6(y)  = sWexp 


<D  ao  cc  n+k+m 

till  l 

rfO  n=0  k=0  r=0 


( ( f2\n/g3°4X2  f/  D ( 

^2D2(o§+o||y  (p-/  \D(a§+a|y  ) \a|+o|  J 


, \ (k+r)!  /k+n+mA 

‘ VK swSisyr 1 r > 


x |y|kwn-rLin 


k+r 


(2D (o|+o]p) 


(2-21) 


with  (f  »9»Oq) 


= 


^a^v,  a3«,  a^;,  > > u 


(a3W,  a^V , a?),  y < 0. 


(2-22) 


A recent  paper  by  Andrews^  ^ gives  an  expression  for  the  den- 
sity function  in  nearly  the  same  situation,  by  means  of  the  character- 
istic function  method.  Comparing  the  development  here  with  that  of 
Andrews,  it  appears  that  following  the  strenuous  development  associated 
with  the  computations  of  the  characteristic  function,  one  still  needs 
to  solve  a convolution  integral  r^Ch  as  our  (2-19)  above.  The  different 
final  expression  he  obtains  results  simply  from  a different  method  of 


2- 


attacking  this  type  of  integral.  For  the  special  cases,  of  course2  the 
form  of  the  expressions  given  here  is  identical  to  that  of  Andrews. 
Moreover,  the  computed  results  we  show  below  match  his  wherever  appli- 
cable. 

2.1.3  REDUCTION  OF  THE  GENERAL  EXPRESSION  TO  SPECIAL  CASES 


For  practical  applications  and  for  purposes  of  checking  our 
results  with  those  of  previous  authors,  it  is  instructive  to  reduce  the 
general  expression  of  (2-21)  under  certain  assumptions. 

2. 1.3.1  Equal  Noise  Power  with  Even  Spectrum 

An  important  case  is  that  for  which  the  noise  inputs  are  of 
equal  power  and  their  spectra  are  even  about  the  center  frequency  of  the 
band.  Under  these  assumptions,  the  cross-quadrature  correlation 

ottir-iont  v*  ic  70m  D6] 


co 


efficient  r is  zeroL  J and  o|  = a|  = a2,  implying  that  R = R'  = 0, 


sc  that  the  term  X in  (2-21)  is  zero.  Defining  the  input  SNR  for 
channels  1 and  2 by 
S 2 


h? 


A H 


1 " 2a7 


and 


c 2 
h2 

n2  2az 


the  general  expression  of  (2-21)  is  reduced  to  the  following  form: 


For  y > 0, 


Pe(y)  = 


^7  exp 


r^THpJ  'h3  " (Y)  hl\ 


x l 


r=0  L /8y/a2 
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and  for  y < 0 


where 


h3  = 2(T+pT  ^h1+h2+2hlh2  cos  29 ) 
h4  = 2(T-pT  (hl+h2"2hlh2  cos  29 ) 


and  where  20  = e.,-0,,  is  the  difference  in  phase  between  the  narrowband 

! C. 

fiol 

signals  at  the  multiplier  input  [see  Figure  2-1  and  (2~3a)].  Muraka1  J 
presents  results  corresponding  to  a further  specialization  of  this  case 
to  0 = 0 and  h-j2  = h|. 

2. 1.3. 2 No  Input  Signals  and  Equal  Noise  Power 

If  the  correlator  inputs  are  assumed  to  be  only  the  noise  of 
equal  power,  that  is,  h?  = h^  = 0 and  0^=0^=  a2,  the  pdf  (2-21)  is 
further  reduced  to 


1J  expt_  5’4h'p)'  ]* y ’ 0 

y < 0 


(2-23b) 

(2-23c) 

(2-23d) 


(2-24) 
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This  result  is  Identical  to  that  of  LezinL  J,  who  began  analysis  with 
these  assumptions. 

2.1.4  THE  PROBABILITY  DENSITY  FUNCTION  FOR  THE  OUTPUT  OF  THE 
STANDARD  OR  SQUARE  LAW  ARRAY  DETECTOR 
In  the  ensuing  chapters  we  shall  have  use  for  the  density 
function  for  the  standard  or  square  law  array  scheme,  as  we  have 
defined  it  in  Section  1.2.2.  If  the  input  to  the  squaring  device  is 
z(t)  = s(t)  + n(t)  = Z(t)  cos  [wt-4>(t)],  then  the  filter  output  is  given 
by  y = Z2/ 2 . 

As  above,  we  may  write  Z2  = z2  + z2  , where  z and  z„  are 

c s c s 

independent  Gaussian  variates  with  identical  variance  o2  and  means 
Scos0  and  SsinQ  respectively  if  s ( t)  = S(t)  cos  [uit-e(t)]. 

It  is  well  known  that  the  sum  of  the  squares  of  n independent, 
unit-variance  Gaussian  variates  with  means  m.  is  a noncentnl  chi-square 
variate  with  n degrees  of  freedom  and  noncentrality  parameter 

+ m|  + ...  + m2  [23].  Therefore,  x = -%r  is  noncer.tral  chi-square 
with  two  degrees  of  freedom  and  parameter  S2/o2,  and  the  probability 
density  function  of  y is 


" ^2)  y > °*  (2-25) 


where  h2  is  the  SNR  at  the  input  to  the  square  law  device. 


2.1.5  COMPUTED  RESULTS 

Judging  the  symmetric  array,  even  spectrum  case  to  be  one  of 
most  general  utility,  we  shall  present  computations  of  this  probability 
density  function  (2-23),  as  well  as  the  related  density  for  the  square 
law  array  (2-25),  in  order  to  show  the  effects  of  the  parameters  p, 
hi  = hz  = h2,  and  0 = (e^-Gg)^;  a,  quite  obviously,  is  a scaling  para- 
meter. In  Appendix  B the  pdf  (2-23)  is  shown  to  be  equivalent  to 


p60)  * 


. - [*(l-p)hS]"  [*(Up)h?]" 

exp  [-(h3  * H4>]  lQ  lQ  m! 


n! 


exP  rz 


y > 0 


pxc  r 2y  1 r__dfy 

exp  [oTl^Tj  “n  |_a"(l-p^)J 


, y < o, 


(2-26) 


* , » ^ 

where  the  polynomials  we  have  defined,  G (x)  = l (m  "~K)  have  the 

m k=Q  n k. 

unusually  useful  computational  property  = G^  + G^”  . The  form  (2-26) 
results  from  application  of  the  characteristic  function  method  to  the 
special  case  we  are  now  considering. 

For  computation  of  (2-26),  we  have  chosen  a "nominal  case"  of 
specific  values  for  the  parameters:  p=0,  o=l,  h=l , 6=0.  In  Figures  2-2 

through  2-5,  (2-261  and  (2-25)  are  plotted  under  varying  parameters.  It 

2 

is  seen  in  Figure  2-2,  that  h has  the  powerful  effect  of  changing  the 
pointed  curve  of  the  no  signal  case  into  curves  which  begin  to  approach 
the  familiar  Guassian. 
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FIGURE  2-2  MULTIPLICATIVE  CONFIGURATION  PDF;  SNR  VARIED 


FILTER  OUTPUT 

PROBABILITY  DENSITY  FUNCTION 
Multiplicative  Configuration 
(Equation  (2-26)) 
p=0,  h2  =1,  o2=l 
e varied 


FIGURE  2-5  STANDARD  CONFIGURATION  PDF;  SNR  VARIED 
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APPLICATION  TO  THE  ARRAY  MODELS 


All  subsequent  analysis  shall  be  based  on  the  application  of 
the  filter  output  probability  density  functions  for  the  "main"  case  in 
which  the  array  and  the  array  noise  input  processes  are  symmetric,  and 
the  noise  spectrum  is  even  about  the  narrowband  center  frequency.  This 
application  shall  consist,  in  this  chapter,  in  showing  the  explicit 
probability  density  functions  (that  is,  containing  all  the  pertinent 
parameters)  for  either  type  of  array,  including  computed  results. 

In  order  to  be  able  to  show  the  effects  of  inter-element 
correlation  on  the  probability  density,  we  have  taken  pains  in  Appendix  A, 
Section  3,  to  obtain  closed  form  representations  for  the  variances  and 
the  correlation  coefficient  at  the  array  sums.  The  resulting  complexity 
in  expression  is  considerable,  but  we  shall  be  quick  to  illustrate  with 
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2.2.1  THE  EXPLICIT  PROBABILITY  DENSITY  FUNCTION  FOR  THE 
MULTIPLICATIVE  ARRAY  MODEL 

For  the  multiplicative  array  model,  from  (1-5)  it  is  evident 
that  we  have  for  the  array  sum  signal  parameters, 

S]  = S2  = MAD(e/M;  M),  - e2)  = e = iMpsina, 

where  A is  the  signal  input  amplitude,  the  total  number  of  array  ele- 
ments is  2M,  a is  the  relative  bearing,  D(x;n)  = f^Tnf  » ^ = 

Referring  now  to  the  two  models  for  inter-element  correlation 
in  Appendix  A,  we  have  for  the  multiplicative  array, 


(2-27) 
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o2  = Ms2  + 2(M-l)as2 
a 

p " H + 2(M-1)a 


ADJACENT  CORRELATION 
a = correlation  coefficient 


,2  _ M2H2D2(0/M;  M) 
n M + 2{M-lja 

o2  = Ms2  + 2M(M-l)s2^— ?— 

M 4.,  -5  2 EXPONENTIALLY  DECREASING 

e He  (1-e  ) 

P bZ+2 (M-l ) ( b+e~  b-l7  CORRELATION 

Correlation  between 

^2  _ b2'-iH2D2(9/M;  M)  elements  m and  n « e”b|m-n| 

b2+2(M-l)(b+e"b-l) 

. A 

where  s2  = noise  variance  at  hydrophones 
H2  = SNR  at  hydrophones. 

Tn  order  to  keep  the  notation  under  control,  let  us  define  for  either 
model  of  correlation  the  two  parameters 
B = a2/Ms2 
C = PB. 


(2-28) 


(2-29) 


(2-30) 


Thus  we  may  conveniently  resort  to  the  white  noise  case  by  setting 
B = 1 and  C = 0.  With  these  definitions,  we  may  substitute  into  (2-23) 
to  obtain  the  explicit  probability  density  function  for  the  multiplica- 
tive array  filter  output: 
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Pg(y)  s P6(y;  a,  3.  M,  s,  H,  a or  b) 


1 _J  -2y  2MH2D2  6-Ccos20  , MH2D2sin26 
BM?  .Ms2<B+C)  ' 8 B2-C2  B 


'MsHDcosQ  B-Cy  T /4HDcos6v 

\ zJj  B/  r\  S(B+CJ 


)Lr/W^|ic\  y>0 


(2-31) 


Jjexp  — ^ 


2y 2MH2D2  B-Ccos29  + MH2D2cos29  B-C 


BMs2  [ms2(B-C)  B B2-C2 


v/MsHDsine  B+c\rt  /4HDsin9/=A  /-MH2D2cos20  B-c\ 

k\^TT) !^Ah  *9* y °* 


(2-32) 

To  illustrate  the  use  of  this  expression,  we  shall  give  four 
examples.  To  obtain  the  white  noise  case,  we  use  B = 1 and  C = 0,  with 
the  result 


p6(y)  = P6(y;  3.  m,  s,  h) 

= exp  {-2y/Ms2  - 2MH2D2  + W2D2sin20) 
Ms^ 


oo 

X | ^^te^Jir|W^^Lr(.MH2D2sfn2e)>  y > 0 (2-33) 


= -V  exp(2y/Ms2  - 2MH2D2  + MH2D2cos29) 


x ^ J Ir^4HDS3'n9>/C^Lr(-MH2D2cos29) , v < 0.  (2-34) 
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For  a second  example,  let  us  consider  a 10-element  multipli- 


cative array  with  adjacent  correlation.  Then  M = 5,  and  for  y > 0, 
P6(y;  a,  3,  5,  s,  H,  a) 

- 1 ovn  -2y  10H2D2  l£L6a-acos29 

5s4(l+1.6a)  5s  u+2.6a)  1+1 .6a  l+3.2a+l .56a^ 


5H2D2sin6  1+2 


1+1. 6a 


oo 

+2. 6a  V/2.5sHDcos8  l+.6a  \r 
l + .6a  ZA  y—  1+1. 6a/ 


w T /4HDcos6*^\,  /-5H2D2sin20  l+2.6a\ 
1r\F(l+l76a)/Lr\  1+1. 6a  l+.6a/ 


For  our  third  example,  let  us  find  the  density  for  an  array 
whose  noise  environment  has  been  determined  to  be  exponentially  decreas- 
ing, in  the  manner  in  which  we  have  defined  this  concept,  and  for  which 
experiments  yield  the  information  b = 1.  Thus  B = .74  and  C = .15. 

We  have  for  this  density,  if  y > 0, 


p6(y;  o,  3>  s,  H,  b-1 ) 


2MHCDC  .26+.74M-.15cos26 


Ms2(.26+.74M)  Ms2( .41+.74M)  . 26+.74M  .045+.34M+.55M2 


pp  p CO 

MjrD'  si n^O  .41+.74M  V/MsHDcos0  .11+.74MY 
.26+.74M  .11  + .74M  Zl  .26+. 741-1/ 

r=0  2/y 


T [4HDcos6/y  ].  /-MH 
r l sT74T  + . 74MJI  Lr  \ . 


-MH2D2sin29  .41+.74M\ 
.26+.74M  .11+.74M/’ 
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As  a fourth  example,  let  us  find  what  is  called  the  "boresight 


case,  where  a = 0.  Then  e = 0 and  D = 1 , so  that 


P6(y;  o,  0,  M,  s,  H,  a or  b) 

■ exp[jif<y/MsW)j  I^Tt)  y 5 0 (2‘37> 


BHs‘ 


exp 


2y MH_/?  Hi 

- ' B+CT  " B ) 


Ms  (B-C) 


, y < 0. 


(2-38) 


2.2.2  THE  EXPLICIT  PROBABILITY  DENSITY  FUNCTION  FOR  THE  STANDARD 
ARRAY  MODEL 

For  the  standard  array  model,  from  (1-4)  we  have  for  the  array 
sum  amplitude,  S = NAD(o'/N;  N),  where  we  have  used  e'  = iNBsina,N  = total 
number  of  hydrophones.  As  in  (2-28,  2-29,  2-30)  we  may  write  for  the 
array  sum, 

a2  = Ns2  + 2 (N-l )as2 

? o ? ADJACENT  CORRELATION  (2-39) 

.2  _ N H D (8  7N;  N) 
n N + 2 (N-l )a 


a2  = Ns2  + 2N(N-l)s2(--e?-  — ) 

h2  = b2NH2D2(9  '/N;  N) 
b2+2(N-l)(b+e’b-l ) 

2 2 

and  we  may  also  define  B'  = a /Ns  . 


EXPONENTIALLY  DECREASING 

(2-40) 

CORRELATION 


(2-41) 
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The  white  noise  case  corresponds  to  B'  =1. 

With  these  definitions,  we  substitute  in  (2-25)  to  rind  the 
explicit  probability  density  function  for  the  standard  array  filter 
output: 


psq(y)  = Psq(y;  O.  3,  N,  s;  H,  a or  b) 


(2-42) 


2.2.3  COMPUTED  RESULTS 

Results  for  the  array  probability  densities  (2-31,  2-3?)  and 
(2-42)  are  presented  in  Figures  2-6  through  2-14. 

For  the  multiplicative  array  model,  the  nominal  case  of  10 
elements  (M  = 5),  white  noise  (a  = 0),  unit  noise  variance  (s2  = 1), 
zero  bearing  (a  = 0),  half-wavelength  array  spacing  (B  = u),  and  zero 
dB  input  SNR  (H2  =1)  was  chosen.  In  Figure  2-G,  H2  is  varied,  yield- 
ing curves  very  similar  to  those  of  Figure  2-2,  although  we  are  now 
dealing  with  an  order  of  magnitude  higher  values  of  the  filter  output, 
since  each  array  sum  is  contributed  to  by  five  array  elements,  and  the 
normalization  we  have  chosen  is  at  the  elements. 

The  number  of  array  elements  (2M)  is  varied  in  Figure  2-7, 
showing  especially  well  the  effect  of  this  parameter  on  the  mean  and 
variance.  It  is  evident  that  the  number  of  elements  does  not  signifi- 
cantly change  the  shape  of  the  distribution,  as  does  the  SNR.  However, 
for  lower  input  SNR,  as  in  Figure  2-8,  a more  dramatic  effect  can  be 
observed. 
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FIGURE  2-7  MULTIPLICATIVE  ARRAY  PDF;  NUMBER  OF  ELEMENTS  VARIED 


FIGURE  2-9  MULTIPLICATIVE  ARRAY  PDF; 
ADJACENT-ELEMENT  CORRELATION  VARIED 
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FIGURE  2-14  STANDARD  ARRAY  PDF;  BEARING  VARIED 


The  effect  of  adjacent  inter-element  correlation  is  displayed 
in  Figure  2-9.  It  is  seen  that  correlation  among  the  array  elements, 
and  therefore  between  the  array  sums,  increases  the  variance  of  the 
filter  output  and  decreases  the  value  of  y at  which  the  peak  of  the 
distribution  occurs. 

Figure  2-10  demonstrates  a phenomenon  of  great  interest  to 
us  in  Chapter  4.  The  bearing  is  seen  to  act  as  a location  parameter 
for  the  distribution  of  the  filter  output.  The  manner  in  which  the 
bearing  slightly  affects  the  SNR,  through  the  directivity  function,  can 
also  be  observed  in  this  figure.  In  Figure  2-11  it  is  shown  that  for 
low  SNR,  bearing  ceases  to  be  a location  parameter  in  the  usual  sens0, 
since  the  peak  of  the  distribution  remains  at  y = 0. 

In  Figures  2-12  through  2-14,  the  distribution  of  the  filter 
output  of  the  square  law  or  standard  array  model  is  seen  to  behave  with 
parameter  variation  in  much  the  same  way  as  that  of  the  multiplicative. 
Here  we  have  used  H2  = 1,  s2  = 1,  6 = ",  a = o,  a = 0,  and  N = 10  as 
the  nominal  case.  In  Figure  2-10,  the  SNR  is  varied,  while  in  Figure  2-13 
and  2-14  :he  array  inter-element  correlation  and  the  bearing  are  varied, 
respecti\ely.  It  is  very  noticeable  that,  although  we  have  used  the  same 
number  of  elements  as  in  the  multiplicative  case,  the  magnitude  of  the 
filter  output  values  is  much  greater.  In  fact,  there  is  a factor  of 
four  involved,  as  can  be  understood  by  comparing  (1-4)  and  (1-5)  for 
N = 2M. 
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2.3  MEAN,  VARIANCE  AND  SNR  AT  THE  FILTER  OUTPUT 

2.3.1  MULTIPLICATIVE  ARRAY 

For  the  narrowband  multiplicative  configuration,  we  have  from 
Section  2.1 , 

y = (s1  + n^)($2  + n2)  - (terms  of  frequency  2u) 
'^lS1(n2ccos6l+n2ssin9l)  T S1S2cos(6)-e2)  + S2(nlccos02+n)ssine2) 

+ "lcn2c  + nlsn2s]  f2-43) 


so  that  E(y)  = po^  + ^^200520  = a^P  + h,h2cos2e)  (2-44) 

and  Var(y)  = jcy^O+p^-r2)  + ^c^S^fpo^e+rs^e)  + ifa^+oijs^) 


_ 1 _2_2n 
- SU-ji^Li 


- r 


+ ^:r*-j  r«2 1 pcus^6+rsi  nZG ; + n^  + 


..  . 

n2J  U-4b) 


where  9 = j (e^^)  and  we  have  used  (1-7)  and  the  well  known  moment 
properties  of  zero-mean  Gaussian  random  variables  ([16],  Section  7.3.2): 

E(uvwx)  = E(uv)E(wx)  + E(uw)E(vx)  + E(ux)E(vw)  (2-46) 

E(uvw)  = 0 (2-47) 

with  (u,  v,  w,  x)  being  any  combination  of  (n^c,  n^  , n2C,  n2s). 
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For  our  main  case  of  the  symmetric  array,  even  noise  spectrum, 
we  have  0l  = o2  = 0.  S1  = S^  = oh  *2,  and  r = 0.  Therefore,  for  the  main 


case, 

E (y ) = o2(p  + h2cos  2e) 

Var(y)  = la^[l  + p2  + 2h2(l  + p cos  2e)]. 

Further,  using  the  notation  of  (2-29),  2-30,  2-31),  we  obtain  expli- 
citly 

E(y)  = Ms2[C  + MH2D2 ( 0/M ; M)  cos  2e] 

Var(v)  = 5 M2s4[B2  + C2  + MH2D2 (0/M) (2B  + C cos  2e)]. 

We  may  use  these  expressions  to  show  the  SNR  at  the  filter 

r*i  fii 

output,  using  the  conventions  of  Middleton1  , Section  5.3: 


SNR  = 
0 


Var(y) 


2M2H4D4cos229 

B2  + C2  + MH2D2(2B  + Ccos20) 


2.3.2  SQUARE-LAW  ARRAY 

For  the  narrowband  standard  or  square  law  array  configuration, 
the  filter  output  is  a noncentral  chi-square  random  variable,  as  noted 
in  Section  2.1.4.  Since  the  mean  and  variance  of  x2(x;  2,  a2)  are  given 
([23])  by  E(x)  = 2 + a2,  Var(y)  = 4 + 4a  , we  have 

E(y)  = a2(l  + h2)  (2-53) 

Var(y)  = °4(1  + 2h2).  (2-54) 

Using  the  notation  of  (2-39,  2-40,  2-41),  we  write  explicitly, 


E(y)  = N s 2 [ B ' + NH2D2 (0 ' / N ; N)] 

Var(y)  = B'N2s4[B'  + 2NH2D2(o'/N;  N)]. 


(2-55) 

(2-56) 
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The  filter  output  SNR  then  becomes 


SNRo  . ff(y)Y;r$  H-oir 


n2h4d4 


B*  (B*  + 2NH2D2 


2.3.3 


COMPARISONS 


(2-57) 


It  is  interesting  to  compare  the  moments  and  SNRs  which  we 
have  just  shown.  We  shall  do  so  for  the  white  noise  case  in  which 
also  6=0  {boresight  case).  Denoting  means  by  E and  variances  by  V, 
we  have 


E$q  = Ns 2 ( 1 + NH2) 
Vsq  * N2;4(1  + 2NH2) 


n¥ 


1 + 2NH 


E = M2s2H2 
m 

V = iM2s4(l  + 2MH2) 
m 


SNRm 

m 


2M2H4 
1 + 2MH2 


(2-58) 


The  comparisons  we  wish  to  show  are  tabulated  as  Figure  2-15,  in  which 
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standard  and  multiplicative  arrays  having  the  same  total  number  of 
elements,  that  is,  N = 2M;  (b)  the  two  arrays  having  roughly  the  same 
beamwidth,  represented  by  N = 4M.  It  is  seen,  for  example,  that  the 
multiplicative  array  SNR  is  less  than  or  equal  to  that  for  the  standard 
array  in  all  these  cases. 


★ ★★★★★ 

Additional  references.  The  method  we  have  used  to  obtain  the  filter 
output  probability  density  functions  is  by  no  means  the  only  one. 
Perhaps  the  most  thorough  method  is  that  which  takes  into  account  the 
filter  transfer  function  of  a realizable  filter,  and  treats  the  filter 
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input  as  a time  series.  In  the  literature,  this  approach  is  often 
attributed  to  Kac  and  Siegert^1^,  and  usually  is  found  applied  to  the 
square  law  detector  [41-43].  An  exception  is  the  work  of  Lampard^^, 
who  studies  a multiplier/filter  configuration.  The  curves  in  MarcuirJ-^ 

f4 1 1 

and  in  Emerson1-  J are  especially  comparable  to  those  we  have  generated. 

Others  not  mentioned  in  the  text  who  have  used  the  same  method 
as  we  to  analyze  the  multiplier/filter,  though  less  generally,  are  listed 
for  reference  [44-46].  A thorough  treatment  of  the  output  SNR  of  such 

[471 

systems  is  that  of  Green1-  . In  these  works,  very  similar  probability 
density  curves  were  obtained  from  sometimes  very  dissimilar  expressions. 

For  those  who  wish  to  pursue  the  mathematical  statistics 
aspect  of  this  problem,  we  provide  a listing  [48-52]  of  what  we  con- 
sider profitable  reading  in  this  area. 
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CHAPTER  THREE 
DETECTION  PERFORMANCE 

Having  now  the  probability  density  functions  for  the  array 

processor  filter  outputs,  we  are  in  a position  to  evaluate  what 
f5l 

Tucker1-  J described  as  "...  the  difficult  question  of  whether  the 
better  signal /noise  performance  [of  one  detector  over  another]  is  really 
indicative  of  a higher  probability  of  detection."  That  is,  we  are  able 
to  calculate  directly  the  probabilities  of  detection  for  the  standard 
and  multiplicative  array  models,  as  functions  of  the  SNR  and  of  corres- 
ponding probabilities  of  false  alarm. 

3.1  DETECTION  CRITERIA 

The  above  assertion  relies  on  various  definitions  and  assump- 
tions. Our  analysis  assumes  the  usual  Neyman-Pearson  arrangement^26^, 
shown  in  Figure  3-1,  in  which  the  value  of  the  filter  output  is  subjected 
to  a "test."  If  the  filter  output  (y)  is  greater  than  a threshold  (T), 
then  the  decision,  "signal  present,"  is  made.  If  y < x,  the  decision 
"no  signal"  is  made. 

Since  due  to  the  noise,  y is  a random  variable,  this  decision 
strategy  is  susceptible  to  two  kinds  of  error.  One,  called  "false 
alarm,"  occurs  when  a threshold-crossing  due  to  noise  only  causes  a 
spurious  "signal  present"  decision.  The  other  error  lies  in  missing  a 
weak  signal  because  a threshold-crossing  did  not  take  place. 

The  problem  in  this  type  of  strategy  is  to  choose  the  threshold 


such  that  an  acceptable  probability  of  error  is  maintained.  Often  there 
is  a high  cost  attached  to  a false  alarm  (for  example,  an  expensive 
weapon  fired),  so  the  usual  practice  is  to  consider  the  threshold  as 


Various  types  of  detection  systems  then  are  evaluated  in  terms 


! 


of  the  probability  of  detection  (P^)  achieved  for  a given  threshold, 
as  a function  of  the  SNR  when  a signal  is  present.  Therefore  the  trade- 
off between  Pp^  and  P^  is  of  great  interest,  and  it  is  standard  practice 
to  exhibit  this  tradeoff  in  the  form  of  "receiver  operating  character- 
istics," in  which  P^  is  plotted  as  a function  of  Pp^  and  the  input  SNR. 
3.2  CALCULATION  OF  FALSE  ALARM  PROBABILITY 


Since  a false  alarm  occurs  when  the  filter  output  y exceeds 
the  detection  threshold  t when  there  is  no  signal  present,  the  probability 
of  false  alarm  is  given  by 
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3.2.1 


MULTIPLICATIVE  MODEL 


Usinq  (2-26)  for  p,(vlH=0),  we  have  for  the  false  alarm 
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probability  for  the  multiplicative  array  model. 


00 


f V exP 

--?y- 

0^ 

o(i+p) 

(t>0) 


= i(l+p)  /dxe 


-x 


u2(l+p) 


PpA  = -Hl+p)  exp 


■2t 


K(1+p) 

or,  to  use  the  notation  of  Section  2.2.1, 
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3.2.2  SQUARE  LAW  MODEL 

Using  (2-26)  for  pgq(y),  we  obtain  for  the  false  alarm 
probability  for  the  square  law  or  standard  array  model, 

Pfa  = f <iy  ^2  exp(-y/°2)  fr*0) 

= exp(-T/cr2)  (3-4) 

or,  in  terms  of  Section  2.2.2, 

PFA  = exp(-T/B'Ns2).  (3-5) 

3.2.3  COMPUTED  RESULTS 

Figure  3-2  illustrates  the  simple  manner  in  which  the  false 
alarm  probability  varies  with  the  threshold.  It  is  clear  that  the 
standard  array  model  requires  a higher  detection  threshold  than  the 
multiplicative  with  the  same  number  of  elements  to  achieve  a similar 
false  alarm  rate.  However,  this  statement  has  to  be  qualified  further 
for  general  application,  since  the  models  do  not  take  into  account  the 
attenuations  and  gains  of  actual  systems,  and  therefore  the  magnitudes 
of  the  filter  outputs  may  not  correspond  in  the  same  way  as  they  do 
here. 

Nevertheless,  the  mean  and  variance  of  the  filter  output  due 
to  noise  alone  tends  to  be  smaller  in  the  multiplicative  configuration, 
allowing  a relatively  lower  threshold  (compare  equations  (2-51)  and 
(2-56),  (2-52)  and  (2-57)  for  N = 2M). 


3.3  CALCULATION  OF  DETECTION  PROBABILITY 

Detection  of  a sigral  which  is  present  occurs  when  the 
filter  output  exceeds  the  detection  threshold.  Therefore,  the  proba- 
bility of  detection  is  given  by 


PD  =/dy  I H>°)  • 


(3-6) 


3.3.1  MULTIPLICATIVE  MODEL 

Using  (2-26)  for  Pg(y),  for  the  multiplicative  array  model 

we  have  for  the  probability  of  detection, 
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To  obtain  the  result,  we  take  the  integral  in  (3-7)  separately: 
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where  r(n,x)  is  the  incomplete  Gamma  function.  Using  [13],  formula 
8.352.2,  we  continue: 
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where  we  have  taken  b = 2/o2(1  + p)  and  c = 4/o2(l-p?). 

Substituting  (3-9)  in  (3-7),  we  find  for  h^h^^h, 
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or,  in  the  notation  of  Section  2.2.1, 
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3.3.2  SQUARE  LAW  MODEL 


Using  (2-25)  for  PSq(y)»  the  probability  of  detection  for 
the  square  law  or  standard  array  model  becomes 

03 

PD  = /dy^  exp(-h2-y/a2)  IQ(2h^/o) 

T 

00 

= J'x  exp(-h2-£x2)  Iq ( v^hx)  dx  (3-12) 

*^r/a 

or  PD  = Q(v^h,  ^f/o),  (3"13) 

where  Q is  Marcum's  Q-function.  In  terms  of  Section  2.2.2,  the  proba- 
bility of  detection  is 

PD  = Q(HD*WF  , J 2x/B  Ns2) . (3-14) 

3.3.3  COMPUTED  RESULTS 

The  behavior  of  the  probability  of  detection  at  the  filter 
outputs  as  a function  of  threshold  and  array  input  SNR  is  shown  in 
Figure  3-3  for  N = 2M  = 4 and  in  Figure  3-4  for  N = 2M  = 10. 

It  is  seen  that  in  general,  the  standard  array  model's 
detection  probability  is  higher  than  that  of  the  multiplicative  for 
the  same  threshold,  subject  to  the  qualifications  cited  in  Section  3.2.3. 
This  is  due  to  the  fact  that,  for  the  same  number  of  array  elements,  the 
bulk  of  the  standard  probability  density  is  concentrated  about  a higher 
mean  value  than  for  the  multiplicative  case.  Another  factor  is  that  the 
detection  probability  for  the  square-law  or  standard  processor  is 
always  unity  for  zero  threshold. 
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SYSTEM  COMPARISON:  RECEIVER  OPERATING  CHARACTERISTICS 


In  the  previous  sections  we  have  exan  Ined  separately  the 
probabilities  of  false  alarm  and  of  detection  for  the  two  array  models 
we  are  considering.  It  was  seen  that  the  multiplicative  model  has  a 
more  desirable  false  alarm  rate  for  the  s tme  value  of  detection  threshold, 
while  the  square-law  or  standard  model  yields  a higher  probability  of 
detection.  However,  this  apparent  tradeoff  is  blurred  by  the  unsuit- 
ability of  the  detection  threshold  as  a basis  for  comparison. 

A more  informative  method  of  comparing  such  systems  is  the 
construction  of  "receiver  operating  characteristics,"  which  eliminates 
the  threshold  altogether  by  converting  the  parametric  relations 


PFA  ~ PFA^ 

PD  = Pd(t,SNR) 


(3-15) 


into  the  relationship  PQ  = ^(Ppp ,SNR) . (3-16) 

Because  with  the  receiver  operating  characteristics  (ROC) 
representation  two  different  threshold  detectors  are  more  directly 
comparable,  a performance  index  such  as  "minimum  detectable  signal"  can 
be  used  to  pronounce  one  the  better  detector  with  some  confidence  in 
the  generality  of  this  kind  of  statement.  (The  minimum  detectable  sig- 
nal is  defineo  as  the  SNR  required  to  insure  a given  probability  of 
detection  for  a certain  false  alarm  rate,  also  given.)  Helstrom^^ 
cites  additional  uses  of  the  ROC;  for  example,  determination  of  para- 
meters for  Bayes  and  minimax  criteria  detection  scheme^. 

The  ROC  curves  for  our  two  models,  for  the  white  noise,  bore- 
sight  case,  are  given  in  Figure  3-5  for  N = ?M  = 4 and  in  Figure  3-6 
for  N = 2M  - 10.  The  fact  that  the  standard  array  processor  ROC  curves 
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are  consistently  above  those  of  the  multiplicative  indicates  that  the 
standard  array  model  is  a better  detector  in  this  case.  They  are 
quite  close  together,  however,  and  the  minimum  signals  for  the  two 
configurations  would  appear  not  to  differ  by  more  than  1 dB. 

The  more  usual  presentation  of  detection  probability  as  a 
function  of  SNR  is  given  in  Figure  3-7  for  four-element  arrays  and  in 
Figure  3-8  for  ten-element  arrays,  for  various  values  of  false  alarm 
probability.  These  curves,  of  course,  contain  the  same  information 
as  Figures  3-5  and  3-6,  and  are  given  here  for  the  convenience  of 
those  who  are  more  familiar  with  this  presentation. 

Arndt^^  compared  SNR  and  angular  resolution  (beamwidth) 
properties  of  various  array  configuration,  including  linear  arrays  of 
the  type  we  are  studying.  Assuming  that  probability  of  detection  is 
proportional  to  output  SNR,  he  is  willing  to  declare  that  additive 
(conventional  or  "standard")  processing  is  the  better  detection  strat- 
egy, although  for  the  price  of  3 dB  in  SNR,  split-beam  multiplicative 
processing  obtains  superior  bearing  "estimation"  (tracking)  capability. 
Arndt  concludes  that  this  "price"  is  too  high  where  detection  is  para- 
mount, although  he  does  not  consider  the  false  alarm  rates  of  the  two 
configurations. 

Our  results,  while  not  directly  comparable,  suggest  that 
more  thorough  examination  may  reveal  that  the  cost  in  detection  proba- 
bility (vs.  false  alarm  rate)  is  less  than  indicated  by  the  difference 
in  SNR. 
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CHAPTER  FOUR 


MAXIMUM  LIKELIHOOD  BEARING  ESTIMATION 

In  Chapter  3 we  have  generated  characterizations  of  signal 
detection  performance  based  on  the  probability  density  functions  of 
Chapter  2.  At  this  stage,  we  have  provided  the  analytical  equipment 
to  answer  a good  many  of  the  questions  raised  in  the  literature  regard- 
ing the  relative  merits  of  the  two  array  configurations  we  are  consi- 
dering. However,  before  making  a general  comparison,  we  should  like 
to  include  another  important  consideration. 

In  addition  to  detecting  the  presence  or  absence  of  a signal, 
we  desire  to  know  the  bearing  of  the  signal  source  precisely.  The 
purpose  both  of  steering  the  array  and  of  seeking  to  improve  its  direc- 
tivity by  summing  the  outputs  of  several  hydrophones,  is  to  obtain  a 
filter  output  whose  value  can  be  depended  upon  to  show  a peak  when  the 
array  is  steered  in  the  direction  of  cue  signal  source.  The  precision 
of  this  technique  is,  of  course,  related  to  the  tracking  (noise) 
conditions  and  also  to  such  factors  as  the  accuracy  of  the  steering 
and  peak-detecting  operations.  (Often  the  "peak-detector1  is  a sonar 
operator, ) 

The  precision  with  which  bearing  is  determined  can  be  calcu- 
lated in  some  fashion  for  a given  realization  of  an  array  system. 
However,  what  we  are  interested  in  discovering  here  is,  given  the  array 
processor  configuration,  what  is  its  ultimate  capability  in  establishing 
bearing?  An  answer  to  this  question  could,  for  example,  motivate  a 
designer's  decision  to  improve  the  performance  of  an  array  subsystem 


or  component,  if  the  configuration's  ultimate  capabilities  were  not 
nearly  being  realized. 

In  Chapter  5 we  show  calculations,  based  on  statistical  con- 
cepts, which  bound  the  minimum  bearing  error  attainable.  In  this 
chapter  we  derive  forms  for  bearing  estimators,  whose  performance  we 
shall  later  compare  with  the  theoretical  bounds. 

4.1  MAXIMUM  LIKELIHOOD  ESTIMATION 

From  (1-4)  and  (1-5),  we  know  that  the  array  filter  output, 

if  there  is  no  noise,  behaves  in  a known  way  as  a function  of  the 

bearing,  so  that  we  can  write  y(t)  = y ( t ; a).  Assuming  we  have  only 

this  filter  output  to  work  with,  we  could  determine  the  bearing  by 

inverting  this  function:  a = a(y).  However,  since  in  general  there  is 

noise  present,  the  result  of  this  operation  would  have  to  be  considered 

* 

an  estimate  of  the  bearing,  written  a , and  illustrated  in  Figure  4-1. 

And  there  would  be  error  associated  with  the  estimate.  Much  attention 
is  given  in  statistics  to  the  problem  of  finding  estimators  which  yield 
the  most  accurate  (least  mean  square  error  or  minimum  variance)  value 
of  the  parameter  being  measured.  Since  we  have  the  probability  density 
function  at  the  filter  output,  we  may  make  use  of  the  results  of  the 
statistical  theory  of  estimation. 

According  to  Cramer1  , the  most  important  general  method  of 
finding  estimates  is  the  maximum  likelihood  method.  Briefly  stated, 
this  method,  as  applied  to  our  problem,  consists  of  considering  the 
maximum  likelihood  bearing  estimate  (a)  to  be  that  bearing  for  which  the 
known  value  of  the  filter  output  (y)  is  the  most  likely  value.  That  is, 

a = | a:  p(y;  ah  p(y;  a)  L (4-1) 
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Therefore,  we  may  determine  the  maximum  likelihood  estimate  by  finding 
the  bearing  which  maximizes  p(y:  a): 


a - {a(y):  3p^;  a?  ° °=  0,  ^ p0)-  (4‘2) 


The  importance  of  this  method  lies  in  the  property  that  [29] 
under  fairly  general  assumptions,  the  maximum  likelihood  method  will 
produce  asymptotically  efficient  estimates,  that  is  estimates  whose 
variance  approaches  the  theoretical  minimum  for  certain  conditions. 

We  shall  not  be  computing  the  secor.j  derivative  indicated  in 
(4-2)  since  we  have  seen  in  Chapter  2 that  the  bearing  (suitably 
restricted  in  value)  acts  somewhat  like  a location  parameter,  and  thus 
we  expect  that  the  solution  to  the  maximum  likelihood  equation  3p/9ct  = o 
to  also  be  the  value  of  u which  causes  the  given  value  of  y to  corres- 
pond to  the  peak  of  the  distribution.  Although  for  small  SNR  this 
location  influence  of  a changes  in  that  the  peak  value  remains  at  y = 0, 
we  shall  assume  that  the  solution  of  the  maximum  likelihood  equation 
continues  to  be  a maximum. 


4.2  THE  MAXIMUM  LIKELIHOOD  BEARING  ESTIMATE  FOR  THE 

STANDARD  ARRAY 

Using  (2-25),  the  probability  density  function  of  the  filter 
output  of  the  standard  or  square  law  array,  we  have  the  following 
equation  for  the  maximum  likelihood  bearing  estimate: 
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where,  from  Section  2.2.2,  we  have  the  notation 


h ^ h(a)  HD(i-gsina;  N). 


Equation  (4-5)  corresponds  to  a special  case  of  the  maximum  likelihood 
equation  for  the  estimates  of  the  noncentrality  parameter  of  a non- 

r °3i 

central  chi-square  variate,  as  given  by  MeyerLJ  J. 

4.2.1  THE  HIGH  SNR  CASE 

As  shown  in  Figure  4-2,  the  ratio  of  Bessel  functions  ori  the 

left  hand  side  of  (4-5)  approaches  unity  for  large  values  of  the  argu- 

2 

ment,  satisfied  by  large  SNR  (H  ).  That  is, 


1^ (x)/I0(x)  =1,  x > 3. 


Thus  for  H^D^y  > 2.25(B')^s^  or  h^y  > 2.25o^,  we  obtain 


or,  from  section  2.2.2, 


h = /y/o. 


y = s^N^H2D^(T0sina;  N) 


= s^sin^(iHiBsin&)/sin^(26sina)  , 


-H4- 


an  implicit  form  which  we  shall  find  useful  in  Section  4.4.  From 
Appendix  C we  have  the  approximate  inverse  or  explicit  forms 


# < h'Hs 


(4-10) 


= 0,  > NHs 


where  the  approximation  assumes  i(N-l)B  sin  a <<  AZ.  Since  the  beam- 
width  is  given  at  most  by  ||N8  sin  a | <tt , this  assumption  is  a valid  one. 

It  is  instructive  to  note  that  (4-9)  is  equivalent  to  taking 

y(£)  = [E(y)  - E(y|H  = 0)f=(\  (4-11) 


that  is,  the  maximum  likelihood  estimate  for  the  high  SNR  case  can  be 
interpreted  as  the  operation  which  assumes  that  the  filter  output  is  its 
mean  value,  less  the  1 ffects  of  noise. 

4.2.2  THE  LOW  SNR  CASE 

For  small  values  of  the  argument  in  (4-5),  the  ratio  of  the 
Bessel  functions  can  be  approximated  by 

I1 (x)/IQ(X)  = .584x  e"-26x,  x < 3.  (4-12) 

lhus  for  H2D2y  < 2 . 25 ( B 1 )2s2,  we  obtain 

A 

.584  jvy  exp  = c//y  (4-13) 
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h = -5-=  ln(^) 

,52/y  v o ' 


or,  from  section  2.2.2, 


D(^N3sina;  N)  = 


.52NH/y 


alj_  1n/ 

m\N+2  N-l  a 


As  for  (r-10),  we  use  (D-9)  to  write 


, , , 48(1 -D) 

a(y)  = * arccos  1 jo 

32(N2-1) 


(4-14) 


(4-15) 


i i,  48  r, 

= + arccos  1 — * — 5 1 - 

p (N-l ) 


fN-* 2 (N-l  )a]  s ^ 1 . 1 68y/ 
/52NH  y ln  N+2{S-l)a 

(4-16) 


THE  MAXIMUM  LIKELIHOOD  BEARING  ESTIMATE  FOR  THE 
MULTIPLICATIVE  ARRAY 

For  the  multiplicative  array  model,  the  maximum  likelihood 


equation  may  be  written 


_ / \ ia=c.  „ a=a 

3P6(y»  a,  3P6  0 3P6 

"is - 80  d^r  = *M3c0Sa~lG 


(4-17) 


Assuming  that  cos  a / 0,  we  have,  using  9 = 0(a), 


(4-18) 


4-8 


From  [13],  formulas  8.972.1  and  9.212.1,  we  observe  that 


e\J(-x)  = (rJl')lF1<r+k+,.|c+l;x>- 


(4-W) 


Substituting  this  expression  in  (2-23),  we  have  for  the  probability 
density  of  the  filter  output  for  h-j  = h2  = h 

QO 

p6(y;  a)  = g(6)  = a1exp[-a2h2(l-pcos26)]  a!J(hcos0)r 

r=0 

x Ir(a4hcos0)  ^(r+1,1;  a5h2sin20), 


(4-20) 


where  we  have  tried  to  simplify  the  notation  somewhat  by  writing 

aj  = a~2exp [ - 2y/a2 ( 1 + a) ] , a2  - 2/ ( 1 - p2 ) , a3  = (l-p)o/2/y,  a^  = 4/v/a(l+p), 

and  a5  - (1+p)/(1-p). 

Performing  the  required  differentiation  of  (4-20)  with  respect 
to  e,  we  obtain 

g’(9)  = a1exp[-a2h2(l-pcos20)] 

00 

x -2a2h[h ' (l-pcos20)+phsin20]  ^ a3(hcose)rlr-(  )^F^(  ) 

1 « r=0 

+ (h'cos0-hsin6)  ^ ra3<hcos0)r_1  I^(  ) ^F^(  ) 
r=0 

00 

+ -/a4(h'cos0-hsine)  ^ a3(hcos0)r[Ir+1  ( )+Ir.-,(  J^Fjf  ) 


2aghsin0(h'sin0+hcos6)  ^ (r+1  Ja^hcose)*  Ir(  ) 

x ^F^(r+2,2;agh2sin20)  (4-21 y 
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I 


In  (4-21)  we  have  used  h'  = dh/de  and  have  suppressed  the  arguments  of 
functions  whenever  they  are  the  same  as  in  (4-20).  By  combining  summa- 
tions and  noting  that  rlr(x)  * x[Ir_^(x)  - Ir+^(x)],  we  may  reduce  (4-21) 
to 

00 

g ' (0)  = a^exp[-a2h2(l-pcos20)]  ^ (a3hcos0)r 

r=0 

x -2a2[h‘(l-pcos20)  + phsin20]Ir(a/,hcos0)1F^ (r+1 ,1  ;a5h2sin20) 

2 2 

+ ai,(h'cos0-hsin0)lr_^(a4hcos0)  jF^r+l  ,l;aj-h  sin  0) 

(4-22) 

+ 2aj-(r+l )hsin0(h'sin0+hcos0)Ir(a^hcos0)^F1 (r+2,2;a5h2sin20)  . 

A slightly  more  compact  form  may  be  realized  by  using  (419)  to  return  to 
the  Laguerre  polynomials: 

00 

g'(0)  = a1exp[-a2h2(l-pcos20)+a5h2sin20]  ^ (a3hcos0)r 

r=0 

x -2a2 [h ; (l-pcos20)+phsin2O]Ir(a4hcos0)Lr(-a5h2sin20) 

+ a4(h'cos0-hsin0)lr_1  (a/.hcos0)Lr(-a5h2sin20) 

«■  2a5hsin0(hsin0+hcos0)Ir(a4hcos0)L|>(-a5h2sin20)  . (4-23) 

Our  task  is  to  solve  the  equation  g 1 ( 0 ) = 0,  for  e = e(y). 

However,  this  we  have  not  been  able  to  do  in  general.  Rather,  we  next 
show  a small  SNR  approximation  to  the  solution,  and  for  higher  SNR,  a 
different  approach. 
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4.3.1  SMALL  SNR  CASE 


For  h sufficiently  small,  we  may  find  an  approximate  solution 
by  taking  the  first  term  of  the  series.  Thus  g‘(e)  = 0 becomes 


0 * -2ag[h' (l-pcos20)+phsin20]Io(a4hcose) 

+ a^(h'cos0-hsin0)I^(a^hcos0) 

+ 2aghsin0(h'sin0+hcos0)I^(a^hcos0) 


14-2*} 


Dropping  the  hh'  terms  and  supplying  the  expressions  represented  by  a^ 
and  a^,  we  obtain 

, Ar:  I,[4hcos0v57a(l+p)] 

— '—x  [h'(l-pcos20)+phsin20]  = 'cosB-hsinO)— 

1-p2  p'  Io[4hcos0v/y/a(l+p)] 

(4-25) 


For  convenience  we  use  a different  approximation  to  the  ratio  of  the 
Bessel  functions  than  (4-12).  From  Figure  4-2  we  see  that  if  x is  quite 


small,  then 


Mx) 


yxi = 2 


* T , X < 1. 


(4-26) 


With  this  approximation,  (4-25)  can  be  manipulated  to  yield  the  implicit 


estimator  form, 


A A A 


o / 1+p \ h1 (l-pcos20)+phsin20  . 

y " a a 1 1 i /\  ^/\ 

hcosG  v ; h'cos0-hsin0 


(4-27) 


In  the  terms  of  Section  2.2.1, 
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y 


(4-28) 


2 3/2  ~ ^ a. 

= s f /B+C  \ Dl(9/M;M)(B-Cco$2e)+CD(6/M:M)sin2e 

HD(e/M;M)  'B'C'  D'(8/M;M)cose-D(e/M;M)sine 


where  D'(0/M;M) 


dD  = 1 

de  “ Msin(0/M) 


* sTn0  Ccos0-Dcos(0/n)].  (4-29) 

For  white  noise  (B  = 1 and  C = 0),  the  implicit  estimator  form  (4-28) 
becomes 

A O A A A 

..  _ s O' s cose-Dcos(e/M)  /4-30) 

HD(D'cose-Dsine)  HD  cos2e-Dcosecos(e/M) 

For  small  e (within  the  array  major  lobe),  we  may  approximate: 


O A0  A /NO  O 

sZ  1-0  /2-D(l-0  /2M  ; 

= 

A AA 

1-D-0V2 

(4-31) 

y ~ ^ ^9  ^9  9 . 

HD  1-20  -D(l -0V2)(  1-0  /2M  ) 

HD 

/\  aA 

1-D-30V2 

or 

, A /N 

y = sV3HD. 

(4-32) 

From  this  we  obtain  the  explicit  estimator  form  for  the  white  noise, 
small  SNR  case, 


a (y)  = iarccos 


48 


p(M-T) 


0-w)] 


4.3.2  A NONCENTRAL  CHI-SQUARE  APPROXIMATION 


(4-33) 


For  SNR  other  than  very  small,  rather  than  work  with  (4-23) 
we  shall  follow  the  procedure  of  approximating  the  density  function  of 
the  filter  output  y by  one  which  has  a closed  form,  and  ihen  deriving 
the  maximum  likelihood  bearing  estimate  for  the  approximate  distribution. 
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Over  a wide  range  of  values  for  SNR,  we  have  found  that  a good 
approximation  to  the  probability  density  function  (2-23)  for  the  multi- 


plicative array  model  is  the  density  function  of  an  equivalent  non- 
central chi-square  random  variable,  chosen  such  that  the  mean  and  variance 
of  the  approximate  distribution  are  the  same  as  the  original. 

Let  x = d + cy/o2  be  a noncentral  chi-square  variate  with  two 
degrees  of  freedom  and  noncentrality  parameter  A2.  Then  we  have 


Py(y)  = cpx(cy/a2+d)/o2  = exp[-lcy/a2+d+A2)/2] 


2a 


x IQ(Xv/cy/a2+d),  cy/a2+d  > 0. 


(4-34) 


From  Section  2.3.1, 


2 2 

E(y)  = a2(p+h2cos29)  = [E(x)-d]  = U2+XZ-d) 

c c 


14-35) 


4 4 4 2 

Varly)  - ^ [l  + p2+2h2(l+  cos28)]  = ^Var(x)  = • (4-36) 

c 

In  terms  of  c we  can  solve  (4-35)  and  (4-36)  to  yield 


A2  = ~ [l+pt+2ht(l+pcos2e)3  - 1 


.2,0,2 


d = 1+— g [l+p2+2h2(l+pcos20)3  - c(p+h2cos29) 


(4-37) 

(4-38) 


An  illustration  of  this  approximation  is  given  by  Figure  4-3. 
Again  assuming  cos  a t 0,  (4-18)  becomes 

A 

18=0 


98 


= 0. 


(4-39) 


For  Py(y)  as  given  in  (4-34),  the  maximum  likelihood  equation  then  is 


FIGURE  4-3 
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A A A A.  f A A. 

( d ' +2XX ' ) IQ ( Xv  cy/o  +d ) 

A I JJ  A A 7—  A A AA  / O a" 

= I-|(Xvcy/a  +d)  [2XVcy/o  +d  + Xd1/  vCy/a +d  ] 


or,  using  x = cy/o2+d. 


IT(X/t) 

I0(x*6T) 


A A A 

d1  + 2XX1 
2X'/x  + Xd'/^x 


Taking  the  ratio  of  the  Bessel  functions  to  be  unity,  we  obtain  the 
quadratic  equation, 

A A AA  AA 

2X'x  - (d'+2XX' )/x  + Xd'  = 0, 


with  the  solutions 

A A A O A A A A Art 

v$T  = d'/2X',X  or  cy/c ? + d = Cd'/2X ' X\ 


Putting  c = 4 into  (4-37)  and  (4-38),  the  first  solution  can 
be  expressed 


_ 

a2  ' 16 


h'[1-(l-p)cos26]  + h(l-p)sln2B 

A A ' A 

h'U‘tPCOs26)  - hpsinZB 


3 + 2p(2-p) 


+ 4fi  [l-(1-p)cos20]. 


Note  the  slight  resemblance  of  this  solution  to  (4-27).  The  second 
solution  has  the  much  simpler  form 

y2  = o n cos20  - o^(i-p) 

= Ms2[MH2D2cos26+C-B/2]. 


(4-40) 


(4-41) 


(4-42) 


(4-43) 


(4-44) 


(4-45) 


I 


1 


This  solution  also  has  the  form,  or  nearly  so,  of  E(y)  - E(y|H=0). 

For  this  reason,  and  for  its  simplicity,  we  shall  restrict  our  atten- 
tion to  it.  An  additional  motive  for  preferring  (4-45)  is  that,  as  the 
SNR  becomes  very  large  and  P6(y)  approaches  the  normal  distribution, 
(4-45)  corresponds  to  the  well  known  maximum  likelihood  estimate  for  that 
case  [57]. 


To  obtain  an  explicit  form  a(y)  for  the  estimate,  we  first  note 
that  (4-45)  can  be  written 

D2(0/M;  M)cos20  = = --rfJT  + I - c)  (4-46) 

M sin  (0/M)  MH  vMs*  L ' 


sin220 


4A 
Sin  0 


4M2sin2(8/M)  M2sin2(e/M) 


= D2(0/M;  2M)  - M2sin2(0/M)D4(o/M;  M) 


O 

r»t-  / . /ii 

* u (0/n;  cv\) . 


(4-47) 


Then,  using  (D-9)  and  recalling  that  0 = {M&  sin  a,  we  find 


a(y)  = iarccos 


48 

82(4M2-1) 


/l  /JL+  B 


4.4  DISTRIBUTIONS  OF  THE  ESTIMATES 


(4-48) 


The  estimate  a(y),  being  a transformation  of  a random  variable, 
is  itself  a random  variable,  whose  probability  density  function  is  given 
by 


A A A 

PA(a;  a)  = py[y(a);  a]  |dy/da|. 


(4-49) 


1 
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For  example,  for  the  multiplicative  array,  high  SNR  case,  the  estimate 
(4-45)  has  the  density  function 


p^(<x;  a)  = Pg[cr^h2cos20-cr2(i~p) ; a]2G2h|h'cos20-hsin20|0‘ , (4-50) 


where  we  continue  to  use  the  shorthand  notations  h = h(0),  0 = 0(a). 

Often  at  this  point  the  analyst  will  attempt  to  provide  an 
expression  for  p^(a;a)  in  the  form  p^(a;a)  = f(a-a).  However,  since  the 
expressions  for  the  density  functions  are  so  complex,  this  objective 
requires  much  approximation,  and  there  is  little  to  be  gained.  We  shall 
be  satisfied  here  to  note  that  these  densities  are  readily  computed  from 
(4-49).  Instead,  this  effort  will  be  put  into  showing  the  means  and 
variances  of  the  estimates  for  the  high  SNR  cases  so  that  their  asumptotic 
behavior  may  be  examined  in  relationship  to  the  theoretical  bounds  in  the 
next  chapter. 

4.4.1  MEAN  AND  MEAN  SQUARE,  THE  STANDARD  ARRAY  ESTIMATE 

The  high  SNR  bearing  estimate  for  the  standard  array,  given  by 
(4-10),  has  the  mean  value 

nW 

E(a)  = a"2  J dy  arcsin[K\/l-/y/NHs]exp(-h2-y/a2)I0(2hv^7o) 

0 

1 

= k J dx  arcsin[K’/l-v6<]exp(-h2-kx)lQ(2ht/H) 


= 2k  j dx  x arcsin[KvT^exp(-h2-kx2)lQ(2hxA) , (4-51) 

G 
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r 


U 


il 


ri 

: 1 


= 1 EL  . k2/n2. 


where  we  employ  K = and  k = = h /D  (ipslnnj  N).  Approximately 


1 


A 2 2 

we  have  E(a)  = 2kKe"h  j dx  <Tkx  x/MT  IQ(2hx»^) 

0 


-h2  V V (kh2)n  (-k  f h,  „2n+2m+l 


* 2kKe~"  £ £ 
n=0  m=0 


) f 

w -it  Jdx  x 


/TiT 


(4-52) 


n=0  m=0 


where  B is  the  beta  function,  B(u,v)  = r(u)  r(v)/r(u+v).  By  a simple 
change  of  summation  index  we  may  write 


2 JH  ,,.L2xn  , ,.\m-n 


EW  - ^ B(3/2,  2«H-2) 


m=0  n=0 


* 2kKe‘h‘  £ 3(3/2.  2m+2)  Lm(h2} 


(4-53) 


m=0 


By  a similar  process  we  find  that 


E[(i)2]  = 2kK2e'ht  ^ B(2‘  2(11+2  > Lm(h2)- 

m=0 

The  beta  functions  B(3/2,  2n+2m+2),  as  shown  in  Figure  4-4, 


(4-54) 


can  be  approximated  by  the  geometric  progression 


B(3/2,  2n+2rr, - 2)  = B(3/2,2)x( .415) 


n+m 


(4-55) 
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Substituting  this  approximation  in  (4-52),  we  obtain 


EM  - j J L4|»££t^a£B(V2.  2) 

n-C  m=0 

* 2kKB(3/2,  2)exp(-h2-.415k)  In(2h/3l5k) 


or  E(o)  * kKB(3/2,  2)exp[-(hV375k) £]/v  ir/74T5kh, 


(4-56) 

(4-57) 


where  we  have  used  the  large  argument  relation  for  the  Bessel  function, 
Iq(x)  = ex//2ir x . Substituting  for  K,  k,  and  h,  we  find 

EM  - -r#=  (4  )V2expr-NH2(D-.644)2/B-3.  (4-58) 

BvWM)  V B ; 

A 2 

In  a similar  manner  we  approximate  E[(ct)  ] by  using 

B(2,  2n+2m+2)  = B(2,  2)x( .335)n+m.  (4-59) 

This  choice  of  approximation  is  illustrated  in  Figure  4-5.  From  (4-54), 
the  mean  square  of  the  bearing  estimate  becomes 


E[(i)2]  = 2kK2e”h^  B(2,  2)  Lm(h2) 

m=0 

= 2kK2B(2,  2)exp(-h2-.335k)  IQ(2h/335k) 
or  E[(i)2]  = kK2B(2,  2)exp[-(h-/335T,2]A/rl71!frh  . 


(4-60) 

(4-61) 


Substituting  for  K,  k,  and  h,  the  expression  for  the  mean  square  is 
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(4-62) 


E[(«)2]  ■ 


2.97 


42  J/2 


32»fi(N2-l ) 


(%)  exp[-NH2(D-.579)2/B']. 


These  approximations  (4-62)  and  (4-58)  are  useful  only  for 
smaller  values  of  5NR  (k  = NH2/B'  < 1). 

4.4.2  MEAN  AND  MEAN  SQUARE,  THE  MULTIPLICATIVE  ARRAY  ESTIMATE 
For  the  multiplicative  array  model,  the  high  SNR  bearing 
estimate  we  have  developed  is  given  by  (4-48),  which  can  also  be  written 

a(y)  = arcsin  K[W(y/o2+£-p)/k]^2  , 0 < y/cT+s-p  < k 

• arcsin  K,  y/a2+J--p  < 0 (4-63) 

* 0,  y/o2+J-p  > k, 

where,  in  analogy  with  the  last  section  we  here  use  K = 
k = MH2/B  = h2/D2(}Bsina;  m). 

To  obtain  the  mean  and  mean  square  of  a(y),  we  shall  use  the 
noncentral  chi-square  approximate  density  (4-34),  with  the  constant  c 
chosen  to  be  4.  The  mean  of  the  estimate  is  then 


E(o) 


-L 


o2(o4) 


-o2d/4 


dy  arcsinK  py (y) 


o2(k+p-£) 

J dy  arcsin 


oV*) 


k[i-^24-p)a)V2  Py(y). 

(4-64) 


The  first  integral  is  approximately 


4-22 


4(p-i)+d  4(p-l)+d 

/ 4y  K(a2/4)  py[a2(x-d)/4]  = | / dx  exp[-(x+X2)/2] 

0 0 

■ K - KQ[x,  \/4(p-i-)+d]»  (4-65) 


Q being  Marcum's  Q-function. 

The  second  integral,  making  the  same  approximation  as  for  the 

first,  is 

1 

K /dx  kaVl-V^ Py[cr2(kx+p-£)] 

0 

1 

*=  2Kk  / 6x\ft-Jx  exp[-(4kx+4p-2+dt'-X2)/2]I0(xAkx^4p-2+d) 
0 (4-66) 

1 _ 2 

= 4Kk  exp(-X2/2-2g)  /dx  x/V”  e"2kx  IgCZXv'kx2*^ 

0 


W CO 

.4* 


n=0  m=0 


x f dx  x2mfl/MT(l+kx2/g)n  , (4-67) 

*0 

where  we  use  the  notation  g = p - 5 + d/4.  Employing  [13],  formula 
3.259.1,  we  reduce  (4-67)  to 


4Kk  exp(-x2/2-2g)  £ £ £ (JXfjW,  2,a+2rf2) 


m=0  n=0  r=0 
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CO  00  CO 


4Kkexp(-*2/2-2g)  J 8(3/2, 2n»2r+2) 


vW/.  ,r 


m=0  n=0  r=0 


« « m 


,2\n+r 


4Kkexp(-X2/2-2g)  J ^ I {^TKT  fcijwdf8*3/2-2"*2) 

m=0  n=0  r=0 


00  00 


■ 4Kkexp{-XZ/2-2g)  £ £ 8(3/2, 2m+2)  !.Jj(X2/2),  (4-68) 

m=0  n=0 

a development  analogous  to  (4-53).  Thus  for  the  mean  value  of  the 
multiplicative  array  bearing  estimate  we  have 


E(a)  = K - KQ(X,  2/5)  + 4Kkexp(-X72-2g) 


00  00 


1 1 B(3/2-2m+2)  i>2/2)-  <4-69) 
m=0  n=0 


Going  through  a very  similar  process  we  find  for  the  mean 

square, 

E[(a)2]  = K2  - K2q(X,  2Sq)  + 4K2kexp(-X2/2-2g) 

* 1 i (9X!i(^)P  o^2>-  (4-7o) 

m=0  n=0 

The  notation  used  in  (4-69)  and  (4-70)  is  summarized  by 

K = 4.9/{W4MZ-l  k = MH2/B 

X2  = [B2+2C2+4MH2D2(&sinct;  M)(B+Ccos29)]/B2 

d * 2 + X2  ■ 4MH2D2(i3sina;  M)cos29/3-4C/B  (4-71) 

g - d/4  - .5  + C/B  - X2/4  - MH202(iBsina;  M)cos29/B. 
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m 


4.4.3 


COMPUTED  RESULTS 


Specific  cases  of  the  probability  density  functions  of  the 
estimates  (4-10)  and  (4-48)  are  given  in  Figures  4-6  and  4-7  for  several 
values  of  the  actual  bearing,  computed  according  to  (4-50).  At  first 
sight  what  seems  remarkable  about  these  curves  is  that  the  bearing  esti- 
mates at  which  the  peak  values  are  located  do  not  correspond  closely  to 
the  actual  bearing.  However,  this  in  itself  is  not  a reliable  indication 
that  the  mean  values  diverge  from  what  is  desirable  --  an  unbiased 
estimate  — since  there  is  a discrete  probability  that  the  value  of  the 
estimate  is  zero  in  each  case.  For  the  multiplicative  array  bearing  esti- 
mate, there  is  also  a discrete  probability  that  the  estimate  equals  its 
maximum  value,  arcsinK,  not  shown  in  the  figure. 

More  revealing  are  the  results  for  the  mean  values  of  the  esti- 
mates, shown  in  Figures  4-8  and  4-9  versus  actual  bearing,  with  a straight 
line  representing  an  unbiased  estimator  drawn  in  for  reference.  The 
trend,  as  the  SNR  increases,  toward  u.tbiasedness  is  clearly  seen.  This 
is  comforting,  since  the  theory  of  maximum  likelihood  estimation  predicts 
this  pehnomenon.  These  computations,  based  on  (4-53)  and  (4-69),  required 
on  the  order  of  50  terms  of  each  summation  for  H2  = 1,2  and  100  terms  for 
V-  = 4. 

The  curves  drawn  with  dashed  lines  in  Figure  4-9  represent  the 
series  term  in  (4-69),  so  that  the  importance  of  the  Q term  is  seen. 

It  appears  from  our  limited  computations  that  the  standard  and 
the  multiplicative  array  estimators  on  the  whole  perform  equally  well 
with  respect  to  bias.  A kind  of  "end  effect"  is  observed  in  either  case 
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FIGURE  4-6 
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FIGURE  4-7 


FIGURE  4-8 


FIGURE  4-9 


when  a approaches  its  maximum  value.  The  difference  is  that  the 
standard  array  estimator  behaves  smoothly  in  this  regard,  while  the 
multiplicative  estimator  displays  a more  drastic  reaction  to  this 
limitation,  imposed  by  the  formulation  (4-63)  of  the  estimator.  It 
should  be  remembered  that,  if  our  attention  is  mainly  upon  bearings 
within  the  major  lobe  (a  < arcs  in  .2  for  the  standard,  <*  < arcsin  .1 
for  the  multiplicative),  then  the  important-  region  is  that  near  a = 0. 

Computations  or  mean  square,  variance,  and  mean  square  error 
are  shown  in  Figures  4-10  and  4-11.  i.^m  the  curves  it  is  evident  that 
with  respect  to  tr.ese  measures,  the  stancard  array  estimator  performs 
better  than  the  multiplicative  for  increasing  bearing. 

Further  consideration  of  estimator  variance  and  mean  square 
error  is  given  in  the  next  chapter,  where  we  compare  them  with  theoreti- 
cal bounds. 
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FIGURE  4-11  ESTIMATOR  MEAN  SQUARES  AND 
VARIANCES,  MULTIPLICATIVE  ARRAY 
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CHAPTER  FIVE 


BEARING  ESTIMATION  PERFORMANCE 


As  mentioned  in  the  introduction  to  Chapter  4,  in  addition 
to  characterizing  the  detection  performance  of  the  two  array  models  we 
are  studying,  we  also  desire  to  evaluate  the  bearing  estimation  capa- 
bilities of  these  configurations.  That  is,  we  should  like  to  know  with 
what  potential  precision  the  filter  outputs  of  these  array  processing 
models  can  be  operated  upon  to  yield  a value  for  the  bearing  of  the 
signal  source. 

In  the  previous  chapter  we  derived  several  forms  for  maximum 
likelihood  bearing  estimators,  based  on  the  probability  density  functions 
at  the  filter  outputs  which  were  found  in  Chapter  2.  In  this  chapter 
we  shall  examine  the  performance  of  these  estimators  in  terms  of  the 
relationship  between  the  theoretical  minimum  estimator  variance  and  the 
actual  variance  of  the  given  estimator,  as  a function  of  SNR. 

Although  there  are  other,  more  exact  bounds  [34,38],  we  shall 
be  making  use  of  the  Cramer-Rao  bound  on  the  variance  of  an  estimator: 


Var(a)  *[^WI‘/E 


3lnp(W;  g) 


3a 


(5-1) 


where  p(w;  a)  represents  the  distribution  cu  which  the  estimate  of  the 
parameter  is  based  [16],  and  where  certain  regularity  conditions  on  the 
probability  density  [39]  hold. 


/ 


Bounds  will  first  be  calculated  based  on  the  Gaussian  distri- 


butions seen  by  the  array  elements,  and  also  those  gaussian  distributions 
holding  at  the  array  sums.  Then  hounds  based  on  the  filter  output  distri- 
butions will  be  calculated.  These  latter  bounds  we  shall  refer  to  in 
a general  way  as  noncentral  chi-square  bounds. 

5.1  CRAMER-RAO  BOUNDS  FOR  GAUSSIAN  DISTRIBUTIONS 

If,  analogous  to  Section  2.1.1,  we  annotate  the  N samples 
{wk>  of  a signai-plus-noise  process  by  the  vector  w and  the  corres- 
ponding signal  components  { x^}  by  x_,  the  joint  distribution  for  the 
Gaussian  case  is 

P(w;  a)  = — L— — exp  [-Hw-xJ'K-1  (w-x)],  (5-2) 

(2ir)  /detK^ 

where  a is  the  parameter  of  interest  and  = E[ (w-x) (w-x)']  is  the 
covariance  matrix  of  the  samples,  which  is  symmetric.  For  this  expres- 
sion of  the  distribution,  we  have 

,nd  psM' 2 - (H)  k-i  (w-xKw-xj-K-yii) 


(5-3) 

(5-4) 


so  that  t^e  expectation  of  (5-4)  is  given  by 


(S)kn1eC(--)(-"^,] 


(5-5) 


5-2 


A -iirflitMj 


The  form  (5-10)  is  equivalent  to  that  shown  by  Seidman  [34]  in  the  same 
case.  The  error  of  a bearing  estimate  a(w;  a)  based  on  the  array  input 
noise  distribution  for  white  noise  is  bounded  by 


Var  (a)  > 6[|^(a)]2/N(N2-1)HVcos2a. 


1 5-12) 


5.1.2  APPLICATION  TO  THE  ARRAY  SUMS 

For  the  standard  array  model,  which  has  the  single  array  sum 
(1-4),  the  quadrature  components  are 


v - n + NAD(£gsina;  N)cos[£{N-l )3sincc] 
c c 


w = n + NAD(igsina;  N)sin[i(N-l )gsina], 


where  n and  n are  modeled  as  independent  Gaussian  variates  with 

2 

variances  cr2  = iiB's  (Section  2.2.2).  Thus  (5-5)  in  this  case  is 


(5-13) 


N2aVcos2cx 

4o2 


| j D'cos[i(N-l)gsina]  - (N-l )Dsin[i(N-l )gsina]J' 


+ |D'sin[-f(N-l)gsina]  + (N-l  )Dcos[£(N-l)gsina]  |2] 

(5-14) 


_NH^cos2_a  [(d,}2  + (N.1)2D2]> 


so  that  an  estimate  a(wc>ws;  a)  based  on  the  standard  array  sum  has  the 
bound 


Vas  (a)  > 28 1 [|^£(a)]2/HH262cos2a[(D' )2  + (N-1)2D2]. 


(5-15) 


For  the  two  array  sums  of  the  multiplicative  array  model,  the 
quadrature  roise  components  have  the  covariance  matrix  given  by  (1-7). 
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Also,  from  (1-5)  we  have 


i 

. 


xlc  a Scose^  x1$  = Ssine^ 
x2c  e Scose2,  x2s  = Ssine2 
where  S = MAD(i£sina',  M) 

6]  = -HM-l)$sina,  02  = -H3M-l)3sina 


(5-16) 


(5-17) 


For  what  we  have  called  our  main  case,  the  quadratic  form  (5-5)  is 


o2(l-p2) 


[/3xlc\2  + 

/9xls  f t 

/3x2c  f + 

/8x2s  \2 

l 8a  / 

\ 3c  / 

\ 3a  / 

l 3a  / 

■*!&)(£)  ♦&-)(£)] 


(6-18) 


M2A232cos2g 

4a2(l-p2) 


[D1  cosO-j -(M-l  jDsinB^]2  + [D'sin9^+(M-1  jDcose^2 


+ [D'cos62-(3M-l)Dsin02]2  + [D'sin92+(3M-l)Dcos62]2 


- 2p[D'cos0-j-(M-l)Dsin6^][D'cose2-(3M-l)Dsine2] 

- 2p[D'sin9^+(M-l )Dcos0^][D'sin92+(3M-l )Dcos92] 

m!&osjx  I 2(D')2+[(M-1)2+(3M-1)2]D2 
2B(l-p2)  1 

-2pcos29[(D' )2+(M-l ) (3M-1 )D2] -4pMDD ' s i n20 I , 


(5-19) 


(5-20) 


where  20  = 0-|  - e2-  If  the  phase  reference  is  made  the  center  of  the 
array,  that  is  if 
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.t.i_my.)ijipiyn-_  , M i,)y^y^rT-ry-  .T-"  . i?jq  vrryfTTHxr- 


0!  = -e2  = £(2M-l)psina, 


(5-21) 


then  (5-5)  becomes 


2 2 2 

m P c0!  n 2[(D')2+(2M-l)2D2](l-pcos29) 

2B(l-p2) 


2 2 2 

MH.i_cos_a_  £(D')2+(2M-l)2D2](B-CcosM3sina).  (5-22) 

B -t 


Therefore,  an  estimate  a(w,  a)  based  on  the  two  array  sums  has  an  error 
bound, 

Varfa)  * (B2-C2)[|^E(a)]2/MH2B2cos2a[(D,)2+(2M'l)2D^](B-CcosM3sina) . 

(5-23) 

5.1.3  COMPARISON  OF  BOUNDS  AT  ARRAY  SUMS 

For  unbiased  estimates,  that  is  for  E(a)  = a,  we  may  compare 
the  minimum  estimator  variances  of  the  two  array  configurations  by  com- 
paring (5-15)  and  (5-23).  Let  the  minimum  variances  of  the  standard 


md  multiplicative  arrays  be  denoted  a2.  and  a2.  respectively. 

min,s  min,m  r J 


Then,  for  the  same  number  of  elements  (N=2M), 


min.s  _ v N/2 


(D;„)2+(N-1)2D2  , D2 


N/2  * N/2 


- (K)2  * (N-1  )2  D2  D2 


(5-24) 


min.m  v N' 


N 


'N 


where  Dj,  = DQ  3 sin  a;  K).  Since  for  the  angles  of  interest, 

we  have 


_2  2 
^mln.s  > min.m 


(5-25) 


at  the  array  sums,  with  equality  holding  at  zero  bearing.  For  example, 
for  N = 10,  8 = ir,  and  a = .05,  the  ratio  is  .948/. 808  = 1.24. 
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5.2  CRAMER-RaO  BOUND  FOR  NONCENTRAL  CHI-SQUARE  DISTRIBUTIONS 
Let  x be  a noncentral  chi-square  random  variable  with  two 
degrees  of  freedom  and  noncentrality  parameter  X2.  Then  the  probability 
deny.'ty  function  of  x is 

p(x)  * (l/2)exp[-(x+A2)/2]  I 0(x<5).  (5-26) 


where  X = x(a),  a a parameter  to  be  estimated.  For  this  density  func- 
tion, 

MWx) 


3lnp(x;  a) 
3a 


-XV  + 


T 

i0(WSO 


X'/x 


(5-27) 


so  that  E 


(*J?)‘ 


of  o [ l-i  (X»^() 

= (X')2(x2  - 2XE.V5T  1 


I0(X^ 


+ E 


IT(X^x) 


12 


li0(*£) 


(A')? 


A2  - <f"2 


00 

/dxtf  e"X/^2I-j  (X»^x) 


+ (l/2)e"^/2  /xdx  e'x/2  i (X^) 


I„(A/x) 


(5-28) 


From  Gradshteyn  and  Ryzhik  ‘[13],  formulas  6.643.2,  9.220.2)  and 
Middleton  ([16],  formula  A. 1.19b)  we  have  'or  the  first  integral. 


/*> 

dxv^T  e”x/2  IT(XA)  = 2X1F1(1,2;  X2/2)  = - -J-  (l-eW2).  (5-29) 

0 
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For  the  second  integral  we  h?ve  approximately, 

/dx  x e"x/2  -1-  I,(X./x)  = f dx  x e"x^2  I,(X»fi) 

•4  !«(**)  1 -4 


(5-30) 


|XVST  1F1 (5/2,2;  X2/2) 


■ X^5T[^1F1(1/2.2;  X2/2)  + (1+X2/2)1F1  (3/2 ,2 ; X2/2)] 

2 

*•  £ v$T  eX  /4[(3+X2)I0(x2/4)  + (1+X2) I-, (X2/4)] . (5-31) 


aguin  making  use  of  [13]  and  [16]. 

Therefore,  (5-28)  becomes 


/31np\2 

1157 


. (x')2|x2  * 4e-x2/2  0-//2> 


* i JS  e'X  V4[(3+X2)l0(X2/4)  + (1+X2)Iv(X2/4)]|. 


(5-32) 


for  large  arguments  of  the  Bessel  functions  (or  large  \)  we  may  write. 


(X')2(2X2-2+4e"A  /2)  = 2(X' )2(X2-1 ) . 


(5-33) 


5.2.1  CRAMER-RAO  BOUND  FOR  THE  STANDARD  ARRAY  FILTER  OUTPUT 

For  the  standard  array  filter  output,  which  is  distributed 
as  a noncentral  chi-square  variable,  we  have 
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(5-34) 


A2  ■=  2h2  ■=  2NH2D2(}gs1n»;  N)/B‘. 


Thus  the  denominator  of  the  Cramer-Rao  bound  expression  (5-1)  is 

21 


= 4(h')2(2h2-l) 


= 2N3H232cos2a(D')2(2NH2D2-B')/(B?)2. 
From  (4-53)  we  may  write,  using  [13],  formula  8.971, 


3_ 

3a 


E(a)  =-2Kk2DD'NBcosae'kE)£  B(3/2,2m+2)[Ln+Li[)_1] 

m=0 


= -2Kk2DD'Bcosae‘kDN  B(3/2,2m+2)  L^kD2), 

m=0 


2 ? 

since  h =kD  , so  that 


[|^  E(i)]2  = 4K2k4D2{D')2N2B2cos2a  e'2kD 


2 B ( 3/2 , 2m+2 ) ^(kD2) 


m=0 


The  Cramer-Rao  bound  for  the  standard  array  filter  output 
is  then  the  ratio  of  (5-37)  to  (5-35),  using  i to  denote  the  series 

Var(a)  > 2K2N3H6D2(  I }2e-2ND2H2/B'  /(B • )2(2NH2D2-B * ) 


v 24N2H4(  £ )V2Nd2h2/B'  /(B')282(N2-1). 


(5-35) 


(5-36) 


(5-37) 


s 


(5-38) 


For  white  noise  (B 1 =1)  this  expression  becomes 

Var(a)  » 24N2H4(  2 w)2e"2ND2n2  /p2(N2-l),  (5-39) 

and  when  the  bearing  is  zero  also, 

Var(i)  > 24N2H4(I  w 0)V2'Ih2 /p2(N2-l).  (5-40) 

5.2.2  CRAMER-RAO  BOUND  FOR  THE  MULTIPLICATIVE  ARRAY  FILTER  OUTPUT 
For  the  multiplicative  array  filter  output,  which  is  distri- 
buted approximately  as  a noncentral  chi-square  variable,  we  have  from 
(4-71)  and  Section  4.3.2, 

X2  = 1 + 2p2  + 4h2(l+pcos23) 

= [B2  + 2C2  + 4MH2U2(&sincr,  M)(B+CcosM3sk  i)]/B2.  (5-41) 

The  "x"  in  (5-26)  corresponds  here  to  4y/'o2+d,  where  y is  the  filter 
output.  The  expression  (5-33)  may  be  used  here  if  we  can  say  that 
3x/3a  = v>,  as  assumed  in  (5-27)  and  in  the  following  analysis.  But 

I?  * M = ST  [3  • 4P  + 2P2  + ^ - 4h2(l-p)cos20] 

* 4M2H2D3cosa  {D‘ [B-(B-C)cos20]  + D(B-C)sin2e)/B2.  (5-42) 

For  white  noise  (B  = 1,  C = 0)  i >is  relation  is 

|^-  = 8M2H2D3cosasin6(D'sin9+Dcos0) > 0 as  a > 0.  (5-43) 

Thus,  though  we  shall  use  (5-33),  we  do  so  with  the  understanding  that 
the  solution  involves  this  additional  approximation  when  the  bearing 
is  nonzero. 
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Proceeding  to  apply  (5-33),  then,  we  have 


= 2(A')2(A2-1)  = 2[(A2)7.2A]2(A2-1) 

e [4M2H2Dgcosa ( D 1 B+D  ‘ C cos 2 9 -DCs i n29 ) 1 2 
B2  + 2C2  + 4MH“D2 ( B+Ccos29 ) 

n [2C2  + 4MH2D2(B+Ccos29)]/2B2. 


(5-44) 


For  white  noise, 

= 32M^ti6D4g2cos2a(Dl  )2 
1 + 4MH2D2 

Using  (4-69),  we  may  write 


|^E(i)  = -KA'|rQ(A,  2,^g)  - 4Kkexp(-A2/2-2a)  AX' 


(5-45) 


2 


Y 

L.t 

m-0  n=0 


(-k  r 


B(3/2,2m+2)  L 


n+1 


(^>n- 

(5-46) 


where  g = d/4  + C/B  - .5.  Representing  the  series  by  ee,  we  have 

• - K2(A')2[|^q(A,  2/g)  + 4kAe,xp(-A2/2-2g)ZZ]2,  (5-47) 

and  the  Cramer-Rao  oound  for  the  multiplicative  array  filter  output 
may  be  expressed,  using  (5-44), 


Va r(a)  > [fxQ(^»  2/g)  + 4kAexp(-\2/2-2g)Zz]2,  (5-48) 
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where  the  notation  is  given  by  (4-71).  For  small  angles,  the  Q term 
vanishes,  yielding 

Var(a)  £ 8K2kVexp(-A2-4g)(ZZ)2/(X2-l).  (5-49) 

For  white  noise  we  obtain  for  small  angles, 

Var(a)  * exp[-4-4MH2D2(2-cos29)](EE)2.  (5-50) 

eV4M"-l  1+4MH  F 

5.3  COMPUTED  RESULTS 

Computations  of  the  various  bounds  on  estimator  error  derived 
in  this  chapter  are  shown  in  the  accompanying  figures. 

In  Figure  5-1,  we  show  diagrammatical ly  how  these  bounds 
correspond  to  the  density  functions  at  different  points  in  the  array 
models.  The  bound  at  the  array  inputs  (5-12)  is  based  on  the  joint 
probability  density  function  of  the  outputs  of  the  array  elements,  and 
take  into  account  the  basic  geometry  of  the  array.  The  bounds  (5-15) 
and  (5-23)  are  based  on  the  density  functions  of  the  array  sums,  and 
reflect  knowledge  of  the  directivity  gained  by  the  summing.  And,  the 
bounds  (5-38)  and  (5-48)  are  baser  on  the  filter  otuput  probability 
density  functions,  and  incorporate  the  effects  of  the  nonlinear  processing 
involved  in  detection. 

This  procedure  of  computing  bounds  at  several  points  in  a 
system  is  somewhat  unusual,  but  was  undertaken  for  two  reasons:  (1)  after 

the  summations,  to  provide  a basis  of  comparison  where  the  two  array  systems 
begin  to  differ;  (2)  after  the  filters,  partly  out  of  curiosity  and 
partly  to  obtain  a measure  of  the  effects  of  the  nonlinear  processing. 

The  effort  was  worthwhile  if  only  for  the  interesting  discussions  the 
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results  have  stimulated  thus  far,  points  which  we  shall  mention  in  the 
following  presentation  of  the  numerical  results. 

Figure  5-2  displays  the  bounds  (5-15)  and  (5-23)  in  comparison 
with  (5-12)  for  estimators  assumed  to  be  unbiased,  that  is  for  which 
E(a)  = a,  for  H = 1 and  H =2,  versus  actual  bearing.  The  arrays  are 
assumed  to  have  ten  elements  and  half-wavelength  spacing.  There  are  two 
significant  features  of  this  figure.  The  first  is  that  for  small  bearings 
the  bounds  computed  from  the  uensity  functions  at  the  array  sums  are 
lower  than  that  computed  at  the  array  inputs.  In  fact,  if  we  compare 
(5-15)  and  (5-12)  for  zero  bearing  we  find  that 

Bound  at  sums  _ 1^  N+l_  n ,r  ci\ 

Bound  at  inputs  " 3 N*1  * “ ’ ' 


where  N is  the  number  of  elements.  Thus  for  N = 2 we  see  that  the  bounds 
are  indeed  equivalent,  so  that  for  other  values  of  N the  comparison 
should  be  valid.  On  the  other  hand,  given  the  manner  in  which  the  bear- 
ing is  "encoded"  or  embedded  in  the  signal  terms  and  recalling  the  results 
obtained  for  estimator  means  in  Chapter  4,  we  cannot  push  this  comparison 
ve^y  far  since  it  appears  likely  that  estimates  of  bearing  based  on  the 
summations  will  in  general  be  biased. 

The  second  feature  to  note  in  Figure  5-2  is  that  the  bound 
(5-23)  for  the  multiplicative  array  bearing  estimate  lies  below  that  for 
the  square-law  or  standard  array  (5-15).  This  has  already  been  noted  in 
Section  5.1.3,  and  implies  that  the  accuracy  of  the  multiplicative 
configuration's  two  summation  (split-beam)  approach  is  potentially 
better  for  bearing  estimation.  Also,  the  increasing  error  for  higher 
angles  typifies  the  cost  of  beamforming. 
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Further  comparisons  among  these  bounds  are  illustrated  in 


! > 


Figure  5-3,  versus  SNR.  Again,  the  position  of  the  bound  corresponding 
to  the  array  sums  of  the  multiplicative  array  model  relative  to  that  of 
the  standard  model  would  seem  tc  imply  that  the  multiplicative  array 
promises  better  performance  for  nonzero  values  of  bearing,  as  noted  in 
Section  5.1.3. 

The  somewhat  difficult  to  compute  noncentral  chi-square  bounds, 
(5-38)  for  the  standard  array  estimator  and  (5-48)  for  the  multipli- 
cative, are  shown  versus  true  bearing  in  Figure  5-4.  Since  estimator 
bias  is  included  in  these  bounds,  these  curves  seem  to  be  informing  us 
that  the  combination  of  the  nonlinear  processing  and  the  particular 
estimator  forms  we  have  selected  promise  much  improved  performance  over 
that  achievable  at  the  array  sums,  for  this  range  of  bearings.  However, 
it  may  be  also,  as  Seidman  implies  in  similar  cases  [34],  that  the 
Cramer-fao  bounding  technique  yields  a loose  (too  small)  bound  for  the 
error  over  the  range  of  bearings  and  SNR  in  which  we  are  interested. 

Estimator  error  curves  from  Figures  4-10  and  4-11  are  super- 
imposed in  Figure  5-5  versus  true  bearing.  From  their  shape,  it  is 
obvious  that,  these  curves  are  "close  relatives"  to  the  post-filter 
bounds  in  Figure  5-4.  However,  it  is  also  manifest  that  these  esti- 
mators do  not  achieve  tne  error  performance  implied  by  the  post-filter 
bounds.  Rather,  the  actual  error  curves  are  more  closely  bounded  by 
the  bounds  of  Figure  5-2.  This  occurrence,  we  believe,  indicates  the 
looseness  of  the  post-filter  bounds. 
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Judging  by  the  manner  in  which  the  actual  estinator  errors 
behave  with  increased  SNR  (decreasing  not  nearly  as  "fast"  as  the 
bounds),  we  may  observe  that  the  particular  forms  for  the  estimators, 
though  obtained  by  means  of  the  maximum  likelih^u  method,  are  not 
asymptotically  efficiert.  That  is,  from  the  data  we  have  obtained, 
the  error,  though  decreasing  with  higher  SNR,  does  not  converge  toward 
the  lower  bound.  Again,  we  remark  that  this  may  be  due  to  the  limita- 
tions of  the  bounding  technique,  but  quite  likely  also  the  estimator 
forms,  being  approximations,  are  therefore  suboptimum. 


CHAPTER  SIX 


CONCLUSIONS 


With  this  chapter  we  conclude  this  work  and  take  stock  of  what 
has  been  done.  First,  the  work  is  summarized,  with  particular  emphasis 
upon  what  are  considered  to  be  the  primary  results.  These  comments  are 
then  qualified  or  put  into  perspective  by  a more  critical  review  of  the 
assumptions  and  methods  employed.  Finally,  suggestions  for  improvements 
and  extensions  of  the  work  are  given. 


6.1 


SUMMARY 


We  began  by  describing  a "debate"  in  the  literature  over  the 
relative  merits  of  what  were  termed  the  standard  and  the  multiplicative 
processing  of  signals  intercepted  by  linear  arrays  of  receiving  elements 
(hydrophones,  for  definiteness).  Our  conception  of  this  debate  was  to 
that  a considerable  degree  the  issues  are  clouded  by  insufficient  charac- 
terization of  the  random  processes  involved,  in  that  some  of  the  more 
important  performance  measures  by  which  the  two  configurations  ought  to 
be  compared  are  probabilistic  in  nature  (e.g.,  probability  of  detection). 
The  goal  of  this  work,  therefore,  was  set  to  provide  at  least  a partial 
remedy  to  this  situation  by  applying  some  of  the  analytical  techniques  of 
statistical  communication  theory  to  the  problem,  with  a view  toward  system 
comparison. 

Because  a kind  of  "gain/bandwidth"  tradeoff  exists  for  analyses, 
rather  simple  narrowband  models  of  standard,  or  square  law,  and  multipli- 
cative array  processing  systems  were  chosen  so  that  attention  could 
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be  given  to  as  rigorous  an  analysis  of  the  performance  of  these  systems 
as  our  abilities  would  permit.  In  this  fashion  we  were  enabled  to  include 
as  generalizations  the  various  array  parameters,  such  as  the  number  of 
elements,  the  inter-element  spacing,  and  the  noise  covariances. 

A rather  general  probability  density  function  for  the  filter  out- 
put of  the  narrowband  array  processor  was  derived  by  the  direct  method. 

This,  we  believe,  is  our  basic  result  in  the  sense  that  the  remainder  of 
the  work  consists  of  application.  The  degree  of  generality  is  such  that 
the  effects  of  all  the  model  parameters  can  be  studied.  Moreover,  alternate 
computational  forms  were  given  for  the  main  case,  in  which  the  array  was 
considered  symmetric  and  the  noise  assumed  to  have  a spectrum  even  about 
the  bandpass  center  frequency. 

Theorizing  a typical  Neyman-Pearson  characterization  of  signal 
detection  at  the  filter  outputs  of  the  models,  we  were  able  to  use  the 
probability  densities  to  calculate  detection  and  false  alarm  probabilities 
and  receiver  operating  characteristics.  From  the  numerical  results 
obtained,  we  observe  that  the  standard  detector  or  processor  is  only  very 
slightly  better  than  the  multiplicative  in  terms  of  signal  detection,  so 
that  we  would  assume  that  the  better  beamforming  property  of  the  multipli- 
cative array  gives  it  the  edge  in  this  category. 

The  more  difficult  area  of  application  of  the  probability  density 
functions  was  in  exploring  the  bearing  estimation  capabilities  of  the 
arrays.  We  fend  ourselves  unable  to  solve  the  general  maximum  likelihood 
equation  for  other  than  very  small  SNR.  However,  our  exact  knowledge  of 
the  distribution  of  the  filter  output  made  it  possible  to  select  as  a good 
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approximation  the  somewhat  more  tractable  noncentral  chi-square  density 
function,  with  a linear  function  of  the  multiplicative  array  filter  out- 
put as  its  argument.  Forms  of  maximum  likelihood  estimators  of  bearing 
were  shown  for  both  high  and  low  SNR  cases  for  both  array  models. 

Further  treatment  was  given  to  these  estimators  for  the  high 
SNR  cases,  in  the  form  of  showing  their  means  and  mean  squares,  as  well 
as  in  displaying  example  computations  of  their  probability  distributions. 
It  was  seen  that  these  estimators  do  approach  unbiasedness,  although  our 
modeling  of  the  multiplicative  array  estimator  apparently  had  a negative 
effect  on  its  performance  in  this  respect.  Cramer-Rao  bounds  on  the 
minimum  estimator  error  were  calculated  for  various  points  in  the  array 
systems,  with  the  bound  based  on  the  density  function  at  the  array  sums 
coming  the  closest  to  the  actual  mean  square  errors  of  the  estimators. 
Although  with  respect  to  mean  square  error,  the  standard  array  estimator 
appeared  to  perform  slightly  better,  it  appeared  also  that  for  some 
ranges  of  bearing  the  multiplicative  crray  estimator  is  better.  Neither 
of  the  two  specific  forms  for  the  estimators  we  selected  was  very 
efficient,  however,  an  occurrence  which  v/e  attribute  to  the  amount  of 
approximation  involved  in  their  derivation. 

6.2  CRITIQUE 

In  keeping  with  our  goal  to  shed  light  on  the  discussion  of 
array  performance,  we  now  offer  some  critical  comments,  pro  and  con,  on 
the  information  which  we  have  generated  to  assist  the  reader  in  assessing 
it. 

We  believe  that  the  derivation  of  the  probability  density 
functions  we  have  shown  represents  an  important  contribution  to  the  know- 
ledge in  this  area,  to  the  extent  that  the  popular  approach  on  which  it 
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Is  based  (utilizing  rather  idealistic  conceptions  of  lowpass  filtering) 
is  faithful  to  the  behavior  of  narrowband  systems.  We  should  like  to 
remind  the  reader  that  several  of  the  parameters  on  which  the  distribution 
is  conditioned,  such  as  bearing  and  signal  amplitude,  might  more  realis- 
tically be  treated  as  random  variables.  Also,  the  validity  of  assuming 

r 

} 

stationary  Gaussian  noise  processes,  although  a popular  procedure,  is 

| 

sometimes  questionable  for  underwater  applications. 

To  our  knowledge,  the  detection  analysis  which  we  have  offered 
Is,  in  its  exactness,  also  a "step  input"  to  the  narrow  field  for  which 
it  applies,  although  we  would  not  be  surprised  to  discover  that  similar 
data  exists  in  various  classified  memoranda. 

} 

Most  comparative  detection  analyses  which  have  come  to  our 
attention,  besides  those  cited  in  Chapter  1,  rely  upon  detector  output 
SNR  as  a figure  of  ..ierit  (e  g.,  [58-60]).  While  probability  of  detection 
is  indeed  proportional  to  SNR,  we  have  demonstrated  that  comparison  on 
the  basis  of  SNR  may  exaggerate  the  advantage  of  the  system  having  the 
greater  SNR.  However,  this  is  not  to  discount  this  approach  where  the 
comparison  involves  detectors  with  the  same  type  of  output  probability 
density  function,  when  this  property  has  been  shown,  as  in  optimization 
of  a given  configuration  [61,  62]. 

The  amount  of  data  which  we  have  shown  is  small,  but  we  offer 
our  computer  programs  in  Appendix  D as  a cover  for  this  deficiency,  if  it 
be  such. 

The  numerical  results  of  our  investigation  of  maximum  likelihood 
bearing  estimation  from  the  filter  outputs  are  somewhat  sparse,  in  pro- 
portion to  the  amount  of  labor  spent  in  derivation.  However,  we  are 
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satisfied  that  the  several  theoretical  expressions  we  have  shown  are 
useful  accomplishments  in  themselves,  going  by  the  lack  of  assistance 
offered  by  the  literature  in  this  area.  It  is  worth  noting  that  though 
the  feasibility  of  optimally  estimating  bearing  using  the  filter  outputs 
has  been  demonstrated,  a more  reasonable  approach  perhaps  is  to  recon- 
sider the  entire  configuration  for  these  purposes.  Stremler  and  Brown 
[54],  for  example,  show  a bearing  estimating  system  using  a phase  detector. 
Also,  split-beam  and  monopulse  techniques  are  popular  (e.g.,  [53,  58]). 

6.3  SUGGESTIONS  FOR  FURTHER  WORK 

There  are  several  directions  which  improvements  and  extensions 
of  this  work  may  take.  We  suggest  a few  examples: 

--  Direct  application:  using  the  PDF's  and  detection  analysis 

to  evaluate  similar  systems,  such  as  the  split-beam  tracker,  for  which 
the  analysis  either  applies  or  can  be  modified  to  apply. 

--  Further  generalization:  expanding  the  classes  of  array 

geometries;  treating  ncndeterministic  or  unknown  signals;  considering 
joint  estimation  of  bearing,  SNR  (see  [55]),  and  perhaps  range;  treating 
nonstationary/nongaussian  noise  processes. 

--  Improvements:  using  numerical  methods  to  calculate  more 

precise  error  bounds  and  to  simulate  more  exact  forms  of  the  maximum 
likelihood  estimator  (see  [56]);  additional  theoretical  investigations 
to  specify  the  system  comparison  problem  over  the  broader  range  of 
bearings;  a closer  examination  of  the  small  SNR  case. 


APPENDIX  A 


DISCUSSION  OF  THE  ARRAY  INPUT  NOISE  PROCESS 


Let  the  array  input  noise  process  be  denoted  v(x,t)  before 
steering  and  n(x,t)  after  sterring,  where  x represents  distance  along 
the  array,  with  x=0  at  the  center  of  the  2M. elements.  Let  v(x,t)  be 
Gaussian,  zero-mean,  and  stationary;  consequently  so  is  n(x,t). 

The  steering’s  effect  of  the  noise  "sample"  at  the  k: th  hydro- 
phone shall  be  characterized  by  the  relationship 

nk(t)  = vk[t  + (k-M-i)T],  k = 1,2 2M  (A-l) 

A.l  AN  ASSUMPTION 

If  R(t)  is  t'.ie  normalized  autocorrelation  function  of  the 
noise  process  before  steering,  let  us  assume  that  we  may  combine  time 
and  space  correlation  along  the  array  such  that  the  covariance  matrix 
of  the  noise  after  steering  is  given  by 


E[nk(t)nm(t)]  = E[vk[t  + (k-M-i)T]vm[t  + (m-M-i)T] 


= Rkn|[(m'k)T]  = VksmRC(nvk)T], 


(A-2) 


where  the  R^t)  are  the  interelement  crosscorrelation  functions,  the 
{rkm}  are  the  "spatial  correlation"  coefficients,  and  the  { s^ } are  the 
noise  variances  at  the  individual  hydrophones. 

In  other  words,  we  are  supposing  that  the  effects  of  steering 
on  the  noise  process  covariances  can  be  modeled  by  the  "cascading"  of 
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spatial  effects  (present  for  no  steering,  or  T=0)  and  time  delay 
effects,  represented  by  the  autocorrelation  function  of  the  process. 
A. 1.1  IMPLICATIONS  FOR  THE  GENERAL  CASE 

The  covariances  of  the  multiplicative  array  sums  are  given 
by 


M M 

E(n?)  = a?  = l l r.  s.  s R](m-k)T]  (A-3) 

1 1 k=l  m=l  m K m 

(the  sum  of  the  'upper  left  quarter1  of  the  array  noise  covariance 
matrix) , 


M 2M 

E(n,n2)  = p0,c2  ■ I l VksmR[(™-k)T]  (A-4) 


(the  sum  of  the  'upper  right  quarter1  or  'lower  left  quarter1  of  the 
array  noise  covariance  matrix),  and 


E(n|) 


2M 


l 

k=M+l 


2M 


I 

m=M+l 


rkmsksn,R[(m-k)T] 


(A-5) 


A. 1.2  IMPLICATIONS  FOR  THE  NARROWBAND  CASE 

In  the  narrowband  case  we  write  vk(t)  = vck(t)  cos  + 
vsk(t)  sin  ut.  Davenport  and  Root  [12]  show  that  we  may  write 


E[vck(t)vck(t+x)]  = E[vsk(t)vsk(t+T)] 


dfW  (f)  COS  2irfi 
v 


= S 2 Rc(t),  and 


A- 2 


(A-6) 


A . 2 EFFECT  OF  NOISE  PROPK  TIES  ON  ARRAY  SUM  COVARIANCES 

The  following  is  simply  a compilation  of  certain  noise  proper- 
ties, so  far  not  specified,  and  their  consequences  for  the  array  sum 
noise  parameters  we  have  just  related  to  the  input  process:  o-j , o ^ p»  r- 
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Narrowband  noise  with  a spectrum  even  about  the  center  frequency 


For  this  case  R (x)  = 0,  and  R(x)  = R (x)  cos  tux,  so  that  r=0. 

v b V 

Uniform  variance  across  array,  correlation  proportional  to  distance 
Here  s^  = s,  all  k,  and  r^  = f(|m-k|),  so  that  = c^. 

Whi te  noise 

We  have  R ( x ) = 0,  R (x)  = <s(x),  s^  +hat  p = r = 0. 

v J v> 

Spatially  independent  noise 

If  rta = w then » = r = °- 

(Note  the  duality  between  spatial  properties  and  spectral  properties.) 

A. 3 FURTHER  MODELING  OF  THE  SYMMETRIC  NOISE  COVARIANCE  MATRIX 

When  the  noise  across  the  array  is  symmetric  and  either  white 
noise  or  spatially  independent,  we  have  o|  = o|  = o2  = Ms2  and  p = 0. 

But  we  should  like  to  be  a bit  more  general  than  that.  The  following 
artifices  are  developed  to  provide  an  alternative  to  the  white  noise 
assumption. 

A. 3.1  ADJACENT  CORRELATION  ONLY 

If  we  say  that  only  the  adjacent  noise  inputs  are  correlated 
with  one  another,  with  correlation  coefficient  a,  it  is  easy  to  perform 
the  summations  of  A. 1.1  to  find 

a2  = Ms2  + 2(M”1 )as2 , pa2  = as2,  a<l.  (A-13) 

A. 3. 2 EXPONENTIALLY  DECREASING  CORRELATION 

If  we  say  that  the  correlation  coefficient  between  hydrophones 
m and  n is  an  exponentially  decreasing  function  of  the  distance  between 
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them,  that  is  rmnR[(m-n)T]  = e“M,n“nlj  we  have  probably  a more  useful 
representation.  For  small  numbers  of  hydrophones  the  calculation  is 
simple.  For  examples: 

M = 1:  a2  = s , = e’b 

M = 2:  a2  = 2s2(l  + e'b),  pa2  = s2(e'b  + 2e'2b  + e"3b) 

M = 3:  a2  = s2(3  + 4e"b  + 2e’2b), 


pa2  = s2(e"b  + 2e'2b  + 


3e~3b  + 


2e 


-4b 


+e-5b). 


It  is  apparent  that  for  large  M or  for  a general  expression,  unless 

terms  are  discarded  (effectively  what  is  done  in  the  tridiagonal  case) 

then  our  effort  to  achieve  a closed  expression  is  defeated.  A way 
around  this  abrade  is  to  model  the  sums  of  offuiagonal  terms  as  inte- 
grals, asserting  that  the  sum  is  proportional  to  area  of  an  appropriately 
chosen  region  under  the  curve  exp(-b| x-y| ) . In  this  fashion  we  have 

a2  = Ms2  + 2M(M-1  )s2  ^ dx  J*  dy  e'b^x'y^  = Ms2  + 2K(M-l)s2 

(A— 1 4 ) 

pa2  = M2s2  J dx  J dy  e b^x_y^  _ M2s2e"b  ^ ^ . (A-l 5 ) 

Thus  for  b going  to  infinity,  we  have  o2  = Ms2  and  p = 0,  while  for  b 

going  to  zero,  we  have  o2  = po2  = M2s2. 
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APPENDIX  B 


ANOTHER  FORM  FOR  THF.  PROBABILITY  DENSITY  FUNCTION 


The  case  we  are  considering  is  that  in  which  the  cross- 
correlation coefficient  r = 0 and  = a ^ (-*-R  = 0),  so  that  the  term 
X in  (2-21)  equals  zero.  Under  these  conditions,  the  sum  and  difference 
envelopes  Z-j  and  are  statistically  independent,  since  we  may  factor 
their  joint  probability  density  function  (2-14): 


P4tZl,Z2)  =SfeXP 


( ^ + I (S2,l}\  Z*  ( Z^+S^  T (S*Z2\ 

\ KWl  3 e p V ~WJ  0 R7 


21 


1 2 
- 


IL 


•m 

2Z0 

z2 

2 2 

X 

SJ;  2'2h4 

(B-l) 


where  a|  = lo2(l+n),  a?  = i-  a2 ( 1 -p ) , = S ~ ( h^+h^2hi cos  2e)/2(l+p), 

and  h2  = S|/2a2  = (h^+h|-2h^ h^  cos  20 )/2( 1 -p ) . x2(x;2,a2)  is  the  noncentral 

chi-square  distribution  with  two  degrees  of  freedom  and  noncentrality 
parameter  a2.  Thus  the  filter  output  y = -g  (Z|-Z|)  is  the  difference 
between  two  scaled,  independent  noncentral  chi-square  random  variables. 

This  type  of  relationship  was  noted  also  by  Lee  [16]  in  a similar  appli- 
cation. 

The  noncentral  chi-square  distribution  has  many  forms  and 
interpretations  (see  [17-20]).  For  example,  we  may  write 

(t  a2)r 


X2(x : 2.  a2)  = exp(-  i a2)  J 

c r=0 


r! 


x2(x;  2+2r) , 


(B-2) 


where  x2(x;  n)  is  the  (central)  chi-square  distribution  with  n degrees 
of  freedom: 

X2(x;  n)  = \ e'  ^x(x/2)’-n"1  [r(n/2)]'1 


(B-3) 
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B.l  Density  for  y using  characteristic  function 

Using  an  expression  given  by  Jayachandran  and  Barr  [21], 
we  may  write  the  characteristic  function  of  the  distribution  of  y as 

♦y<t>  * «P  C-(h|+h*)3  I i 

y J 4 m=0  n=0  n‘ 

m-l  5 ,n*n-k,  t?  k i m+I  f 


<Y>  & • > 


+ lbs.)  y (iwn-r, 

0-7^  ' 2 ! A>  ■ 1 


7i+io|t)r+1 


(B-4) 


Thus  we  obtain  for  the  density  of  y 
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, , 2 r , , : : ty  O-ehlf  [i  0+p)hJ] 

p6(y)  ■ £ exp  c-(h  +h )]  J-  a-a-  -2-„— 4 

m=0  n=0 


f(mfnn'k)(¥)'k  4ArPy  ■ 2 + *]■ » > » 


m 

i 

k=0 
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m k 

Defining  the  polynomials  gIVx)  = £ (m+""k)  fr  , we  obtain  finally 

m k=Q  n k. 


(B-5) 


p6(y)  = exp  [-(h|+h|)]  J J 


[i(l-p)h§f  [i(Up)^3 
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n=0  m=0 


m! 


n! 


g:  ijirfkrJ.  »< » 


(B-6) 
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B.2  The  polynomials  g[J(x) 

We  list  here  several  of  the  properties  of  these  polynomials. 


Definition:  tj(x)  * iC^)  fr  ; «£(x>  ■ em(x).  $-1 


(B-7) 


k*0 


m n 


Iterative  relationships:  gJJ(x)  = gr  + l G^fx)  * gJVj(x)  + lx)  (B-8j^ 

r=0 


also,  Gj(x) 


k=0 


<f)  CTk(x)  ■ .!.  <J>  - 1 (I)  £t<*>  <B-9> 


V m-k 
(m,n>p) 


k=0 


(n>m) 


k=0 


k'  m-k 


(m>n) 


Differentiation  formula:  G^(x)  = G^^x) 


(B-10) 


n m k 

Addition  formula:  G"(x+y)  = T f-r  G"  . (x) 

mi  . g K ! fli  l\ 


(B-ll ) 


Relationships  with  special  functions: 

The  polynomials  gJJ(x)  may  be  related  to  the  confluent  hyper- 
geometric functions  Y(a,  0;  x)  and  2F0(y»  6;  ; -1/x),  to  the  Laguerre 

polynomials  L^(x),  and  to  the  Whittaker  functions  W.  (x)  by  the 
1 A,U 


following  formulae. 

GJ|(x)  = jjj’l'i'-m,  -m-n;  x)  (B-12) 

m+n+1 

* — Y(n+1 , m+n+2;  x)  (B-13) 

vm 

* ^1-  2r0(-m,  n+1 ; ; -1/x)  (B-14) 

■ sr •*n  *a  Wx)  for  * ■ in  » ■ t1  (B-,5) 

r (-Dm  rm*n’1(x).  (B-16) 

m 


By  using  the  identity  (B-15),  our  expression  (B-6)  for  the  pdf  can  be 

shown  to  be  identical  to  that  in  [63,  p.  63,  eq.  28]. 
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B.3  Equivalence  of  forms 

We  now  show  that  the  form  we  have  just  derived  for  the 
probability  density  function  is  equivalent  to  that  of  (2-23)  under  the 
assumed  conditions 


• • ,m,  n so  « 

i.  i i i rrK) 


^m+n-kx  x a o 


m=0  n=0  m!n!  m ' m=0  n=0  k=0  ' n ' k!m!n! 


^k,m+k.n 
/m+n\  x a b 


r r r /np’rH  x a D 
m=0  n=0  k=0  " k!(m+k)!n! 


amtkxk 


L X fm+kjlk!  b) 


m=0  k=0 


CD  ^ 

cj0<7>  (2^>iFi(ra+1-1;  b) 

JL 

* •"  J0  <7>  V2^>  Lm<-b>- 


using  (Cummer's  transformation  on  the  last  step. 


(B-17) 
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APPENDIX  C 


INVERSION  OF  THE  DIRECTIVITY  FUNCTION 


The  directivity  function,  as  it  has  been  defined  in  Chapter  1,  Is 
D = U((5;  N)  = sinN(5/Nsin(5 

* jj-  [cos(5  + cos3d  + cos5(5  +...+  cos(N-l)0],  N even  (C-l) 

* J [1  + 2cos2d  + 2cos4d  +...+  2cos(N-l)^]t  N odd.  (C-2) 


In  this  appendix  we  derive  an  inversion  of  D,  that  is,  in  a way  to 
express  (5  = d(D;  N) . 

Recalling  that  cosnx  = 1 - 


= 1 - (nx)2/2,  (nx)2  « 12, 


we  may  write  for  D,  N even, 


D * | [1  - d/2  + 1 - 9(5 2/  2 + 1 - 25(52/  2 +...+  1 - (N-l)2(52/2] 

■ in  - 4 []2 + 2,2 + 52  +’--+ 

= 1 - (N2-l)(52/6,  [ (N-l  )0]2  « 12  (C-4) 

where  we  have  used  Gradshteyn  and  Ryzhik  [13].  formula  0.122.2.  For  N 
odd,  we  obtain 


where  once  again  we  have  used  [13],  formula  0.121.2.  Therefore,  for  N odd 
or  even  we  can  write 


02  « 6(1-D)/(N  -1 ) . 


( 06) 


If  we  use  the  first  three  terms  of  the  cosine  power  series 
representation,  we  obtain 

M - (N2-l)02/6  + H4  + 34  + 54  +...+  (N-l )4] 

* 1 - (N2-l)02/6  + ( 3N4- 1 0N2+7 )04/ 360 , [(N-1)0]2  « 30. 


(C-7) 


Application  to  Chapter  4. 

In  Chapter  4,  we  use  i - 0sina/2.  Therefore  we  have  the  approximation 
£2sin2a/4  = 6(1-D)/(N2-1) 


p2(l-cos2a)/2, 


or  a = arcsin 


in  2 1 

W-l 


arccos 


48(1-D) 

e2(N2-l) 


(C-8) 
2 


2 

, (^y-  psina)  «12 
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APPENDIX  D 
COMPUTER  PROGRAMS 

It  is  not  our  purpose  here  to  burden  the  reader  with  a full 
disclosure  of  the  many  numerical  procedures  we  have  employed.  We  shall 
list,  however,  for  the  benefit  of  those  who  can  interpret  them,  the 
major  computer  programs  we  have  written,  so  that  they  may  be  open  to 
scrutiny  and  available  to  those  who  may  wish  to  perform  similar  calcu- 
lations. 

The  programs  are  all  straightforwardly  written  in  the  BASIC 
language.  Those  programs  listed  are  labeled  as  follows: 

(1)  basic  multiplicative  probability  density  function  (2-26) 

(2)  basic  noncentral  chi-square  probability  density  function  (2-25) 

(3)  multiplicative  array  detection  analysis  (3-2,  3-10) 

(4)  standard  array  detection  analysis  (3-5,  3-13) 

(5)  multiplicative  array  estimator  mean  and  mean  square  (4-69,  4-70) 

(6)  standard  array  estimator  mean  and  mean  square  (4-53,  4-54) 

(7)  multiplicative  array  estimator  bound  (5-48) 

(8)  standard  array  estimator  bound  (5-38). 

Some  manual  calculations  were  employed  to  supplement  programs  5 and  7, 
using  [31]. 
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PROGRAM  ONE 


1 REM  BASIC  PROGRAM  F0R  P6(Y) 

5 DIM  G(100,100),H(100),J(100),U(100) 

10  LET  G(0.0)-1 

15  LET  H(0)-1 

20  LET  J(0)*1 

25  LET  U(0)«1 

98  REM  LINES  100  T0  335*  ENTER  RH0,  SIGMA  2,  P0WER  SNR,  AND  THETA | 

99  HEM  COMPUTE  CONSTANTS. 

100  READ  R,S,H,T 

200  PRINT  " PROBABILITY  DENSITY  FUNCTION" 

205  PRINT 

210  PRINT  "RH0«"R,  "SIG2«"S,  "PSNR*"H, "THETA “”T 
215  print 

220  PRINT  "Y","P6(Y)M 

305  LET  Cl-=2*H*(l-R*e0S(2*T))/(l-Rf2) 

310  LET  C2«H*(e0S(T))t2 

315  LET  C3=H*(SIN(T))t2 

320  LET  C4-(l-R)/(l+R) 

325  LET  C5— 2/(S*(l+R)) 

330  LET  C6=?./(S*(1-R)) 

335  LET  C?=4/(S*(l»Rt2)) 

400  READ  Al,A2tA3,L 
405  PpH  Y’=Al  T0  A2  STEP  A3 
410  LET  Y1=C7*ABS(Y) 

415  LET  P-O 

498  REM  LINES  500  T0  545*  CALCULATE  THE  F0LYN0MIALS  C(Y;K,N). 

500  -F0R  N*1  T0  L 

505  LET  G(0,N)=1 

510  LEI'  H(N)«C3*H(N-1)/N 

515  f0r  K-l  t0  L 

520  LET  U(K)»Y1*U(K-1)/K 

525  LET  G(Kf0)=G(K-l,0)  + U(K) 

530  LET  G(K,N)=G(K-lfN)  + G(K,N-1) 

540  LET  J(K)=C2*J(K-1)/K 
542  NEXT  K 
545  NEXT  N 

599  REM  LINES  600  TO  725?  COMPUTE  DENSITY  FUNCTI0N. 

600  IF  Y<0  THEN  7 00 

605  F0R  N=0  T0  L 

610  F0R  K=0  T0  I. 

615  LET  P=P  +(C4t(K-N))*j(K)*H(N)*G(K,N) 

620  NEXT  K 
625  NEXT  N 

630  LET  P=P*EXP(C5*Y-C1)/S 
635  G0  T0  800 
700  F0R  N*0  T0  I 

705  F0R  K=0  T0  L 

710  LET  P-P  4 (C4t(K-N))*j(K)*H(N)*G(K,K) 

715  NEXT  K 

720  NEXT  N 

725  LET  P=P*'EXP(C6*Y-Cl  )/S 
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800  PRINT  Y.P 

805  NEXT  Y 

1020  DATA  .2, 1,1,0 

1040  DATA  -7.5, 7. 5,. 5 

1041  DATA  30 

2000  END 

PROGRAM  TWO 

5 READ  A1,A2,A3 

10  READ  H,S 

15  PRINT  "YV'P(Y)" 

20  F0R  Y=A1  T0  A2  STEP  A3 

25  LET  U *'2*SQR (H*Y/S  ) 

30  G0SUB  300 

35  LET  p*=k*exp(-h-y/s)/s 

40  PRINT  Y,P 

45  NEXT  Y 

50  DATA  0,10,1 

55  DATA  1,1 

CO  STpfP 

300  LET  V=.5*U 

305  LET  K-l 

310  LET  G-l 

315  F0R  N-l  T0  20 

320  LET  G=G*V12/NI2 

325  LET  K=K  + G 

330  NEXT  N 

335  RETURN 

500  END 

PROGRAM  THREE 

5 dim  f(ioo),h(igo),e(ioo),u(ioo),v(ioo) 

10  LET  H(0)-1 

15  LET  U(0)=1 

20  LET  V(0)=1 

30  LET  F(0)=1 

32  READ  M,B,A,H1,S1,A4 

35  LET  E=.5*M*B*SIN(a) 

40  LET  D=SIN(E)/(M*SIN(E/M)) 

45  LET  H-Hl*(H*D)t2/(M  + 2*(M-i)*A4) 

50  LET  R=A4/(M  + 2*(M-l)*A4) 

55  LET  S=M*S1  + 2*(M-i )*A4*S1 

60  LET  T”E 

65  PRINT  " MULTIPLICATIVE  ARRAY  DETECTION  ANALYSIS" 

70  PRINT 

75  PRINT  "H-"M}"BETA="Bj"ALFA«"A;  "C0RREL»"A4 

SO  PRINT  "INPUT  SNR="H1 5 "Si *"S1 

5 PRINT  "RH0-"Rr,SIGt2="Sj"SNR«"H;MTHETA="T 
90  PRINT 

95  READ  L 

105  READ  A1,A2,A3 

110  PRINT  «y","PD","?F" 
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115  LET  Cl«H*(C0S(T))t2 

120  LET  C2 -H* (SIN (T ) )T2 

125  LET  83“(1-R)/(1+R) 

130  LET  C4*2/(l-R) 

135  F0R  Y«A1  T0  A2  STEP  A3 

140  LET  C5-2*Y/(S*(1+R)) 

145  LET  ?1«0 

150  F0R  N=0  T0  1 

155  LET  U (N+l ) =€2*U (N )/ (N+l ) 

160  F0R  M=0  T0  L 

165  LET  V(M+l)=€l*V(M)/(M+l) 

1?0  F0R  K«0  T0  M 

175  LET  F(m)=F(K)*(N+K+i)/(K+l) 

180  LET  E(K~K)-0 

185  F0R  J=0  T0  M-K 

190  LET  H(J+1)=C5*H(J)/(J+1) 

195  LET  E(M-K)=E(M-K)  + H(j) 

200  NEXT  J 

205  LET  P2*=U(N)*V(M)*(C31(M-N))  ' 

210  I£T  P3=F(K )* ( C4f(M-K ) )*E(M-K) 

215  LET  P1*=P1  + P2*P3 
220  NEXT  K 

225  NEXT  M 

230  NEXT  N 

235  LET  P4~-C5-2*H*(l-R*C05('2*T))/(l-Rt2) 

240  LET  P5=.5*(1+R)*EXP(P4)*P1 

245  LET  P6«.5*(l+R)*EXP(-C5) 

250  PRINT  Y.P5.P6 
255  NEXT  Y 

260  DATA  5 ,3. 141 59 ,.00001 

265  LATA  1,1,0 

270  DATA  30 

275  DATA  0,70,5 

280  END 

PROGRAM  FOUR 


5 DIM  U(200),H(200),E(200) 

10  LET  U(0)=1 

15  LET  H(0)=1 

20  READ  N,B1,A5,S1,H1,A4,L,K2,L1 

25  PRINT  " STANDARD  ARRAY  DETECTION  ANALYSIS” 

30  PRINT 

35  LET  T=.5*N*B1*SIN(A5) 

40  LET  D=SIN(T)/(N*SIN(T/N)) 

45  LET  S=N*S1  + 2*(N-l)*A4*Sl 

50  LET  H=Hl*(N*D)t2/(N  + 2*(N-i)*A4) 

55  PRINT  ’tN="Nj',BETA="Blj,,ALPHA="A5 

60  PRINT  "SI ="S 1 , "HI ="H1 , "A  ="A4 

65  PRINT  "S-"S, "H="H, "THETA  =,,,T 

66  PRINT  ML="L,"K2=''K2,"L1="L1 

- 67  PRINT 

68  PRINT  "YV'PD'V’PF" 

69  READ  A1,A2,A3 
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70 

F0R 

Y2-A1  T0  A2  FfEP  A3 

75 

LET 

Y1-Y2/S 

80 

LET 

A-SQR(2*H) 

85 

LET 

B-SQR(2*Y1) 

140 

IF  1 

(A+B)<20  THEN  200 

145 

LET 

K1 =(A4B-20)/ (K2+1 ) 

150 

LET 

A=A-K1 

155 

LET 

B=B«K1*X2 

200 

IF  A>B  THEN  220 

205 

LET 

X=.5*At2 

210 

LET 

Y=,5*Bt2 

215 

G0  T0  300 

220 

LET 

X-.5*Bt2 

225 

LET 

Y«.5*At2 

300  LET  Q1=0 

305  F0R  M=0  T0  L 

310  LET  U(M+l)=Y*U(M)/(M+l) 

315  LET  E(M)=0 

320  F0R  K=0  T0  M 

325  LET  H(K+1)=X*H(K)/(K+1) 

330  LET  E(H)-E(H)  + H(K) 

335  NEXT  K 

340  LET  Q1=Q1  + U(M)*E(M) 

345  NEXT  M 

350  LET  Q2=Q1*EXP(-X-Y) 

405  IF  A<B  THEN  416 

410  LET  Q=Q2 

415  G0  T0  425 

416  Gjfe’JB  550 

420  LET  Q-l  - Q2  + X*EXP(-X~Y) 

425  LET  P1=EXP(-Yl) 

430  PRINT  Y2.Q.P1 

431  NEXT  Y2 

480  DATA  4, 3. 14159,. 00001 

481  DATA  1,1,0 

498  DATA  120,. 93, 30 

500  DATA  0,50,5 

545  ST0P 

550  LET  V=.5*A*B 

555  LET  K=1 

560  LET  G=1 

565  F0R  N=1  T0  LI 

570  LET  G=C*W2/Nt2 

575  LET  K*=K  + G 

580  NEXT  N 

585  RETURN 

600  END 

PROGRAM  FIVE 

5 DIM  U(94,94),L(94,94),E(94),C(94) 

,-10  LhT  U(0,0)-1 

15  LET  L(0,0)=1 

20  LET  B(0)=4/15 


25  LET  B(l)— 32/315 

30  LET  C(0)«l/6 

35  LET  C(l)-l/20 

40  READ  M,B,H1,N2 
50  P0R  R«1  T0  N2 

55  let  B(R+1)-2*(R+1)*(2*R+3)*B(R)/((2*R+3,5)*(2*R44,5)) 

60  LET  C(R+i;-(R+l)*(2*R+3)*C(R)/((R+2)*(2*R+5)) 

65  NEXT  R 

100  LET  K=2*SQR(6/(4*Mt2-l) )/B 
105  PRINT  K 
110  LET  K1=M*H1 

115  F0R  A=.00001  T0  .21  STEP  .01 
120  LET  E-.5*M*B*SIN(A) 

125  LET  D=SIN(E)/(K*SIN(E/M)) 

130  LET  L=SQR(1+4*K1*B*2) 

135  LET  G=.25  + 2*Kl*Dt2*(SIN(E) )f2 
140  LET  G1=SQR(G) 

145  LET  A9=G*Lt2 

150  LET  B9=-2*K1 

155  LET  C9=*5*Lt2 

160  LET  L(1,0H-C9 

165  LET  L(2,0)=l-2*C9+.5*C9t2 

170  P0R  R*=0  T0  N2 

180  LET  U(R+1.0)=B9*U(R.0)/(R+1) 

190  LET  U(R,l)=A9*U(R,0)/(R+l) 

200  F0R  S-l  T0  N2 

210  LET  U(R+l,S)=39*U(R,S)/(R+S+i) 

220  LET  U(R,S+l)=A9*U(R,S)/{(R+S+l)*(S+l)) 

230  LET  L(S+1iR)“((2*S+R+1-C9)*L(S,B)-(R+3}*L(S-1,R))/(S+1) 

240  NEXT  S 

250  LET  L(0,R+l)=(CR+l)*L(0fR)-L(ltR))/C9 

260  LET  L(l,R+l )=((R+2)*L(l ,R)-2*L(2,R))/C9 

270  NEXT  R 

300  LET  HI  =0 

305  LET  M2=0 

310  F0R  R=0  T0  N2 

315  F 0R  S=0  T0  N2 

320  LET  HI -Ml  + U(R,S)*L(R,S)*B(R) 

325  LET  M2=M2  + U(R,S)*L(RfS)*C(R) 

330  NEXT  S 

335  NEXT  R 

340  LET  M1=M1*4*K*K1*EXP(-C9-2*G) 

345  LET  M2=M2*4*Kt2*Kl*EXP(-C9-2*G) 

350  PRINT  A, Ml ,M2,L,2*G1 
355  NEXT  A 

400  DATA  5,3.14159,4,93 

415  END 

PROGRAM  SIX 

5 DIM  U(100),B(100),L(100),C(100) 

.10  LET  U(0)*=1 

15  LET  B(0)«=4/15 

20  LET  B(l)=32/315 
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25  LET  L(0)«1 

30  LET  C(0)“l/6 

35  LET  C(l)-l/20 

40  READ  N,B,H1,N2 
45  LET  K-2*SQR(6/(Nt2-l))/B 
•50  LET  Kl-tf*Hl 

52  LET  U(l)— K1 

55  P0R  A-.00001  T0  .21  STEP  .01 
60  LET  E-.5*N*B*SIN(A) 

65  LET  D»6IN(E)/(N*SIN(E/N)) 

70  LET  H=N*H1*D12 

75  LET  L(l)=l-H 

80  LET  M1=B(0) 

85  LET  >12=C(0) 

90  F0R  Nl-1  T $ N2 

100  LET  U(Nl+l)=-Kl*U(Nl)/(Nl+l) 

102  LET  B(N1+1)=2*(N1+1)*(2*N1+3)*B(N1)/((2*N1+3.5)J'(2*N144.5)) 
105  LET  C(N1+1  )=(N1+1  )*(2*Nl+3 )*C(N1 )/((Ni+2)*(2*Nl+5)) 

110  LET  L(N1+1)=((2*N1+1-H)*L(N1)-N1*L(N1-1))/(N1+1) 

115  LET  Ml=^l  + U(Nl)*L(Nl)*B(Nl) 

120  LET  H2-M2  + U(N1)*L(N1)*C(N1) 

121  NEXT  N1 

122  LET  Ml =2*M1*K*K1*EXP ( -H  ) 

123  LET  M2=M2*2*K*2*K1*EXP(-H) 

125  LET  V1-=M2-M1*2 

130  LET  V2«=M2-2*A*M1+A1‘2 
135  PRINT  A,M1,M2,V1,V2 
140  NEXT  A 

145  DATA  10,3.14159,4 

150  DATA  99 

200  END 

PROGRAM  SEVEN 

5 DIM  U(94,94),L(94,94),B(94) 

10  LET  U(0,0)=1 
15  LET  L(0,0)=1 
20  LET  B(0)=4/l5 
25  LET  B(l)=32/315 
30  READ  M,B,H1,N2 
40  F0R  R«=l  T0  N2 

50  LET  B(R+l)=2*(R+l)*(2*2+3)*B(R)/((2*R+3.5)*(2*R44.5)) 

65  NEXT  R 

100  LET  K=2*SQR(6/(4*Mt2-l))/B 
105  PRINT  K 
110  LET  K1=M*H1 

115  F0R  A=. 00001  T0  .21  STEP  .01 
120  LET  E=.5*M*B*SIN(A) 

125  LET  D=SIN(E)/(M*SIN(E/M)) 

130  LET  L=5QR(l+4*Kl*Dt2) 

135  LET  G-.25  + 2*Kl*Dt2*(SIN(E))t2 

140  LET  Gl-SQR(G) 

' 145  LET  A9=G*Lt2 

150  LET  B9=-2*K1 

155  LET  C9-.5*Lt2 
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n 160  LET  L(l,0)«2-C9 

8 165  LET  L(2,0)-3-3*C9+.5*C9t2 

*'  200  F0R  R-0  T0  H2 

215  LET  U(R+1.0)-B9*U(R,0)/(R+l) 

220  LET  U(Rtl)-A9*U(R,0)/(R+l) 

225  F0R  S-l  T0  N2 

230  LET  U(R+1,S)=B9*U(R,S)/(R4S+1) 

235  LET  U(R,S+1  )=A9*U(R,S)/((R+S+1  )*(S+l)) 

ll  240  LET  L(S+lfR)=((2*S+R+2-'C9)*L(SlR)-(R+S+l)*L(S-lfR))/(S+l‘) 

245  NEXT  S 

246  LET  L(0,R+lW(R+2)*L(0,R)-L(l,R))/C9 

24?  LET  L(l,R+l)=((R+3)*L(lfR)-2*L(2,R))/C9 

250  NEXT  R 

300  LET  Ml-0 

310  F0R  R=0  T0  N2 

315  F0R  S*=0  TP  N2 

320  LET  M1=M1  + U(R,S)*L(R,S)*B(R)*(R+S+1-2*G)/(R+S+1) 

330  NEXT  S 

335  NEXT  R 

340  LET  Ml=Mlf2*Kt2*8*Klt2*L't2*EX?(-Lt2-4*G)/(Lt2-l ) 

350  PRINT  A,M1,L,2*G1 
355  NEXT  A 

400  DATA  5,3.14159,2,60 

415  END 

PROGRAM  EIGHT 

5 DIM  U(100),B(100),L(100) 

10  LET  U(0)-i 

15  LET  B(0)=4/15 

20  LET  Bfl}=32/315 

25  LET  L(o)=l 

40  READ  N,B,H1,N2 

45  LET  K=2*£QR(6/(N^2-l))/B 

50  LET  K1=N*H1 

52  LET  U(l)=~Kl 

55  F0R  A=. 00001  T0  .21  STEP  ,01 

60  LET  E=.5*N*B*SIN(A) 

65  LET  D=SIN(E)/(N*SIN(E/N)) 

70  LET  H=N*Hl*m2 

75  LET  L(l)*=2-H 

80  LET  M1=S(0) 

90  F0R  Nl=l  T0  N2 

95  LET  U(N1+1)=-K1*U(N1)/(N1+1) 

100  LET  B(Nl+l)=2*(Nl+l)*(2*Nl+3)*B(Nl)/((2*Nl+3.5)*(2*Nl+4.5)) 

110  LET  L(Nl+l)=(2-H/(Ni+l))*L(Nl)-L(Nl-l) 

115  LET  M1=M1  + U(Nl)*L(Nl)*B(Nl) 

121  NEXT  N1 

122  LET  Ml*=2*Ml*X*KTt2*D*EX?(-H) 

130  LET  Ml=Mlt2/(2*N*Hl*(2*H-l)) 

135  PRINT  A,M1 

'140  NEXT  A 
145  DATA  10, 3.141 59,. 5, 99 
200  END 
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ABSTRACT 

Analysis  is  made  of  a "multiplicative"  receiv- 
ing array  model  which  is  known  to  have  about  one- 
half  the  beamwidth  of  the  "additive"  or  convention- 
al (square-law)  array  with  the  same  number  of  ele- 
nents.  The  probability  density  functions  for  the 
filter  outputs  of  array  models  are  shown  for  mono- 
chromatic, planewave  signals  in  the  presence  of 
narrowband  Gaussian  noise,  for  broadside  arrays 
with  equal  spacing  and  uniform  gains.  The  number 
of  elements,  SNR,  source  bearing,  and  interelement 
noise  correlation  are  treated  as  parameters  of  the 
r<fistribution,  and  their  influences  are  displayed 
graphically. 

Based  on  these  probability  density  functions 
for  the  square-law  and  the  multiplicative  array 
processors,  we  then  examine  their  performances  as 
signal  detectors  and  as  bearing  estimators. 


INTRODUCTION 

In  judging  the  capabilities  of  sonar  array 
processing,  it  is  a corrrton  practice  to  use  the 
processor  output  SNR  (signal-to-noise  ratio)  as  a 
figure  of  merit  for  detection,  and  the  array  con- 
figuration's beamwidth  as  a measure  of  bearing  es- 
timation performance.  Thus,  based  on  these  cri- 
teria Arndtfl],  for  example,  concludes  that  con- 
ventional or  additive  processing  is  the  best  method 
for  detection  purposes,  while  correlating  or 
multiplicative  processing  is  the  choice  for  obtain- 
ing bearing.  These  conclus  ns  had  also  been 
reached  by  Nelsby  and  Tucker  [2]  with  some  reserva- 
tions. However,  we  wish  to  point  out  that  these 
criteria,  SNR  and  beamwidth.  while  often  adequate 
for  pursuing  the  optimization  of  a given  configura- 
tion (3],  are  not  always  suitable  for  motivating  a 
choice  between  two  different  configurations. 

The  usual  standard  of  detection  performance  is 
the  probability  of  detection  (PD)  achieved  for  a 
given  false  alarm  rate  (PFA)  at  an  assumed  system 
input  SNR  (Neyman-Pearson  model  of  detection). 
Therefore,  even  though  it  is  reasonable  to  assume 
that  the  PO  is  proportional  to  output  SNR  for  each 
of  two  different  systems,  in  general  one  cannot  say 
that  the  PD's  are  comparable  if  the  SNR  is  the  same 
since  there  is  yet  another  parameter  to  be  speci- 


fied: the  detection  threshold,  or  equivalently  the 
the  false  alarm  rate  on  which  the  threshold  is 
based. 

Similarly,  while  beamwith  (resolution; is  rer- 
tainly  a factor,  bearing  estimation  performance  is 
to  be  understood  in  terms  of  minimum  error  attain- 
able from  an  arrangement  in  which  the  bearing  is 
extracted,  or  indirectly  measured,  from  observa- 
tions corrupted  by  noise.  Therefore,  a true  com- 
parison of  two  configurations  with  respect  to  bear- 
ing estimation  requires  an  accounting  of  the  effi- 
ciency with  which  each  smooths  the  noise. 

The  limitations  we  have  pointed  out  are,  of 
course,  more  or  less  understood  by  those  employing 
the  conventional  figures  of  merit.  The  limitations 
are  tolerated  because  of  the  difficulty  in  obtain- 
ing the  precise  quantities  on  which  systems  com- 
parisons need  to  be  based,  quantities  which  are  re- 
lated to  the  statistical  behavior  of  the  system 
outputs.  In  this  paper  we  present  the  probability 
density  function  of  the  detector  output  for  a class 
of  multiplicative  arrays,  and  employ  it  to  calcu- 
late detection  and  bearing  estimation  capabilities 
of  these  arrays  in  comparison  with  conventional  or 
additive  arrays  which  use  a square- law  detector. 
Though  necessarily  we  select  a narrow  class  of  sys- 
tems in  order  to  achieve  the  rigor  we  seek,  both 
the  analytical  approach  and  the  results  (described 
mere  fully  in  [4 J ) lend  themselves  to  extensions, 
thereby  encouraging  more  thorough  system  analysis 
in  the  field  of  sonar  arrays. 

MULTIPLICATIVE  ARRAY  PROBABILITY  DENSITY  FUNCTION 

Consider  the  model  of  a multiplicative  array 
as  shown  in  Figure  1,  in  which  the  2M  omnidirec- 
tional receiving  elements  are  equally  spaced  d 
units  apart  alc.ig  a straight  line.  The  sum  of  one 
half  of  the  array  element  outputs  is  multiplied  by 
that  of  the  other  half,  and  the  product  is  lowpass 
filtered.  We  make  the  following  assumptions: 

SIGNAL:  the  array  intersects  a planewave, 
monochromatic  signal  with  known  amplitude  A and  un- 
known bearing  a,  relative  to  the  perpendicular  of 
the  array. 

This  work  was  supported  in  part  by  the  Office  of 
Naval  Research  under  Contract  N00014-67-A-0377. 
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FIGURE  1 . Multiplicative  array  model. 


NOISE:  the  array  elements  experience  station- 
ary, zero-mean  Gaussian  noise  such  that  the  vari- 
ance at  each  element  is  s2  and,  at  a given  instant 
of  time,  adjacent  elements  only  are  correlated, 
with  coefficient  a. 

Under  these  conditions  the  sum  channels  of 
Figure  1 have  the  signal  terms 

Si(t)  * HADMcos(u)t-0.),  i - 1,2 

where  the  directivity  function  is  given  by  DM  = 
sinHM0sina)/Msin(i6sina),  with  6 = 2-nd/X  and 
20  ■ 01-02  = MBsina.  The  noise  terms  n. (t)  have 
equal  variances  o2  = Ms2  + 2(M-l)as2  and  are  cor- 
related with  coefficient  p = a/[M  + 2(M-l)a]. 

The  filter  output  is  given  by  y ( t)=J-(Z 1 2-Za 2) 
where  Zi  is  the  envelope  of-f(ui  + u2)  and  Z2  is 
the  envelope  of4(ui  - u2),  assuming  ideal  lowpass 
filtering.  That  is,  y is  the  difference  between 
two  independent,  scaled  noncertral  chi-square  ran- 
dom variables  (for  n^  of  unequal  variance,  the  Z. 


art  dependent).  2t  can  be  shown  [4]  that  the 
probability  density  function  of  y in  this  case  is 


pn(y)  ■ pa(y;  a,  e,  i,  s,  h,  a) 
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where  B = o2/Ms2,  C * pB,  and  H2  * A2/2s2,  and 
where  the  I ( ) are  the  modified  Bessel  functions 
of  the  first  kind  of  order  r,  and  the  L ( ) are 
the  r:th  degree  Laguerre  polynomials.  r The  un- 
correlated noise  case  corresponds  to  B = 1 and 
C = 0. 


Results  for  the  multiplicative  array  probabil- 
ity density  function  (1)  are  presented  in  Figures 
2 through  5.  The  nominal  case  of  2M  * 10  elements, 
uncorrelated  noise  (a=0),  unit  noise  variance  (s2 
*1),  zero  bearing  (a=0),  half-wavelength  array 
spacing  (B=it),  and  zero  dB  input  SNR  (H2=l)  was 
chosen.  In  Figure  2 the  SNR  is  varied,  while  the 


FIGURE  2.  Multiplicative  array  PDF  ( equation  1) 
for  input  signal- to-noise  ratio  varied. 
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FIGURE  3.  Multiplicative  array  PDF  (equation  1} 
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FIGURE  4.  Multiplicative  array  PDF  (equation  1) 
for  noise  correlation  coefficient  varied. 

number  of  elements  is  varied  in  Figure  3.  It  is 
seen  that  these  two  parameters  have  similar  ef- 
fects on  the  distribution.  The  effect  of  adjacent 
Interelement  correlation  is  displayed  in  Figure  4, 

FIGURE  5.  Multiplicative  array  PDF  (equation  1) 
for  l>earing  angle  varied. 
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FIGURE  6.  Square-law  array  model. 

showing  how  increased  correlation  increases  the 
variance  of  the  filter  output,  and  decreases  the 
value  of  y at  which  the  peak  cf  the  distribution 
occurs.  Figure  5 demonstrates  a phenomenon  of 
interest,  for  the  bearing  is  seen  to  act  as  a lo- 
cation parameter  for  the  distribution,  indicating 
the  potential  of  using  the  filter  output  for 
bearing  estimation. 

SQUARE-LAW  ARRAY  PROBABILITY  DENSITY  FUNCTION 

Consider  the  conventional  or  square-law  array 
with  N elements,  as  modelled  in  Figure  6.  Using 
the  same  basic  assumptions  as  above,  it  can  be 
shown  that  the  filter  output  y is  a scaled  noncen- 
tral chi-square  random  variable  with  probability 
density  function 

Ps(y)  = Ps(r.  O,  B.N,  s,  H,  a) 

FIGURE  7.  Square-law  array  PDF  (equation  7) 
for  input  signal-to-noise  ratio  varied. 
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where  B*  * [N  + 2(N-l)a]/N 

and  * sin(}Nesina)/Nsin(iBsina). 

In  Figures  7 and  8,  the  probability  density  func- 
tion of  the  filter  output  of  the  square-law  array 
model  is  seen  to  behave  with  parameter  variation 
In  much  the  same  way  as  that  of  the  multiplicative 
Here  we  have  used  H2  = 1 , s2  = 1 , 8 = it,  a = 0,a=0, 
and  N a 10  as  the  nominal  case.  In  Figure  7 the 
SNR  is  varied,  and  in  Figure  8 the  source  bearing 
Is  varied. 


DETECTION  PERFORMANCE 

Having  now  the  probability  density  functions 
for  the  array  processor  outputs,  we  are  in  a posi- 
tion to  evaluate  what  Tuck<_r  [5]  described  as 
"...the  difficult  question  of  whether  the  better 
signal/noise  performance. . .is  really  indicative  of 
a higher  probability  of  detection."  That  is,  we 
are  able  directly  to  calculate  the  probabilities 
of  detection  for  the  two  array  models,  as  functions 
of  SNR  and  of  corresponding  probabilities  of 
false  alarm. 


FIGURE  8.  Square-law  array  PDF  (equation  2) 
for  the  bearing  angle  varied. 


maximum  PFA)  to  be  used  to  pronounce  one  the  better 
detector  with  some  confidence  in  the  generality  of 
this  kind  of  statement.  Helmstrom  [6]  cites  addi- 
tional uses  of  the  ROC. 


Assuming  the  decision  "signal  present"  is 
made  if  the  filter  output  y exceeds  a threshold  t, 
the  probability  of  false  alarm  or  false  alarm  rate 
and  the  probability  of  detection  are  given  by 


PFA  p(y|h=0)dy,  PD  =/  p(y|H>0)dy.  (3) 

For  the  multiplicative  array  model  we  have  (t>0) 
PFA  ■ «1  . C/«)  explFf:^!  («) 


PD  - exp 


-2x 

Ms*(B+C) 


2MH2n;B-Ccos29l 


jLILJY  . 
B UM 


b^T7) 


For  the  square-law  array  model  we  have 
PFA  = exp(-T/3'Ns2) 
and  PD  = Q(HDN/2?7BT,  /«!t/B'Hs2), 
where  Q is  Marcum's  Q-function. 


(5) 


(«) 

17} 


ROC  curves  for  our  two  models,  for  the  uncor- 
related, boresight  (a=0)  case  are  given  in  Figure 
9 for  N=2M=4  and  in  Figure  10  for  N-C'HO.  The 
fact  that  the  square-law  array  processor  curves 
are  consistently  above  those  of  the  multiplicative 
indicates  that  the  square-law  array  model  is  a 
better  detector  in  this  case.  They  are  quite 
close  together,  however,  and  whereas  Arndt  [l] 
gives  an  SNR  difference  of  3 dB,  cur  results  imply 
that  the  minimum  detectable  signals  do  not  differ 
by  more  than  one  dB. 

MAXIMUM  LIKELIHOOD  BEARING  ESTIMATES 


If  the  filter  output  y of  the  multiplicative 
array  model  is  to  be  operated  upon  to  obtain  an 
estimate  of  the  bearing  a,  it  is  shown  in  [4]  for 
uncorrelated  noise  the  maximum  likelihood  bearing 
estimator  has  the  approximate  forms 


a(y)  - iarccos[ 


1 - 


4B 

(F (M2-l ) 


low 


SNR 

(8) 


°<y)  ' 

Also,  the  corresponding  forms  for  the  square-law 
model  are 


It  can  be  shown  that  the  multiplicative  model 
has  a more  desirable  false  alarm  rate  for  she  same 
value  of  detection  threshold,  while  the  sauare-law 
model  yields  a higher  PD.  However,  this  apparent 
tradeoff  is  blurred  by  the  unsuitability  of  the 
detection  threshold  value  as  a basis  for  ccmoari- 
son.  The  more  informative  "receiver  operating 
characteristics"  (ROC),  given  by  PD  - f (PFA;  SNR), 
render  two  different  threshold  detectors  more  di- 
rectly comparable,  enabling  a performance  index 
such  as  "minimum  detectable  signal"  (the  SNR  re- 
quired to  insure  a given  PD  for  a specified 


a(y)  = ‘arccos|1  ' FfF^ryl1  ' 752iVy 1 n ( "Tis7^)] I 

low  SNR  (10) 

«(y)  1 arcsinly^fl  - ^j_),  high  SNR.  (11) 

In  [4]  the  mean  vilues  of  these  estimators  for 
high  SNR  and  for  N = 2M  are  calculated  to  be 

E(«)  » 2kKe'h2^  ^"W^m^L^h2)  (12) 

square  law 
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and  E(a)  » K - KQ(X.2^)  ♦ 2kKexpHA2-2g) 

n 0 multiplicative 

with  K * 4.9/BrfFT,  k “ NH2,  h2  * NH20j*/2. 

X*  « 1 + 2h2,  4g  * X2  - 2h2cos2e,  and  the  B{  , ) 

•re  beta  functions. 

Computations  of  these  mean  values  are  given  in 
Figures  11  and  12,  with  a straight  line  reoresent- 
Ing  an  unbiased  estimator  drawn  in  for  reference. 
The  trend,  as  the  SNR  increases,  toward  unbiased- 
ness--where  E(a)  = a,— is  clearly  seen  for  the 
smaller  values  of  bearing.  On  the  whole,  the  two 
estimators  appear  to  behave  quite  similarly  with 
respect  to  bias. 

BOUNDS  ON  ESTIMATOR  ERROR 

As  a means  of  learning  estimator  performance 
we  have  chosen  the  Cramer-Rao  bounds.  Using  the 
above  expressions  for  the  mean  values  of  the  bear- 
ing estimates  and  working  with  the  filter  cutout 
probability  density  functions,  we  may  calculate  the 
familiar  Cramer-Rao  bounds  on  estimator  variance: 


Var(a)  > [^(a)]2/E{[^|ki5l]*}. 


An  example  of  this  bound  for  each  array  model  Is 
given  In  Figure  13,  in  which  for  small  bearings  the 
square-law  array  is  seen  to  have  a smaller -upptr 
bound  on  variance.  It  should  be  kept  in  mind  that, 
since  there  is  a bias,  the  mean  square  error  will 
be  greater  than  the  variance.  These  curves,  as 
well  as  others  we  have  calculated  for  actual  vari- 
ance and  mean  square  error,  show  a smaller  error 
for  the  square-law  array  for  small  (major  lobe) 
bearings,  while  the  fact  that  the  curves  cross  at 
higher  values  of  bearing  seems  to  imply  that  there 
are  values  at  which  the  multiplicative  array's 
estimator  is  superior. 

CONCLUSION 

He  have  very  briefly  summarized  recent  inves- 
tigations which,  by  direct  calculation,  demonstrate 
that  the  actual  relative  performances  of  square-law 
and  multiplicative  array  processors  in  signal  de- 
tection and  in  bearing  estimation  can  differ  sig- 
nificantly from  those  predicted  using  the  systems' 
SNR  and  beamwidth  as  figures  of  merit. 


FIGURE  9.  Receiver  zpetitlng  cV^racreri's tics C,V=4i 
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FIGURE  11.  Mean  value  of  square-lav  array's 
bearing  estimate  vs.  actual  bearing. 

REFERENCES 

1.  L.  K.  Arndt,  "Responses  of  arrays  of  Isotrop- 
ic elements  in  detection  and  tracking,"  U.  S.  Navy 
Electronics  Laboratory  Research  Report  1456.  21 
April  1967. 

2.  V.  G.  Welsby  and  D.G.  Tucker,  "Multiplicative 
receiving  arrays , "journal  of  the  British  ire,  19, 
pp.  369-382  (June  1959). 

3.  R.  J.  Whelchel,  "Optimization  of  conventional 
passive  sonar  detection  systems,"  Naval  Avionics 
Facility  (Indianapolis)  Technical  Report  1440, 

26  June  1959. 

4.  L.  E.  Miller,  "Signal  detection  and  bearing 
estimation  capabilities  of  multiplicative  array 
processors,"  Ph.D.  dissertation,  The  Catholic 
University  of  America,  Washington,  D.  C.,  1973. 

5.  D.  G.  Tucker,  "Signal/noise  performance  of 
multiplier  (or  correlation)  and  addition  (or  inte- 
grating) types  of  detector,"  (British)  IRE  Mono- 
graph Number  120R,  February  1955. 

6.  C.W.  Hel Strom,  Statistical  Theory  of  Signal 
Detection,  Pergamon  Press,  New  York,  1960. 


FIGURE  12.  noan  value  of  multiplicative  array's 
bearing  estimate  vs.  actual  bearing. 
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tion error  vs.  actual  bearing. 
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of  an  Analog  Cross-Correlator  with  Correlated 
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Abstract — The  probability  density  function  (pdf)  for  the  output  of  an 
analog  cross-correlator  with  correlated  bandpass  inputs  is  derived.  The 
pdf  is  derived  by  a “direct  method”  without  resorting  to  the  “charac- 
teristic function  method,”  which  usually  requires  contour  integrations  in 
a complex  plane  for  inversion  operations.  The  correlator  consists  of 
bandpass  filters,  a multiplier,  and  a zonal  low-pass  filter.  We  treat  the 
general  situation  in  which  the  two  inputs  are  narrow-band  signals  of 
unequal  power  and  of  different  phases.  The  bandpass  input  noises  are 
assumed  to  be  correlated  and  may  have  different  powers.  In  the  Appendix, 
another  derivation  for  the  pdf  is  given  in  the  special  case  of  equal  power 
correlated  noise.  This  derivation  is  based  on  the  fact  that  the  correlator 
output  random  variable  is  the  difference  of  two  independent  noncentral 
chi-square  variables  of  two  degrees  of  freedom.  We  show  that  the  two 
expressions  for  the  pdf  (one  from  the  direct  method  and  the  other  from 
the  characteristic  function  method)  are  indeed  equivalent.  Finally,  we 
discuss  two  major  areas  of  application. 

I.  Introduction 

THE  PROBABILITY  density  function  (pdf)  for  the 
output  of  an  analog  cross-correlator  with  correlated 
bandpass  inputs  is  derived.  The  correlator  consists  of 
bandpass  filters,  a multiplier,  and  a zonal  low-pass  filter. 
We  treat  the  general  situation  in  which  the  two  inputs  are 
narrow-band  signals  ot  unequal  power  and  of  different 
phases.  The  noises  in  the  input  channels  are  assumed  to  be 
correlated  with  unequal  power. 

The  problem  of  obtaining  the  pdf  for  the  output  of  a 
cross-correlator  whh  bandpass  inputs  has  been  a matter  of 
continuing  interes*  [1]— [7],  To  reduce  the  complexity 
inherent  in  the  mathematical  aspects  of  the  problem, 
previous  authors  have  often  restricted  their  considerations 
to  certain  simplifying  assumptions  or  approximations  in 
deriving  a general  expression  for  the  pdf.  However,  Andrews 
[5]  is  an  exception;  he  has  obtained  a general  expression 
of  considerable  complexity  by  utilizing  the  characteristic 
function  method.  In  this  paper  we  show  how  the  pdf  can 
be  derived  in  a direct  fashion  without  employing  the 
characteristic  function  method.  The  resultant  expression 
for  the  pdf  appears  to  be  much  simpler  than  that  of  Andrews. 
In  the  Appendix,  yet  another  derivation  of  the  pdf  for  the 
correlator  output  is  given.  This  derivation  is  based  on  the 
fact  that  the  correlator  output  random  variable,  manip- 
ulated suitably,  a difference  of  two  independent  non- 
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central  chi-square  variables,  each  with  two  degrees  of 
freedom,  Under  this  formulation  for  the  output  variable, 
the  application  of  the  characteristic  function  method  gave 
an  easy  derivation  for  the  pdf.  These  findings  resulted  from 
a study  of  the  detection  capabilities  of  narrowband  multi- 
plicative array  configurations  [8]. 

The  “direct  approach”  we  have  used  to  obtain  the  filter 
output  probability  density  function  is  by  no  means  the  only 
one.  Perhaps  the  most  thorough  method  is  that  which 
takes  into  account  the  filter  transfer  function  of  a realizable 
filter  and  treats  the  filter  output  as  a time  series.  In  the 
literature,  this  approach  is  often  attributed  to  Kae  and 
Siegert  [9],  and  usually  is  found  applied  to  the  square-law 
detector  [10]-[I2],  An  exception  is  the  work  of  Lampard 
[3],  who  studies  a multiplier/filter  combination. 

Other  references  which  use  the  same  “zonal  filter”  model 
as  we  do  to  analyze  the  multiplier/filter,  though  less  gener- 
ally, are  [1],  [2],  [4],  and  [5],  A thorough  treatment  of  the 
output  SNR  of  such  systems  is  that  of  Green  [13],  In  this 
area,  very  similar  pdf  curves  were  obtained  from  sometimes 
very  dissimilar  expressions. 

II.  Tm  Cross-Correlator  and  Formulation  of  the 
Problem 

The  cross-correlator  to  be  studied  is  shown  in  Fig.  1. 
The  inputs  to  the  multiplier  an:  two  narrow-band  processes 
// , ( z ) and  «,(/),  and  their  product  :(t)  is  passed  through  a 
zonal  low-pass  filter  to  yield  the  output  v(t),  as  shown.  The 
processes  i/ , ( z ) and  i/2( z ) are  said  to  con.u-i  of  the  super- 
position of  the  deterministic  signals  .v, (z ).  .v,(z)  and  the 
noise  processes  «,(/),  «,(/),  respectively.  At  the  same  in- 
stant of  time.  /i,iZ)  and  «,(/)  are  assumed  to  be  jointly 
Gaussian  random  variables  with  variances  it, 2 and  a,2, 
correlation  coefficient  />,  and  zero  means. 

The  multiplier  is  assumed  to  be  instantaneous,  so  that  its 
output  can  be  written 

:{t)  = w,(z)  x u2U) 

= f V , ( z ) 4-  Z7, (z )]  x [r,(f)  + Z7;(z )] 

= J[m,(z)  + ;/;(z)J2  - ![»/,( z ) - u2{t )]2 


= [u<z)  + n ,(z)]2  - [v4lz)  + zzjfz)]2 

0) 

where 

■'.1,4(0  = J['|(0  ± v2(z)l 

(la) 

and 

Z7  x 4<  Z ) = }[«,(/)  ± «>(Z  )]. 

(lb) 

This  arrangement  may  be  recognized  as  the  old  “quarter- 
square  multiplier"  idea  used  in  analog  computation. 
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Fig.  1.  Block  diagram  of  cross-correlator. 


For  the  new  variables  we  have  defined,  we  have  the 
following  moments: 


What  is  meant  by  “zon:,l  low-pass  filter”  is  that  the  filter 
output  y is  equal  to  the  multiplier  output  r,  less  the  terms 
of  frequency  2c o.  Thus  if  we  define  the  sum  and  difference 
terms 

z,(/)  A sft)  + n3(t)  = i[M,(t)  + u2(0] 

= Z,(/)  cos  [co/  - 0,(/)] 

= zlc(t)  cos  wt  + z,,(/)  sin  art 

Zj(/)  A Si(t)  + n4(/)  = - «2(0] 

= Z2(t)  cos  [art  - 02(/)] 

- z2c(r)  cos  wt  + z2s(t ) sin  wt 


(4a) 


we  find  that 


£(«3)  = 0 

£(m4)  = 0 (2a) 

£(«j2)  A a32  = i(<7,2  + lpaxa2  -)-  a22) 

i-W)  4 - 2pcr1cr2  + a22) 

E(n3nA)  A £<t3<t4  = |(<7,2  - a2).  (2b) 

Note  that  if  it,  = a2,  then  R = 0;  that  is,  «3  and  n4  are 
uncorrelated. 

Specialization  to  what  we  are  calling  the  narrow-band  case 
consists  in  the  conventional  assumption  that  the  following 
quadrature  representations  apply: 

•*3(0  = S3(t)  cos  [ft)/  - 03(/)] 

MO  = MO  cos  [w,  - 04(/)] 

”3(0  = «3<r(0  cos  + «3s(/)  sin  wt 

= N3(/)  cos  [ft)/  - 03(/)J 

n4(/)  = n4f(/)  cos  co/  + n4s(/)  sin  wt 

= jV4(/)  cos  [ a)/  - 04(/)]  (3) 

where  if 


>(/)  = i[z,2(0  - z22(/ )] 


(4b) 

(5) 


where 


Z1>22  = i[5,2  + £22  + N2  + N22  ± 2SyS2  cos  (02  - 0,) 
+ 2S,  Afj  cos  - 0,)  ± 25,^  cos  (02  - 0,) 
± 2S2/V,  cos  (0,  — 02)  + 2 S2N2  cos  (02  — 02) 
± 2Ar,N2  cos  (0,  - 02)] 

= ^3.42  + NiA2  + 253  4jV3  4 cos  (03  4 - 03  4) 

(5a) 


and 


ian  0i2  — 


S,  sin  0,  + S2  sin  02  + N,  sin  0,  + JV2  sin  02 


S,  cos  0,  ± S2  cos  0,  + V,  cos  0,  ± A’2  cos  02 

(5b) 


= -53,4  sin  03  4 + A'3,4  -in  0,  4 
53.4  cos  034  + N3  4 cos  03  4 


Our  object  is  to  find  the  pdf  of  the  filter  output  y(/)  at  a 
given  instant  in  time.  (From  this  point  on,  reference  to 
time  will  be  suppressed.) 


s,(r)  - S i(t)  cos  [cu/  - 0{( / )] 

»i(0  = N,(t)  cos  [wt  - 0 ,(/ ) ] 
where  i = 1,2  then 

S3,42  = i[S,2  + S22  ± 2S,52  cos  (0,  - 0?)] 
S,  sin  0,  ± S2  sin  02 


(3a) 


Distribution  of  the  Sum  and  Difference  Terms 

If,  to  use  vector  notation,  we  refer  to  the  input  noise 
components  as  .r'  = (Hlf,H1J,«2f,«2j),  sphere  the  prime  (') 
is  used  to  indicate  the  transpose,  then  we  may  express  the 
joint  pdf  of  the  input  noise  components  as 


Po(-0  = 


47i2\  det  K 


= exp  [—  jx'Kx  '*] 


tan  03  4 = 


S,  cos  0,  ± S2  cos  02 
^3.42  = i[N,2  + N22  ± 21V1iV2cos(0l  - 02)] 
N | sir  0,  ± N,  sin  02 


tan  03i 


/V,  cos  0|  ± JV2  cos  02 


(3b) 


where  A^,  the  covariance  matrix,  is  postulated  to  be 
Kx  = cov  [//,<■(/), n,.,(/).«2c(/),«2i(/)] 

= £[.r'.r] 


With  this  narrow-band  representation,  wc  have  for  the 
output  of  the  multiplier 


^ ('‘3  + «3>2  - (V4  + «4)2 


pa  ,<r2  — r<r  t <t2 


\a  2 
TO  j IT  2 
2 


a-, 


ro  ,<r2 

/)(T,(T2 

0 


l_  fG \G 2 pG j fT 2 


0 


(7, 


= ([S,2  + 2 .S', A’,  cos  (03  - 0j)  4-  A',2  - .S42 
- 2S4A'4  cos  (04  - 04)  - Af42] 

4-  (terms  with  frequency  2w).  (4) 

15 


That  is,  the  noise  correlation  at  the  inputs  to  the  multiplier 
is  considered  to  be  such  that,  for  the  same  time  instant 


T(nlcn2c)  = £(«,, h2s)  = paia2 
£(«u«2S)  = ~E(nlsn2c)  = ra  ,<r  2. 


..  ..... 
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We  may  consider  the  sum  and  difference  noise  com- 
ponents n = (n3c,n3s,nic,nis)  as  a linear  transformation 
n = Ax,  where  the  matrix  A is  given  by 


The  quadratic  form  appearing  in  the  integral’s  exponent 
reduces  to 


A = 


i 0 4 0 

0 } 0 } 

4 0 -4  0 

0 i 0 -ij 


Q = J5  [<T4'(Z'2  + Sj2)  + a'2(Zl 1 + s"2) 


We  may  then  write  the  pdf  of  the  variables  n as 


1 


p'(n)  = exp 
4irD 


(6) 


- 2<t3<t4A'S3S4  cos  (04  - 03  - x) 

- 2<t3<t4A'Z1Z2  cos  ($,  -</>,-  x) 

- 2<r4 VZ , cos  (<j>,  - t)  - 2<t31VZ2  cos  (<p2  ~ «)] 

(12) 


where 


where  the  new  covariance  matrix  K„  is  easily  computed  to  be 


K.  = AKrA 


0 


0 


/?<t3<t4  /?'cr3'T4 

- /?'<t3<t4  Ra 3<t4 


Ro3oa  — /?'<x3<r4 

L - R'OiO*  «<r3rr4 


0 


0 


(7) 


D = x det  A;  = <t32<t42[1  - R2  - (A')2].  (7a) 


X2  = R2  + (R')2  tan  v = R'/R 
V2  = ct^Sj2  + ai2X2Si2 

— 2<t3<t4A'S3S4  cos  ( 03  — 0A  + x) 
w2  = <t32s42  + <t42a:2s32 

- 2<t3<t4.YS3S4  cos  (03  - 04  + x) 
a 4S3  sin  fl3  — <t3A'S4  sin  (04  — x) 


tan  f = 


The  variances  a32  and  <r42,  and  R are  as  given  in  (2),  and 
A'<t3<t4  = —\raxa  2. 

It  follows  then  that  the  quadrature  expansions  of  the 
signal  plus  noise  variables  r,  and  r2  as  just  defined  have  the 
density  function 


tan  w = 


<r4S3  cos  03  - <t3AS4  cos  (04  - x) 
<x3S4  sin  04  — <t4A'S3  sin  (03  + x) 


<x,S4  cos  ()4  - <t4A'S3  cos  (0}  + x) 
Using  the  relationship 


• (13) 


1 


P2(=ir.=ir=2,.=j»)  = t-t- e*P  [ ~ i(::  - s)'K„-'{z  - S)] 
4 z'D 


exco''w  = V £„/„(x)  cos  mv,  c„  = L’ 


n = 0 


n = 0 


n > 0 


(8) 


where 

U - s)’  - (:lc  - S3  cos  0„  :is  - S3  sin  </-, 

:2c  - S4  cos  04,  r2j  - sin  04).  (9) 


and  performing  the  integration  with  respect  to  0,  and  </>2, 
ii:  a manner  quite  similar  to  Middleton  [14,  ch.  9]  we 
obtain 


III.  General  Expression  of  the  Probability  Density 
Func  tion  for  tiii  Correlator  Output 


Recall  that  the  filter  output  is  y = 1(7, 
density  function  of  7,  and  Z2  is 


Z2:).  The 


PAZ 


,.Z;)  = | c/0 1 f c/02p3(Z1.0,.Z2.02) 

J (>  J 0 

= Z,Z2  </0,  </r/>2p2(Z,  cos0,.  Z,  sin  0 , . 

* i)  /a 

Z2  cos  02.  Z2  sin  02 ) 

/ 7,  f2*  l’2" 

= CJp  (I0  c/0 , 

4-0  J0  J0 

•exp  [ - i(Z  - j)1  A'„  ‘<Z  - j)] 


p4(Z,.Z,)  = 'l^exp  |~  [cr4:Z,2  + ct32Z22  + U2]j 
V / Z|\  r /flr3lf  Z2\ 

■ lm  \ Cos  m(r  - u’  — x)  (14) 

\ n * 


\ D 

where  wc  have  used 


c 2 = cr42532  + (T.,2.S4:  - 2fT3a4.\  .V,.S4  cos  (fl4  - 03  - v), 

(15) 


(10a) 


and  /„,(  •)  is  the  111th  ordcr  modified  Bessel  function  of  the 
first  kind. 

To  obtain  the  density  for  i\  wc  use  the  transformations 


(10bt 


where 
(Z  - sy 

- (7,  cos  0 , - .S'3  cos  03,  /,  sin  0,  - .S',  sin  03, 

/2  cos  </>,  - ,V4  cos  d4.  Z2  sm  02  --  .S4  sin  04). 


/,  — \ 2r  cosh  u. 

/.  \ 2r  smli  1/.  for  /,  > /2,  or  v > 0 and  u > 0 

and 

7,  = \ — 2r  sinh  u, 

Z2  — \ —2)  cosh  1/,  for  7,  < 7,,  or  r < 0 and  u > 0. 


(11) 

16 


(16) 


- - - 
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Then  we  find  that  the  density  for  y is  given  by  the  expression 

Pf,(y)  = [ dups(y,u) 

Jo 

f X / f 

I du  p4(\2y  cosh  u,  V2 y sinh  u),  y > 0 

I du Pi(\' - 2 y sinh  u,  V —2y  cosh  u),  y < 0 
'Jo 

= — exp  ~ -E*] 

D L 2D  J 


X /*  X 

J]  £m  cos  m(p  - w - x) 
m = 0 J $ 


du  sinh  2ue~‘,<:osh  2u 
• Im(b  cosh  u)Im(c  sinh  u)Im(d  sinh  2 u)  (17) 


with 


a = |>'|(<r32  + <x42)/2Z)  d = trj^riyl/D 


and 


(M  = (<t4Fv/2v/Z),  a j lV\2y/D), 


y > 0 


= (<r3  WV  - 2 y/D,  aA  Vs!  - 2 y/D),  y < 0. 


(18) 


2 n=o  a — o n!  (n  + m)\  k\  (k  + m)! 

■ I*  dxe_0X(x  + 2)n+‘+mx‘+(1/2),"/m[cv/^] 
Jo 

= -e“t  y£n+£m(bl2sl2)2^m  (\d)2k+m 


2 n=o  it  = o r=  o n!  (n  m)\  k\  ( k + m )! 

. | (n  + k + mj  2„  + *+m-,. 

• p(Jxe"‘,xx‘  + (1,2'm+7m[c\^]. 
vo 


and  d,  and  rearranging  terms  for  clarity,  we  have  finally 


Pb(y)  = 


7 exP  i T7T  [2ffo2!>1  + V ‘ 


2D 


<V  + ff4 

x x a-  n + fc  + m / 

• I I I I 

m-0  n = 0 fc-0  r = 0 \. 


ct32  + (T42 


/gg3a.tx 
2D2(o32  + tr42) 


)■ 


. (ILY  ( °>wx2  v / d y 

\2D2/  \d(o}2  + (t42)'  \<t32  + <x42/ 


• £m  cos  m(v  — w — x) 


(k  + r)l 


n!  (n  + m)!  k!  (k  + w)! 


with 


f + ” + "V (— 

\ r / \2  D(g32  + a 4*’)/ 

(21) 


(f-9^o2)  = 


(<t4K,<t3H'\c42), 
(<r3  M/,<r4F,rr32), 


>’  > o 

y < 0. 


(22) 


By  writing  x = 1 4-  cosh  2u,  we  may  change  the  integral  to 

l-e“  ( dx  e~axIm[b\/$(x  + 2)]/m[cv/ixJ/m[d\/x(x  + 2)] 
2 J o 

= le0£  g (b!2sfl)2n+m  ({d)2k+n 


(19) 


Making  use  of  Gradshteyn  and  Ryzhik  [15,  formula 
6.643.4],  we  obtain  the  following  expression  for  (19): 

(£)"  i (“ + B 

£ * "+vw  (ib)2n+m(jd)2k+m 

»-o  k = o ,=  o n)  (n  + m)l  A!  (A  + m)! 


A recent  paper  by  Andrews  [5]  gives  an  expression  for 
the  density  function  in  nearly  ihe  same  situation  by  means 
of  the  characteristic  function  method.  Comparing  the 
development  here  with  that  of  Andrews,  it  appears  that 
following  the  strenuous  develooment  associated  with  the 
computations  of  the  characteristic  function,  one  still  needs 
to  solve  a convolution  integral  such  as  our  (19).  The  different 
final  expression  he  obtains  results  simply  from  a different 
method  of  attacking  this  type  of  integral.  For  the  special 
cases,  of  course,  the  form  of  the  expressions  given  here  is 
identical  to  that  of  Andrews.  Moreover,  the  computed 
results  we  show,  match  his  wherever  applicable. 

IV.  Reduction  of  the  General  Expression  to  Special 
Cases 

For  practical  applications  and  for  purposes  of  checking 
our  results  with  those  of  previous  authors,  it  is  instructive 
to  reduce  the  general  expression  of  (21)  under  certain 
assumptions. 

A,  Equal  Noise  Power  with  Even  Spectrum 

An  important  case  is  that  for  which  the  noise  inputs  are 
of  equal  power  and  their  spectra  are  even  about  the  center 
frequency  of  the  band.  Under  these  assumptions,  the  cross- 
quadrature correlation  coefficient  r is  zero  f 14]  and  a,2  = 
<t22  = a2,  implying  that  R = R'  = 0,  so  that  the  term  X 
in  (21)  is  zero.  Defining  the  input  SNR  for  channels  1 and 
2 by 


(20) 


where  the  C"+r(-)  are  the  Laguerre  polynomials  of  order 
k + r and  parameter  m.  Substituting  in  (17)  and  for  a,  b,  c, 


V 4 


V 

2cr2 


and 


I,  2 A 


?12 

2o2 


17< 
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Fig.  2.  Probability  density  functions  for  cross-correlator  output  ></) 
for  several  values  of  identical  input  SNR's  when  p — 0,  0,  = 02, 
and  <r2  = 1 . 


This  result  is  identical  to  that  of  Lezin  [7],  who  began 
analysis  with  these  assumptions. 

C.  Computed  Results  for  Special  Cases 

Judging  the  equal  power  (<r, 2 = o22  = a2),  even  spectrum 
case  to  be  of  sufficient  interest,  we  have  plotted  several 
curves  for  /i,  = h2  = h in  order  to  show'  the  effects  of  the 
parameters  p,  h2,  and  0 = (0,  - 02)/ 2;  a , quite  obviously, 
is  a scaling  parameter.  In  the  Appendix  the  pdf  (23)  is 
shown  to  be  equivalent  to 

P6()’)  = —2  exp  [-(/i32  + /u2)] 
a 


the  general  expression  of  (21)  is  rtducea  to  the  fo'low'ing 
form  for  _y  > 0: 


. y [Mj  ~ P)v/l"T~p1f  , (2 h,  j 2y  \ 
'=0  L sltyja2  J \ O V 1 + pi 


and  for  y < 0 


f i [id  - p)/.32]m  [id  + PM* 

„To  o ml  n ! 


exp 

[ -2  y 1 

G"  1 

V(l  + p). 

M !V(1  - p2)J 

exp 

. 

[ 2y 

G"  ~4V  1 

L<72(  1 “ P) J 

" L2(l  - p2)J 

(25) 


where  the  polynomials  we  have  defined,1 


c:  (X) 


if 


m + n - k 


)- 

/ kl 


Pe(y)  = exp 


2y 


1 a ( 1 - p) 


/ij2  = (/i,2  + hi2  +•  2/i,/i2  cos  20)  (23c) 

2(1  + p) 

= (/i  | 2 + h22  — 2/t|/i2  cos  20)  (23d) 

2(1  — p) 

and  where  20  = 0,  — 02  is  the  difference  in  phase  between 
the  narrow-band  signals  at  the  multiplier  input  (see  Fig.  1 
and  (3a)).  Muraka  [6]  presents  results  corresponding  to  a 
further  specialization  of  this  case  to  0 = 0 and  h2  = h22. 


B.  No  Input  Signals  and  Equal  Noise  Power 


If  the  correlator  inputs  are  assumed  to  be  only  the  noise 
of  equal  power,  that  is,  /j , 2 = h22  = 0 and  a 2 — a2  = 
ff2,  the  pdf  (21 ) is  further  reduced  to 


p*(y)  = { 


i 

- exp 
(T 


' 2y 

it!(1  + p). 

2y 

V(l  - p) 


y > 0 

y < 0.  (24) 


have  the  unusually  useful  computational  property  G"m  = 
G"m_,  + G"m~'.  The  form  (25)  results  from  application  of 
ti  e characteristic  function  method  to  the  special  case  we  are 
now  considering. 

For  computation  of  (25),  we  have  chosen  a “nominal 
case”  of  specific  values  for  the  parameters:  p = 0,  it  = 1, 
h — 1,  0 = 0.  In  Fig.  2,  the  parameter  h 2 (SNR  going  into 
the  multiplier)  is  varied.  It  is  seen  to  have  the  powerful 
effect  of  changing  the  pointed  curve  of  the  no-signal  case 
into  curves  which  begin  to  approach  the  familiar  Gaussian. 

The  correlation  coefficient  at  the  multiplier  inputs  p is 
varied  in  Fig.  3,  where  it  is  evident  thai  increasing  correla- 
tion produces  a more  “noiselike”  distribution,  that  is,  one 
with  a smaller  apparent  SNR.  Figs.  2 and  3 include  curves 
which  exactly  match  those  given  by  A.ndrews  [5],  although 
the  scaling  is  somewhat  different  due  to  his  inclusion  of  a 
scale  factor  in  the  definition  of  the  multiplier. 

The  phase  between  the  two  multiplier  signal  inputs  0 is 
varied  in  Fig.  4.  The  effect  of  this  parameter  is  seen  to  shift 
the  concentration  of  probable  values  closer  to  zero,  decreas- 
ing the  apparent  SNR.  Note  that  the  cases  0 = rr/6  and 
0 = Tij 3 are  “reflections”  of  one  another. 

V.  Remarks  on  the  Applications  of  the  Results 

The  pdf  for  the  output  of  a cross-correlator  with  bandpass 
inputs  arises  in  many  situations  of  practical  interest.  The 
equivalent  pdfs  (23)  and  (25)  have  been  applied  to  the 


1 A further  discussion  on  the  properties  of  this  polynomial  Gi(.v) 
is  given  in  the  Appendix. 
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Fig.  3.  Probability  density  functions  for  cross-correlator  output  \ii) 
for  several  values  of  noise  correlation  coefficients  when  h2  = 1, 
<rl  = 1.  and  l\  02. 


For  those  who  wish  to  pursue  the  mathematical  aspect 
of  this  type  of  problem,  consult  [23]— [30], 

Application  of  the  pdf  for  the  correlator  output  in 
problems  of  digital  communication  is  also  rewarding.  The 
correlator  model  we  have  considered  corresponds  to  the 
receiver  structure  for  binary  communication  in  which 
the  detection  is  accomplished  with  a noisy  reference  The 
performance  (error  probability)  computation,  when  the 
noises  between  the  input  and  the  reference  channels  are 
correlated,  can  be  computed  using  the  pdf  of  (23)  for 
/t/  = /i2 1 = h2.  Cooper  [2]  has  shown  that  'he  applica- 
tion of  the  pdf  for  the  correlator  outpi.t  to  binary 
orthogonal  communication  systems  produces  the  same 
known  result  of  error  probability  under  uncorrelated  noise 
(p  = 0)  assumptions. 


Acknowledgment 


0.48 


Fig.  4.  Probability  density  functions  for  cross-correlator  output  it r ) 

for  several  values  of  phase  difference  0 = (fl,  - ()-)  2 when  n = 0. 

It2  = I,  and  n2  = I . 

study  of  a class  of  multiplicative  array  configurations  [8] 
often  employed  in  sonar  or  underwater  detection  and 
estimation  applications.  In  this  application,  the  inputs  to 
the  multiplier  are  summations  of  the  outputs  of  array 
e'ements.  Having  an  explicit  formulation  for  the  pdf 
permitted  direct  calculation  of  the  effects  on  the  distribution 
of  parameters  such  as  signal  bearing,  the  number  of  array 
elements,  interelemcnt  correlation,  array  spacing,  and 
input  SNR 

Calculation  of  receiver  operating  characteristics  (ROC) 
based  on  our  results  revealed  that  the  multiplier  .-'filter 
combination  has  the  potential  for  performing  signal  detec- 
tion very  nearly  as  well  ax  the  more  conventional  square  - 
law  array  configuration,  whereas  previous  analysis  basevl 
on  SNR  alone  as  a figure  of  merit  predicted  a 3 dB  supe- 
riority for  the  conventional  system.  This  particular  applica- 
tion has  been  summarized  in  [22]. 


The  authors  wish  to  express  thanks  to  the  referees,  who 
suggested  the  relationship  of  the  polynomials  C/”,(.v)  to 
other  special  functions  and  who  drew  attention  to  [28]-[30], 


Appendix 

Another  Form  for  the  Probability  Density  Function 

The  case  we  are  considering  is  that  in  which  the  cross-correla- 
tion coefficient  r = 0 and  ra,  = ra2  (R  -*  0),  so  that  the  term  X 
in  (21)  equals  zero.  Under  these  conditions,  the  sum  and  dif- 
ference envelopes  Z,  and  Z2  arc  statistically  independent,  since 
we  may  factor  their  joint  probability  density  function  (14) 


/74(Z,./'-2) 


where  ra,:  - (rah  I - />),  ra4:  (rah  I />),//, : = S,:  2ra,:  - 

(// , ’ + hr  -v  2/1,/t,  cos  2d)  2( I 4 />),  and  /;4;  54;  2ra4* 

</ii’  4 lir  -h\h2  cos  2d)  2( I /»}.  /;(.\ ;2,<r ) is  the  non- 
central  chi-square  distribution  with  two  degrees  of  freedom  and 
noneenirality  parameter  a2.  Thus  the  (liter  output  v l(Z,:  - 
Z:;)  is  the  difference  between  two  scaled,  independent  non- 
central chi-square  random  variables.  This  type  of  relationship 
was  noted  also  by  1 ee  [16]  in  a similar  application. 

The  nonccntral  chi-square  distribution  has  many  forms  and 
interpretations  (see  [17]  [20])  for  example,  we  may  write  the 
follow  ing 


/*(.v;2 .<;’)  exp  i -i<r)  ^ J ( \ ; 2 -*  2 /•)  (27) 

r~* n r ! 


where  /*< a ;//►  is  the  (central)  chi-square  distribution  with  n 
degrees  of  freedom 


^ gjH 
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A.  Density  for  y Using  Characteristic  Function 

Using  an  expression  given  by  Jayachandran  and  Barr  [21  ],  we 
may  write  the  characteristic  function  of  the  distribution  of  y as 


also. 


m = «p  [-(V  + v>]  i i (A1" 

m = 0 n = o m i n\ 


z f g;t‘w=  z 

k = 0 \k ) k = 0 

- .1 


GJ  k(x ) 


icr/t)1 


2.\k  + l 


/I  - p\m+l  f (m+  n-  r\  [id  + p)Q 
\ 2 / rk\  m /(I  + /<x42/)r+1/ ' 


(29) 


Differentiation  formula: 
dxk 

Addition  formula: 


c:  (x)  = c:_k(x) 


Thus  we  obtain  for  the  density  of  >• 

p.<»  - 4exp  + V)j  z z 

m = 0 n -0  m ! 


G"Jx  + y)  = £ £ CJ-*<*). 

k = 0 ft  ! 


(34) 


(35) 


(36) 


[id  + ?)/>/]" 


y (m  + n - ks 

\/Lz  _£\'V 

4 y 

formulas: 

2m*  I 

k = o \ n 

K 2 ) 7 

<x2(l  4 p)’ 

i 

-1 

GnJx)  = — -m  - x) 

tn ! 

2 + 2A  , 

y > 0 

J 

+ 1 

4>(w  + |,  m + 

n + 2;  x) 

y lm  + n — r\ 

i (Lt/TV  1 

- -4  y . 

m! 

k\  m ) 

'1  2 ) * | 

o2(i  - pY 

Am  / 

= — : 2^0  -»».  » + >; 

* 1 

m!  \ 

2 4 2rj , 

y < 0.  - 130) 

J 

= V!^*|i2u)' 

f //I 

for  a = — 

Defining  the  polynomials 

m: 

m 

lm  + n - k\  xk 

Gnm  (x)  = Z 

k=0 

( n ) A! 

1 

=■  (-DmG;m'"’"1U). 

Relationships  to  Special  Functions:  The  polynomials  G'm(.v) 
may  be  related  to  the  confluent  hypergeometric  functions 
xV(aJ ? ;.v ) and  ;-l/.v),  to  the  Laguerre  polynomials 


(37) 

(38) 

(39) 


, h = 


we  finally  obtain 

/>„(>•)  = -iexp  [-(/,32  - k4*»  Z Z 

(7  n-'0m*0  W ! 

[id  + pivt 


exp 

- 2y 

Gi 

4.v 

02(l  + p) 

n2d  - p2) 

exp 

2>'  _ 

Gm 

I ~4>' 

o2(l  - p) 

XJn 

[rj2(l  - p2) 

y > (1 

v < 0, 

(31) 

B.  The  Polynomials  G"  (v) 

We  list  here  several  of  the  properties  of  these  polynomials. 
Definition: 

Gi(.v)  A £ + " ~~  k)  f G“(v)  = cJa)  CS  = I. 

v = o \ n / ft ! 


m - n 
2 

(40) 

(41) 

B>  using  the  identity  (40),  our  expression  (31)  for  the  pdf  can 
be  shown  to  be  identical  to  that  in  [30,  p.  63,  cq.  28]. 

C.  Equivalence  of  Forms 

We  now  show  that  the  form  we  have  just  derived  for  the 
probability  density  function  is  equivalent  to  that  of  (23)  under 
the  assumed  conditions 

1 < jil« 

I v lidw 

m=o  n — o m!  n: 
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_ y v1  / W t H — k \ X Cl  P 


Iterative  relationships: 

in  n 

Gj"(  v)  — - V g;  ,(  v)  = G,;  ,(.»■)  4 g:  ‘(.V) 

m:  r n 


(32) 
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using  Kummer's  transformation  on  the  last  step. 
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This  paper  shows  two  new  results  on 
binary  DPSK  communication  system.  The 
paper  is  based  on  a probability  density 
function  developed  under  this  contract. 


ON  THE  BINARY  DPSK  COMMUNICATION  SYSTEMS 
IN  CORRELATED  GAUSSIAN  NOISE 

Ohong  S.  Let  and  Leonard  E.  Miller 
ABSTRACT 

Using  the  recently  developed  probability  density  function  v/e 
obtain  "directly"  the  error  rate  expressions  for  the  binary  differential 
phase-shift  keyed  (DPSK)  systems  when  the  noise  values  at  the  sampling 
instants  in  adjacent  time  slots  are  statistically  dependent.  Tv/o  cases 
are  considered:  One  corresponding  to  equal  SNR  at  each  of  the  two  sampling 
instants,  and  the  other  to  unequal  SNRs.  The  consideration  of  the  former 
case,  together  with  the  assumption  of  unequal  a priori  symbol  probabilities 
Pq  and  P.j  results  in  an  error  rate  expression 

PE=  t {P0-P,)p]e'h2 

where  p is  the  noise  correlation  coefficient,  and  h2  is  the  SNR.  This 
expression  shows  clearly  why  P^  is  independent  of  noise  correlation  when  the 
source  symbols  are  equi -probable  provided  intersymbol  interference  is  assumed 
absent.  We  then  obtain  the  error  probability  expression  in  terms  of  unequal 
SNRs  (at  the  two  sampling  instants)  and  the  correlation  coefficient  p. 

Since  intersymbol  interference  in  a binary  DPSK  system  gives  rise  to  uneciual 
SNRs,  this  expression  provides  a useful  formula  for  estimating  the  system 
performance  under  such  circumstances. 

This  work  was  supported  in  Dart  by  the  Office  of  Naval  Research  under 
Contract  No.  N0001 4-A-0377-0021 . 


I. 


INTRODUCTION 


The  probability  of  error  in  binary  differential  phase-shift  keyed 
(DPSK)  systems  has  been  analyzed  by  several  authors  [lJ-[6].  Different 
methodologies  have  been  employed  by  the  previous  authors  in  obtaining 
the  error  rate  expression.  In  a broad  sense,  the  previous  techniques 
may  be  classed  into  two  categories: '‘indirect"  analogy  technique  and 
"characteristic  function"  method.  By  "indirect  analogy  technique"  we  mean 
that  the  computation  of  error  probability  is  based  on  the  knowledge  of 
other  communication  systems  such  as  binary  orthogonal  systems.  In  such 
an  approach  the  technique  does  not  require  the  probability  density  function 
(pdf)  of  the  decision  variable  because  the  probability  of  error  expression, 
namely  the  probability  distribution  function,  can' always  be  brouoht  to  a 
form  that  is  equivalent  to  that  of  binary  orthogonal  system.  Moreover, 
this  technique  has  always  been  restricted  to  the  case  where  the  noise 
samples  at  adjacent  time  slots  are  assumed  to  be  statistically  independent. 
The  characteristic  function  method  [6],  on  the  other  hand,  was  employer  to 
compute  th,_  error  rate  in  a most  general  situation  to  account  for  the 
correlated  noise  and  the  intersymbol  interference. 

Our  purpose  here  is  to  obtain  the  probability  of  error  expression 
in  a "direct"  fashion  by  using  the  pdf  of  the  decision  variable  developed 
recently  by  the  authors  [7].  We  obtain  the  error  rate  expressions  of  the 
two  cases  of  interest.  The  first  is  that  corresponding  to  equal  SNR  at 
each  of  the  two  sampling  instants,  and  the  other  to  unequal  SNRs.  The 
consideration  of  the  former  with  the  assumption  of  unequal  source  symbol 
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probabilities  will  lead  to  an  error  probability  expression  v/hich  shows 
clearly  how  it  is  related  to  the  correlation  between  two  noises. 

It  also  shows  why  the  error  probability  is  independent  of  noise  correlation 
when  the  symbol  probabilities  are  equiprobable  provided  that  the  inter- 
symbol interference  is  assumed  to  be  absent. 

Finally,  we  consider  the  case  of  unequal  SNRs  at  the  two  sampling 
instants.  Hubbard  [8]  has  shown  that  the  effects  of  intersymbol  inter- 
ference is  equivalent  to  consideration  of  different  SNRs  in  the  two  received 
pulses.  Thus,  the  error  rate  expression  we  obtain  in  terms  of  two  different 
SNRs,  along  with  the  noise  correlation  coefficient,  would  serve  a useful 
formula  for  estimating  the  "worst"  [8]  case  intersymbol  interference  perfor- 
mance. 


II.  ANALYSIS 

In  binary  DPSK  systems,  the  sequence  of  binary  source  symbols 
0 (space)  and  1 (mark)  differentially  phase-modulate  the  carrier  between, 
say,  0 and  n radians.  One  possible  rule  of  differential  phase  coding  is 
that  whenever  the  source  symbol  is  0 the  carrier  phase  remains  unchanged 
from  that  for  the  previous  symbol  (or  bit),  and  if  the  source  symbol  is  i, 
the  carrier  phase  is  shifted  ir  radians  from  that  for  the  previous  symbol . 
Thus,  a change  of  n radians  occurs  in  carrier  phase  whenever  the  symbol  1 
is  produced  by  the  source. 

The  receiver  model  for  the  binary  D°SK  system  to  be  considered 
is  shown  in  Figure  1,  which  is  an  ideal  product  demodulator.  As  shown  in 
the  figure,  the  two  inputs  to  the  multiplier,  u-j(t)  and  u9(t),  represent 
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the  received  waveforms  for  two  consecutive  source  symbols.  Each  input  is 
assumed  to  contain  bandpass  Information-carrying  signals  S^(t)  and  bandpass 
Gaussian  noises  n^(t),  i = 1,2: 

Uj(t)  = S,(t)  + n-j  (t)  (1) 

u2(t)  * S2(t)  + n2(t)  (2) 

Let  us  assume  that  the  "present"  and  the  "preceding"  source  symbols  are 
identified  with  carrier  phases  6-j  and  e2>  respectively.  Then  we  may  write 

S-j(t)  = A cos  (uit-0-| ) (3) 

and 

S2(t)  = A cos  Mt-T)-82]  (4) 

If  we  assume  the  carrier  frequency  to  have  been  chosen  such  that 
&J  = 2rk,  k integer, 

we  have 

S2(t)  = A cos  (ojt-02)  (5) 

The  noises  are  bandpass  stationary  Gaussian  processes  and  hence  we  may 
represent  them  in  the  form  of  well-known  Rice  decompositions  [9], [10]: 

n-j(t)  = x-j  ( t)  cos  «ut  + y-j  ( t)  sin  wt  (6) 

n2(t)  = x-j ( t-T ) cos  [u(t-T)]  + y-j ( t-T ) sin  [oj(t-T)] 

= x-j  (t-T)  cos  cut  + yj(t-T)  sin  u>t 

= x2(t)  cos  ojt  + y,(t)  sin  i>t  (7) 

t 
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What  Is  noted  in  (6)  and  (7)  is  that  the  noise  processes  n^(t)  and  n2(t) 
are  not  “two  different”  processes  but  rather  n2(t)  is  a translated  version 
of  the  noise  process  n-j(t). 

Our  objective  is  to  include  in  the  analysis  the  influence  of 
correlation  between  n-j(t)  and  n2(t).  Thus  we  have 

E[n1(t)n2(t)]  = E[n1(t)n1 (t-T)]  = a2p(T)  (8) 

and 

o’  * E[n|(t)]  = E[xf(t)]  ■ E[y|(t)];  10,2,  (9) 

and  p(T)  is  the  normalized  correlation  coefficient  of  noises  n-j(t)  and 
n2(t). 

The  pdf  of  the  Decision  Variable 

The  decision  variable  y ( t)  at  time  t is  the  ideally  lov/pass 
filtered  version  of  the  product 

x(t)  = u](t)xu2{t)  = [S1(t)+n1(t)]x[S2(t)+n2(t)] 

= [A  cos'  (ojt-e^ )+n-j (t)3x[A  cos  (<ut-e2)+n2(t)].  (10) 

We  now  need  to  know  the  pdf  of  y ( t ) at  the  zonal  lowpass  filter  output 
(Figure  1).  The  problem  of  obtaining  the  pdf  of  the  lowpass  filter  output 
for  the  system  model  that  fits  our  situation  was  solved  by  Miller  and  Lee  [7] 
in  greater  generality,  including  the  case  of  correlated  noises  of  unequal 
powers.  Our  need  here  is  a special  case  of  the  problem  treated  in  [7],  and 
from  (23)  in  [7]  we  have  the  pdf  expressions  as  follows: 
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For  y > 

p(y;p) 


0 with  noise  correlation  coefficient  p: 
7 exp  [-  -02(-]£p7  “ h3  ’ ^ hS] 


\ N^T  j y 

r=0  L TW^2  -I  r ' 0 ^1+p'  r 
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and  for  y < 0 wifh  noise  correlation  coefficient  p: 

P(y;p)  - ^ exp  - (^)  h|  - h|] 


where 


h5  “ irkr  thf+hl+?hi  h2 cos  (0r®2^ 

hi  ^ 2 ( 1 -p ) thi+h2’2hlh2  cos  (°re2^ 


(lib) 


(12a) 


(12b) 


Ir(  ) is  the  modified  Bessel  function  of  the  first  kind  of 
order  r 

Lr(  ) is  the  Laguerre  polynomial  of  order  r 
and  h?  is  the  signal-to-noise  ratio  (SNR)  in  channel  i,  i=l,2. 


It  is  evident  that  the  pdf  given  in  (11)  is  that  of  a general 
form  in  that  the  two  inputs  to  the  multiplier  (Figure  1)  have  different 
signal  power  levels  and  the  noises  are  correlated.  In  another  sense,  how- 
ever, the  pdf  given  in  (11)  is  a "restricted"  form  in  that  the  noise  inputs 
are  of  equal  power  and  their  spectra  are  even  about  the  center  frequency  of 
the  band  (symmetrical  bandpass,  see  Section  IV  of  [7]).  With  the  above  pdf 
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we  will  be  able  to  obtain  the  error  probabilities  of  binary  DPSK  system 
under  two'  different  assumptions  of  practical  interests.  Namely,  the  error 
rates  for  equal  and  unequal  input  power  level,  respectively. 

It  has  been  known  [5],  [6]  that  when  the  source  symbols  are  equi- 
probable  and  the  intersymbol  interference  is  assumed  absent,  the  error  prob- 
ability of  a binary  DPSK  system  is  independent  of  nois^  correlation.  How- 
ever, when  the  prior  probabilities  of  the  source  sysmbols  are  unequal,  the 
error  probability  is  found  to  depend  very  much  on  the  noise  correlation.  The 
absence  of  intersymbol  interference  corresponds  to  a situation  where  the 
two  input  powers  to  the  multiplier  a^e  equal.  It  may  be  of  interest  to  know, 
then,  in  what  manner  the  error  probability  is  affected  by  the  noise  correlation 
In  fact,  in  a laboratory  simulation  for  DPSK  system,  the  binary  symbols 
are  not  always  generated  in  an  equi -probable  fashion,  and  the  simulation 
results  are  not  always  "poorer"  than  the  theoretical  prediction  of  the 
ideal  situation.  As  will  be  noted  subsequently,  the  simulation  results 
that  are  "better”  than  the  ideal  theoretical  predict ion^-^e  readily 
reconcilable  in  view  of  the  error  rate  expressions  to  be  derived  in  tne 
following  section. 

Also  of  interest  is  the  error  probability  expression  for  the 
binary  DPSK  system  for  unequal  input  rower  levels.  A practical  situation 
which  corresponds  to  this  assumption  is  that  of  intersymbe!  interference, 
since  the  attribute  of  intersymbol  interference  is  the  unequal  pulse  levels 
fer  the  two  successive  symbols.  It  will  be  explicit  in  the  expression  to 
be  derived  later  that  in  such  case  the  error  performance  depends  upon  the 
noise  correlation  regardless  of  the  distribution  of  the  source  symbols. 
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III.  THE  PROBABILITY  OF  ERROR  WHEN  THE  SUCCESSIVE  PULSES  HAVE  EQUAL  POWER 
The  pdf  that  Is  applicable  to  the  particular  situation  where  the 
two  input  channels  to  the  multiplier  (Figure  1)  have  equal  SNRs  is  obtained 
when  we  substitute  into  (11)  the  conditions 


h^  = h2  = h2  = A2/ 2a2. 
We  then  obtain: 
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In  binary  DPSK  system,  the  binary  decision  is  based  on  the  compari- 
son of  the  decision  variable  with  "zero"  threshold.  The  reason  is  that 
when  noises  are  assumed  to  be  absent  at  the  input,  the  decision  variable 
is  given  by 

y(t)  \t  cos  (0r0z)  os) 

Since  e-j^  either  0 or  +*  in  DPSK  system,  only  the  sign  of  y(t)  is 
significant  in  making  a binary  decision. 

30^ 
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Let  us  assume  that  the  reference  (previous)  phase  Og  is  0. 

Then  8^  is  either  0 or  it  radians,  depending  upon  whether  the  "next"  symbol 
(present  symbol)  is  0 (space)  or  1 (mark),  so  that  e^-e.,  = 0 or  e-j-Gg  =ir- 
If  f^is  it  instea  then  e i -Gg  = 0 or  e-j-eg  = -*»  again  depending  on  whether 
"zero"  or  "one"  i r.  transmitted, respectively.  From  (14b)  it  can  be  easily 
verified  that 

p(y;p|e1-e2  * *)  * p(y;p |eT-e2  = -*)  (16) 

Thus,  the  conditional  pdf's  are  obtained  as  follows: 

p(y;P|space)  = p(y  ;P  j e-, -e2  = 0)  (17a) 

p(y;P|mark)  = p(y;p  |e-,-e2  = ir)  (17b) 

The  density  functions  of  (17)  with  the  normalization  of  u - y/a2  are  given 
by 


p(u;P|space)  = 


and 


exP(y~)  exp(-h2)  ; u<0 


^ ^,exp(-  w > j0  ^ 


u<0 


(18a) 


p(u*rf)|mark)  = 


“O*1  £>1  [^]r‘rC 


r=0 


2u 


( exp{-  jj^-)exp(-hz)  ; u>0 


7- 
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(18b) 


From  (18a)  and  (18b)  we  observe  the  "symmetry"  properties  to  be 

p(u;p | space)  = p(-u;-p|mark)  (19) 

The  conditional  pdf  (19)  is  plotted  under  various  parameter  conditions  in 
Figures  2-3.  In  Figure  2 the  variations  of  the  conditional  pdf  for  different 
correlation  coefficients  are  shown  for  fixed  identical  input  SNRs  of  the  two 
multiplier  inputs.  Note  the  dramatic  influences  of  the  correlation  coeffici- 
ents for  each  input  SNR.  In  Figure  3 we  have  shown  the  variations  of  the  pdf 
for  different  SNRs  for  fixed  noise  correlation  coefficients. 

To  compute  the  probability  of  error  we  first  obtain  the  conditional 
error  probabilities: 

P(e | space ) = probTy<0 ! ©i ~©2  = 0}  (20) 

and 

P(e|mark)  * prOb^^>0 1 el ~e2  = ±1t}  (21 ) 

To  compute  (20),  we  need  only  obtain  from  (14) 
p(y;p|e1-02  s o,  y<0)  = p'(y;p|e1-02  = 0) 
which  can  be  shown  to  be 

n (y;p|i>re2 " 01  * T2  exp  ^T^p)  ' 2hZ(T^  H h2(fe,] 

■ * ^ exp  'h2]  (22) 

In  a similar  manner,  to  compute  (21)  we  need  to  obtain  from  (14) 


P+(y.p)8,-82  ■ »)  ■ Jr  exp  [-  - 2h*  {^)  ♦ h*  (j^)] 


Using  (22)  and  (23),  we  compute  the  conditional  probabilities  of  error, 
and  we  get 

P(e|space)  = J°  p_(y;P|are2  = 0)  dy  = ^ e'h2  (24) 

and 

P(e|mark)  .f  p+(y;p|ei-e2  = *)  dy  --  ^ e-h2  (25) 

Assunrng  that  the  a priori  probabilities  of  source  symbol  0 (space)  and 
1 (mark)  to  be  Pg  and  P^,  respectively,  we  obtain  the  unconditional  proba- 
bility of  error  to  be 

PE  = Po^eisPace)  + P-]  P (e  jmark ) 

- ? P + (P,-P0)p]  e~^2  (26) 

When  the  source  rymbol  probabilities  are  equiprobable,  namely  when  Pg  ~ P1 » 
(26)  reduces  to  the  well  known  classical  expression  P£  = j exp(~h2)  for  the 
binary  DPSK  error  rate,  where  h2  is  the  input  carrier-to-noise  power  ratio 
(SNR). 


From  the  above  expressions,  we  observe  the  known  fact  [5], [6] 
that  the  error  probability  is  independent  of  noise  correlation  when  the 
source  symbols  are  equi-probable.  It  is  interesting  to  note  the  simplistic 
manner  in  which  these  variables  affect  the  ideal  classical  error  rate 
exp-ession. 

IV.  THE  PROBABILITY  OF  ERROR  WHEN  THE  SUCCESSIVE  PULSES  HAVF 

UNEQUAL  POWER 

Our  purpose  here  is  to  compute  the  probability  of  error  in  the 
binary  DPSK  system  when  the  two  multiplier  inputs  have  different  SNRs.  As 


- //  - 


alluded  to  in  the  earlier  discussions,  this  case  corresponds  to  that  of 
intersymbol  interference.  Our  aim  is  to  show  that,  under  such  conditions, 
the  error  rate  expression  is  not  independent  of  the  noise  correlation  p(T) 
even  if  the  binary  source  symbols  are  assumed  to  be  equi-probable. 

Here  we  make  use  of  the  general  form  of  the  pdf  given  in  (11). 
The  method  of  analysis  parallels  that  considered  in  Section  III. 


Let  us  assume  that  the  reference  phase  is  0 and  the  transmitted 
source  symbol  is  1 (mark).  Then  e-j -02  = it.  The  conditional  probability  of 
error  is  thus  computed  from 


P(e;p,h1,h2|mark)  = J"  p+(y;p,h1  .h2)e1-e2  - n)dy  (27) 


where  the  pdf  in  the  integrand  is  given  by  (11a)  with  -e2  = it.  The  inte- 
gration indicated  by  (27)  then  results  in 


P(e;p,h^ ,h2| mark ) 


= ^-  exp  - 


(h1+h2)2  (hrh2)s 
4 2 ( 1 +p ) 


X 

r=0 


ihel 

Ji 


(h-,-h2) 


• /87^ 


[' T (hl+h2>2]  Fr<»sVh2>  <2f» 


where 


Fr(p;h1,h2)  4 


-«yy-r/2 


Ir(2e.y)dy 


and  where 


„ A 2 
o^TT+pT 

& hrh2 

8 = ^y 


(29a) 


(29b) 

(29c) 


•34^ 


Using  the  formulas  9.220.2  and  6.643.2  in  [11],  one  can  show  that  (29a) 
reduces  to 


Fr(p»hi ,h2)  = p-  0r  a ^ -jF^ (1  »r+l ;e2/2a) 


(30) 


where  ) Is  the  confluent  hypergeometric  function. 


Putting  (29b)-(30)  into  (28)  we  obtain 
P(e;p,hi,h2|mark)  = exp 


[- 


(h1+h2)2 


rh2)zl 

TwrJ 


(h 
2(l+p) 


X k 

r=0 


r[t  (i^)(hrh2)2]  iFi  [1>r+1i  4(i+p j ] 


x L. 


[-  1 <^<Vh2>2] 


(31) 


It  can  be  shown  that  the  other  conditional  probability  of  error  is  exactly 
the  same  as  (31)  except  the  sign  change  in  p: 


P(e;p»h-j,h2|space)  = P(e;-p  ,h1  ,h2|mark) 


(32) 


When  h^hg  is  imposed  on  (31)  and  (32),  the  results  reduce  to  the  expression 
we  have  already  obtained  in  Section  III  [see  (24)  and  (25)].  That  is. 


when  h^  = h|  = h2 


and 


P(e;p,h1,h2|mark)  = e"h2 


P(e;P,h1,h2|space)  = e_h2 


(33a) 


(33b) 


35< 


The  probability  of  error  for  the  case  of  unequal  input  SNRs  for 
equi- probable  binary  DPSK  system  is  then  given  by 

PE  “ 1 CP(e;p,h-, ,h2 jmark)  + P(e;p,h1  ,h2|space)] 

B ^ [P(e;p,hph2  |mark)  + P(e;-p»h1  ,h2jmark)]  (34) 

where  it  is  noted  that  the  conditional  error  probability  (31)  is  all  that 
Is  needed  to  compute  (34). 
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Figure  1:  The  receiver  block  diagram  for  binary  DP5K 

communication  system. 
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Fiaure  2a:  pdf  of  the  normalized  (u =y/a2)  decision  variable 

(cross-correlator  output)  with  noise  correlation 
coefficients  as  parameters  when  h2  = SNR  = 0.5. 
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Figure  2b: 


pdf  of  the  normalized  (u -y/a2)  decision  variable 
(cross-correlator  output)  with  noise  correction 
coefficients  as  parameters  when  h2  = SNR  = 2. 


p(u;/>!  space)  = p(-u;-/o|  mark ) 


Figure  2c:  pdf  of  the  normalized  (u=y/o2)  decision  variable 

(cross-correlator  output)  with  noise  correlation 
coefficients  as  parameters  when  h2  = SNR  = 4. 


Figure  3a: 


pdf  of  the  normalized  (u=y/CT2)  decision  variable 
(cross-correlator  output)  with  h2  = SNR  as 
parameters  for  p = 0 (a  "classical"  case). 


p(u;pl  space)  = p(-u;-plmark) 
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Figure  3c:  pdf  of  the  normalized  (u=y/o2)  decision  variable 

(cross-correlator  output)  with  h2  e SNR  as 
parameters  for  p = 0.4. 
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